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1.1 Submicron Structures 
Laboratory 

The Submicron Structures Laboratory at MIT 
develops techniques for fabricating surface 
structures with linewidths in the range from 
nanometers to micrometers, and uses these 
structures in a variety of research projects. 
These projects of the laboratory, which are 
described briefly below, fall into four major 
categories: development of submicron and 
nanometer fabrication technology; nanometer 
and quantum-effect electronics; crystalline 
films on non-lattice-matching substrates; and 
periodic structures for x-ray optics, spectros¬ 
copy and atomic interferometry. 


1.2 Microfabrication at 
Linewidths of 100 nm and 
below 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation 
Grant ECS 87-09806 

Project Staff 

Martin Burkhardt, James M. Carter, William Chu, 
Kathleen Early, Reza A. Ghanbari, Yao-Ching Ku, 
Alberto Moel, Hui Meng Quek, Or. Mark L. 
Schattenburg, Professor Henry I. Smith, Siang- 
Chun The, Anthony Yen 


1 Princeton University, Princeton, New Jersey. 
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A variety of techniques for fabricating struc¬ 
tures with characteristic dimensions of 0.1 
/tm (100 nm) and below are investigated. 
These include: x-ray nanolithography, holo¬ 
graphic lithography, achromatic holographic 
lithography, electron-beam lithography, 
focused-ion-beam lithography, reactive-ion 
etching, electroplating, and liftoff. Develop¬ 
ment of such techniques is essential if we are 
to explore the rich field of research applica¬ 
tions in the deep-submicron and nanometer 
domains. 

X-ray nanolithography is of special interest 
because it can provide high throughput and 
broad process latitude at linewidths of 100 
nm and below, something that cannot be 
achieved with scanning-electron-beam or 
focused-ion-beam lithography alone. We are 
developing a new generation of x-ray masks 
made from inorganic membranes (Si, Si 3 N 4 , 
and SiC) in order to eliminate pattern dis¬ 
tortion. To achieve multiple-mask alignment 
compatible with 50 nm linewidths, we will 
fix the mask-sample gap at 4 /tm, translate 
the mask piezoelectrically, and detect align¬ 
ment to < 10 nm by a dual-grating interfero¬ 
metric scheme. Phase-shifting x-ray masks 
should permit us to achieve 50 nm linewidths 
at such gaps. In previous studies we showed 
that a pi-phase-shifting mask improves 
process latitude by increasing the irradiance 
slope at feature edges. For linewidths below 
50 nm we bring the mask membrane into soft 
contact with the substrate by electrostatic 
means. A variety of techniques are used to 
pattern the x-ray masks including e-beam 
lithography, focused-ion-beam lithography 
(FIBL), holographic lithography and sidewall 
shadowing. 

Since the early research on x-ray lithography, 
it was believed that photoelectrons, released 
when x-ray photons are absorbed, severely 
limit the resolution of x-ray lithography for 
linewidths below 100 nm. For example, it 
was estimated that using the Al* x-ray (A = 
0.834 nm, E = 1.5 keV), it would not be 
possible to reliably achieve linewidths below 
100 nm. We explored this issue exper¬ 
imentally and found previous estimates to be 
in error. Figure 1 shows that a 30 nm wide 
line on an x-ray mask is faithfully replicated 
with no measurable linewidth change using 
any one of three x-ray wavelengths, Al* 



Figure 1. Replication in PMMA of a 30 nm-wide gold 
absorber line with (a) C* (A = 4.5 nm), (b) 
CulU = 1.3 nm), and (c) AI*(A = 0.83 nm). 


(0.834 nm), Cu L (1.33 nm), or C* (45 nm). 
The implication is that for linewidths down 
to ~ lOnm, diffraction and not photoelectron 
range is the factor that will limit resolution. 

A new type of achromatic holographic litho¬ 
graphy (AHL) was developed that enables us 
to achieve 100 nm-period gratings (50 nm 
nominal linewidth). Previously, 200 nm was 
the minimum period achievable. The AHL 
scheme is illustrated in figure 2, and the 
exposure results shown in figure 3. This 
technology will be used to make gratings for 
x-ray spectroscopy and atom beam 
interferometry, and to fabricate new classes 
of quantum-effect electronic devices. 



Figure 2. Schematic of the achromatic holographic 
configuration. 
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Figure 3. Scanning electron micrograph of 100 nm- 
period grating lines after exposure, shadow evaporation 
of nickel and reactive ion etching in an oxygen plasma. 


1.3 Improved Mask Technology 
for X-Ray Lithography 

Sponsors 

Semiconductor Research Corporation 
Contract 87-SP-080 
Hampshire Instruments Corporation 

Project Staff 

Yao-Ching Ku, Kenneth Lu, Professor Henry I. 
Smith, Lisa Su, JoAnne Gutierrez, Flora S. Tsai 

In order to utilize x-ray lithography in the 
fabrication of submicron integrated elec¬ 
tronics, distortion in the x-ray mask must be 
eliminated. Distortion can arise from stress 
in the absorber, which is usually gold or 
tungsten. Tungsten is preferred because it is 
a closer match in thermal expansion to Si, 
SiC and other materials used as mask mem¬ 
branes. However, W is usually under high 
stress when deposited by evaporation or 
sputtering. Earlier, we demonstrated that for 


a given type of substrate, zero stress (i.e., 
less than 5 x 10 7 dynes/cm 2 ) can be achieved 
by controlling the sputtering pressure to 
within one-tenth of a militorr. We are now 
developing a computer-controlled system for 
monitoring in situ, during deposition, the 
stress in sputtered W on x-ray mask mem¬ 
branes. Stress is determined from the reso¬ 
nant frequency of the membrane. 

We have also investigated a number of mate¬ 
rials as mask membranes including: Si, Si 3 N 4 , 
SiC, and laminates of Si 3 N 4 /poly Si/Si 3 N 4 
and Si0 2 /Si 3 N 4 /Si02. The strongest mem¬ 
branes were Si rich Si 3 N 4 A 1.2 fzm thick 
membrane of this material can sustain a full 
atmosphere pressure differential across a 
span of 20 mm. There is a concern, 
however, that Si 3 N 4 will undergo stress 
changes during the heavy x-ray irradiation 
anticipated in 1C production. We are investi¬ 
gating this issue as well as novel approaches 
to enhance the durability of Si membranes. 


1.4 Study of Electron Transport 
in Si MOSFETs with 
Deep-Submicron Channel 
Lengths 

Sponsor. 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 

Project Staff 

Professor Dimitri A. Antoniadis, Gregory A. Carlin, 
Paul G. Meyer, Ghavam Shahidi, Professor Henry 
I. Smith, Siang-Chun The 

Electronic conduction in sub-100 nm 
channel length Si MOSFETS is being studied 
in order to observe non-stationary transport 
effects, i.e., effects arising from rapid spatial 
variation of the electric field, and the 
resulting disparity between the actual carrier 
temperature and the temperature that would 
correspond with the given field at steady- 
state conditions. In simple terms, non- 
stationary effects are observed when transit 
times across the very short channels become 
comparable to, or shorter than, energy relaxa¬ 
tion times. Previous work has produced 
NMOS devices with channel lengths down 
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to 60 nm and gate oxides as thin as 2.5 nm. 
The devices were fabricated with a combina¬ 
tion of x-ray nanolithography and optical 
projection lithography. Electron velocity 
overshoot and a reduction of hot-electron 
effects, as measured from the substrate 
current, were observed in these devices. 

Our focus has now shifted to the fabrication 
of self-aligned deep-submicron MOSFETs. 
These devices will use polysilicon gates and 
cobalt salicidation (self-aligned silicidation). 
The aim is to reduce parasitic source and 
drain resistances and gate-source overlaps. 
This will permit a clearer observation of the 
non-stationary effects, e.g., in both substrate 
and gate currents, and will enable us to push 
device switching speed. 

Deep-submicron self-aligned silicided NMOS 
devices have been successfully fabricated 
using optical projection lithography followed 
by partial photoresist ashing to define small 
features. Devices with effective channel 
lengths down to 230 nm were fabricated 
using this technique. As in the non-self- 
aligm'd device process, both boron and 
indium were used as implant species to 
control punchthrough and threshold voltage. 
The sheet resistance of the source-drain 
regions was reduced by silicidation from 110 
to 13.6 Ohms per square, and source-drain 
contact resistances on the devices were 
reduced from 52 to 10.7 Ohms. 

The self-aligned NMOS process is under¬ 
going further development to include x-ray 
nanolithography for gate polysilicon pat¬ 
terning. This has included the development 
of inorganic x-ray mask technology to avoid 
distortion and achieve precise alignment. A 
simple optical scheme to align x-ray masks to 
substrates to less than 1 has also been 
developed. The absorber on the x-ray mask 
will be patterned using a combination of 
optical lithography and focused-ion-beam 
(FIB) lithography. In addition, an effort is 
being made to further improve consistency 
and control in the salicidation process. 

Work is also in progress to develop a corre¬ 
sponding deep-submicron self-aligned 
PMOS process. This should give us 100 
nm-channel-length CMOS circuits fabricated 
by a technology compatible with commercial 
mass production. 


1.5 Studies of Electronic 
Conduction in One-Dimensional 
Semiconductor Devices 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation 
Grant ECS-85-03443 

Project Staff 

Professor Dimitri A. Antonaidis, Stuart B. Field, 
Professor Marc A. Kastner, Udi Meirav, Samuel L. 
Park, John H.F. Scott-Thomas, Professor Henry I. 
Smith 

Sophisticated processing techniques and 
advanced lithography have allowed us to 
enter what we believe is a fundamentally 
new regime in the study of electronic con¬ 
duction in one-dimensional systems. A 
slotted-gate MOSFET structure (figure 4) 
was used to produce an electron gas at the 
Si/Si0 2 interface beneath the gap in the 
lower-gate. This was done by biasing the 
upper gate positively, while keeping the 
slotted gate just below threshold. Fringing 
fields around the lower gate confined the 
electron gas to a width substantially nar¬ 
rower (~ 25 nm) than the distance sepa¬ 
rating the two halves of the slotted gate (~ 
70 nm). The slotted gate was produced 
using x-ray nanolithography and liftoff. It 
was composed of refractory metals to allow a 
subsequent high temperature anneal. This 
anneal removed damage created by the 
e-beam evaporation of the refractory metal, 
so that the electron gas had a mobility of 
15,000 cm 2 / V-sec at 4.2K. The electrical 
conductance of the 1-D gas was measured 
as a function of the upper gate voltage for 
temperatures less than 1K, and a surprising 
series of periodic oscillations was seen in the 
conductance (figure 5). Changing the gate 
voltage can be thought of as changing the 
number of electrons per unit length of elec¬ 
tron gas. Since the conductance is thermally 
activated, the oscillations reflect a periodic 
change in the activation energy of the elec¬ 
tron gas as the electron density is changed. 
Computer simulations solving Poission's 
equation and the single particle Schrodinger 
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Figure 4. (a) Schematic cross section and (b) top 

view of the slotted-gate device. The inversion layer, 
formed by the positively biased upper gate is confined 
by the lower gate. The thermal oxide and refractory 
metal lower gate are both 30 nm thick, and the chem¬ 
ical vapor deposition oxide is 45 nm thick. The width 
of a narrow inversion layer is exaggerated in (b). 

wave equation strongly suggest that the 
electron gas is dynamically one-dimensional 
when the oscillations are most strongly seen. 
That is, the electrons are in the lowest 
quantum energy level of the potential well 
created by the fringing fields of the slotted 
gate. 


A one-dimensional electron gas will form a 
charge density wave (CDW) when the 
electron-electron interaction is the dominant 
energy. This is the same regime that the 
device is in when the oscillations are seen. 
In a CDW, electrons are spaced periodically. 
This minimizes the total energy of the elec¬ 
tron gas. Impurities in the device can pin the 
CDW. Two impurities in the conducting 
channel would in turn strongly pin and 
weakly pin the CDW as the electron density 
of the gas is changed. This corresponds to 
commensurability and incommensurability of 
the CDW with the spacing of the impurities, 
and results in oscillatory behavior of the con¬ 
ductance with electron density. The period 
of the oscillations varies randomly from 
device to device, with no dependence on 
length. The period changes in the same 
device when it is heated to room temperature 
and then cooled. This strongly suggests that 
mobile impurity atoms are doing the pinning. 
It should be stressed that these oscillations 
are seen with no magnetic field, implying 



8.40 8.45 8.50 8.55 


V c (Volts) 



Figure 5. Top panel: Conductance G versus gate 
voltage Vg for a 2 /<m-long inversion layer. Bottom 
panel: Fourier power spectrum of the top panel data. 
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that the phenomenon is fundamentally dif¬ 
ferent from phenomena requiring Landau 
quantization (such as the Quantum Hall 
Effect). Also, the effect seems to require a 
many-body theory (i.e., one incorporating 
electron-electron interactions). 

1.6 Surface Superlattice 
Formation in Silicon Inversion 
Layers Using 0.2 /im-Period 
Grating-Gate Field-Effect 
Transistors 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
U.S. Air Force - Office of Scientific Resea:ch 
Grant AFOSR-88-0304 

Project Staff 

Professor Dimitri A. Antoniadis, Phillip F. Bagwell, 
Professor Marc A. Kastner, Professor Terry P. 
Orlando, Professor Henry I. Smith, Anthony Yen 

We have been studying distinctly quantum 
mechanical effects in electrical conduction 
using the silicon grating-gate field-effect 
transistor (GGFET). The Si GGFET is a dual 
stacked-gate MOS-type structure in which 
the gate closest to the inversion layer 
(bottom gate) is a 200 nm period grating 
made of refractory metal. A Si0 2 insulating 
layer separates the grating gate and the 
inversion layer from a second continuous 
aluminum gate (top gate). Two types of 
electronic devices are possible in this geom¬ 
etry: (1) a lateral surface superlattice (LSSL) 
GGFET «n which carrier conduction is per¬ 
pendicular to the grating lines; and (2) a 
quasi-one-dimensional (QID) GGFET which 
confines carriers to move in narrow inversion 
strips beneath the grating lines. 

Here we describe our observation of a transi¬ 
tion from conduction through ~ 300 multiple 
parallel quantum wires to conduction in a 2D 
electron gas (standard MOSFET) using the 
QID GGFET. We can induce such a dimen¬ 
sional transition by varying the two gate 
voltages of the GGFET. When a strong mag¬ 
netic field is applied perpendicular to the 
plane of the inversion layer at temperatures 
below 4.2K, we observe a large negative 


transconductance whose onset voltage marks 
the transition from 1 D to 2D conduction. 

This is shown in figure 6, where the transi¬ 
tion from 1 D to 2D behavior occurs at a top- 
gate voltage around zero volts. The 
significant observation is that with a mag¬ 
netic field perpendicular to the plane of the 
channel, adding more electrons to the 
MOSFET actually reduces the device con¬ 
ductance at the ID to 2D transition point. 
The oscillations below the transition point 
are due to magnetoelectric subbands, while 
those above the transition are due to 
Shubnikov de Haas oscillations as the 
Landau levels pass through the Fermi level. 
This striking dimensional transition in a mag¬ 
netic field has never been observed before. 


Family of Magnetic Fields (Tesla) 



Figure 6. Plot of the conductance of a quasi-ID 
grating-gate MOSFET as a function the top gate 
voltage, for several values of the perpendicular mag¬ 
netic field. At a top gate voltage around zero volts 
there is a transition from 1 D to 2D conduction with an 
attendant strong negative transconductance. 

The origin of the unusual magnetoresistance 
effects that we have observed is still unclear. 
At the magnetic fields where the effect 
becomes prominent, the cyclotron radius is 
much smaller than the width of the quantum 
wires, so a semiclassical explanation 
involving an interplay of the cyclotron radius 
and wire width is probably inappropriate. 
Also, the two terminal magnetoresistance 
measurements significantly complicate inter¬ 
pretation of the data. A noninvasive four- 
terminal measurement would be ideal, but it 
is not clear that such a measurement is pos- 
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sible. Finally, the exact way in which the 
current feeds into the source contact from 
the bulk device may also be important. At 
higher magnetic fields where the magneto¬ 
conductance oscillations begin to appear, 
conduction takes place through edge states. 
These edge states may also play a role in 
giving rise to the effect. 


1.7 Study of Surface 
Superlattice Formation in 
GaAs/GaAIAs Modulation 
Doped Field-Effect Transistors 

Sponsor 

U.S. Air Force - Office of Scientific Research 
Grant AFOSR-88-0304 

Project Staff 

Professor Dimitri A. Antoniadis, Martin Burkhardt, 
William Chu, Professor Jesus A. del Alamo, Reza 
A. Ghanbari, Khalid Ismail, Professor Marc A. 
Kastner, Professor Terry P. Orlando, Professor 
Henry I. Smith, Anthony Yen 

We have used the modulation-doped field- 
effect transistor (MODFET) as a test vehicle 
for studying quantum effects such as electron 
back diffraction in a GaAs/AIGaAs material 
system. In a conventional MODFET the 
current transport is modulated by a contin¬ 
uous gate between source and drain. In our 
studies, we have used Schottky metal 
gratings and grids for the gate, as illustrated 
in figure 7. Such gates produce a periodic 
potential modulation in the channel. 

The grid was produced by x-ray nanolitho¬ 
graphy and liftoff. The x-ray mask of the 
grid was produced by two successive x-ray 
exposures, at 90 degrees to one another, 
using a master mask that was fabricated via 
holographic lithography. The latter yields 
coherent gratings over areas several centime¬ 
ters in diameter. A new technique was 
developed that yields grating and grid pat¬ 
terns only in the channel region between 
source and drain. This has simplified the 
overall process and enhanced its reliability. 

The MODFET is normally on, that is, a nega¬ 
tive gate bias of about — 0.2 V must be 
applied to pinch off conductance from 
source to drain. As the gate bias is used 


G*AS/G»A1A» MODFET Latercl Surface 
Superlattice (botb fraticf and grid) 
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Figure 7. Schematic cross section of a grid-gate 
MODFET device. Contacts to the grid are made by 
pads off to the sides of the conduction channel. 


above this threshold point, the height of the 
periodic potential modulation is reduced and, 
simultaneously, the Fermi energy is raised 
(or, equivalently, the electron wavelength is 
reduced) in the 2D electron gas residing at 
the AIGaAs/GaAs interface. When the elec¬ 
tron wavelength phase-matches the periodic 
potential, electron back-diffraction occurs 
provided the inelastic length (i.e., the coher¬ 
ence or phase breaking length) is longer than 
the grating-period. Such back diffraction is 
manifested by a drop in the conductance. A 
stronger back diffraction effect is observed in 
the case of a grid because true minigaps are 
formed. The measurements of conductance 
modulation of grating and grid-gate 
MODFETs agrees with theoretical pred¬ 
ictions. In the grid gate devices it was also 
possible to observe negative differential 
resistance which might be due to sequential 
resonant tunneling. 

We plan to decrease the periodicity of the 
gratings and grids by a factor of two, to 100 
nm period. For devices with such fine grating 
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periodicity the superlattice effect should 
become more pronounced and observable at 
higher temperatures. We will also conduct 
magnetotransport measurements with devices 
of 100 and 200 nm periodicity. 

1.8 Study of One-Dimensional 
Subbands and Mobility 
Modulation in GaAs/AIGaAs 
Quantum Wires 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
U S. Air Force - Office of Scientific Research 
Grant AFOSR-88-0304 

Project Staff 

Professor Dimitri A. Antoniadis, Phillip F. Bagwell, 
William Chu, Professor Jesus A. del Alamo, Reza 
A. Ghanbari, Khalid Ismail, Professor Marc A. 
Kastner, Professor Terry P. Orlando, Professor 
Henry I. Smith, Anthony Yen 

In order to study one-dimensional conducti¬ 
vity in the AIGaAs/GaAs modulation-doped 
structure without the conductance fluctu¬ 
ations normally associated with single micro¬ 
scopic systems, we fabricated arrays of 100 
parallel quantum wires (MPQW). The meas¬ 
ured Hall mobility in the 2-D electron gas at 
the AIGaAs/GaAs interface was 250,000 
cm 2 /Vs at 4.2K. The multiple quantum wires 
were produced by x-ray nanolithography, 
liftoff of Ti/Au, and ion etching of the 
AIGaAs. The charge concentration in the 
quantum wires could be increased by 
applying a positive bias to a back-side 
contact or by illumination. (There was no 
persistent photocurrent in these samples). 

Figure 8 shows the measured drain-source 
current versus backgate bias for the quantum 
wires (solid line) and for an equivalent 2-D 
device (dashed line). To explain the struc¬ 
ture, we have developed a semi-classical 
model to calculate the conductivity of these 
devices. As seen in the dotted line in figure 
8, our calculations match the measured data 
quite well. 

Based on our calculations, we conclude that 
the observed mobility modulation is associ¬ 
ated with populating the quasi-one-di- 


< 



Figure 8. Drain-source current (l D s) as a function of 
substrate bais (V S ub) for the multiple parallel quantum 
wires (solid line) and scaled down 2D channel (dashed 
line). The inset shows the theoretical density of states 
in a Q1 D wire at O K. The dotted line is an analytical 
calculation. 

mensional subbands. The negative 
differential transconductance regions arise 
when higher subbands begin to be popu¬ 
lated. Currently, we are working on fabri¬ 
cating new versions of these devices that 
allow for control of the confining potential, 
and hence the position of the subbands, by 
use of a top gate. It is hoped that these 
devices will verify the validity of our model. 

1.9 Arrays of 

Field-Effect-Induced Quantum 
Dots 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
U S. Air Force - Office of Scientific Research 
Grant AFOSR-88-0304 

Project Staff 

Professor Dimitri A. Antoniadis, Martin Burkhardt, 
William Chu, Professor Jesus A. del Alamo, Reza 
A. Ghanbari, Khalid Ismail, Professor Marc A. 
Kastner, C.T. Liu, Professor Terry P. Orlando, Pro¬ 
fessor Henry I. Smith, T.P. Smith, Professor Dan 
Tsui 

A metal grid on a modulation-doped 
AIGaAs/GaAs substrate (depicted in figure 


14 RLE Progress Report Number 132 








Chapter 1. Submicron Structures Technology and Research 


9a) produces a two-dimensional periodic 
potential modulation at the AIGaAs/GaAs 
interface via the Schottky effect. If a gate 
electrode is attached to the grid, the potential 
can be further modified with an external 
voltage source. By changing the gate 
voltage from positive to negative values, the 
potential seen by the electrons located at the 
AIGaAs/GaAs interface can be varied from 
uniform (in which case the electrons behave 
as a 2-D electron gas) to weakly coupled 
zero-D quantum wells (figure 9b) to isolated 
zero-D quantum dots (figure 9c). We have 
made such structures, with spatial periods of 
200 nm in both orthogonal directions using 
technology similar to that described in 
Section 1.7, but now the grid gate occupies 
an area of several square millimeters. The 
isolated quantum dots and the attendant 
zero-dimensional electronic subbands were 
examined in collaboration with D. Tsui at 
Princeton University using far-infrared (FIR) 
cyclotron resonance. Transitions between 
the discrete energy levels in the quantum 
dots were observed as a function of magnetic 
field. Results were in agreement with a the¬ 
oretical model. 



E 



Figure 9b. Depiction of potential seen by electrons at 
the AIGaAs/GaAs interface for weakly coupled 
quantum dots. 



Figure 9c. Potential for the case of isolated quantum 
dots. 


In collaboration with T.P. Smith and K. Ismail 
at IBM we applied a magnetic field perpen¬ 
dicular to the plane of the grid gate and 
measured magnetocapacitance as the gate 
voltage was swept, changing the electronic 
system from two-dimensional to coupled 
quantum wells, to insolated quantum dots. 
The normal oscillations due to Landau level 
filling are strongly suppressed when the 
Fermi level lines up with the minibandgaps. 
Thus, these measurements give important 
confirmation of the minband structure we 
observed earlier in transport. In the weakly- 
coupled-quantum-dot regime the magneto¬ 
capacitance oscillations show evidence of 
fractal behavior when the number of flux 
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quanta per unit cell of our grid is a rational 
number (i.e., 1/2, 1/3, 3/5, etc.). 

We are currently fabricating a new set of 
grid-gate MODFETS, using an improved fab¬ 
rication process and will study their trans¬ 
port, capacitance, and absorption properties 
as a function of magnetic field. 


1.10 

Planar-Resonant-Tunneling 
Field-Effect Transistors 
(PRESTFETs) 

Sponsor 

U.S. Air Force - Office of Scientific Research 
Grant AFOSR-85-0154 

Project Staff 

Professor Dimitri A. Antoniadis, William Chu, 
Khalid Ismail, and Professor Henry I. Smith 

Encouraged by the observation of negative 
differential resistance in the surface- 
superlattice MODFET (Section 1.7), we have 
proceeded to fabricate and test a double¬ 
barrier MODFET shown schematically in 
figure 10. Under the gate electrodes, the 
energy level diagram is as sketched in figure 
11. With proper biasing, resonant tunneling 
should be observable, hence the name of the 
device PRESTFET (Planar-Resonant Tun¬ 
neling Field-Effect Transistor). 

Three types of measurements were performed 
(at 4.2 K) on a PRESTFET device in which 
the gate electrodes were 60 nm long and 
separated by 60 nm: (1) gate electrodes are 
connected together, the drain-to-source 
voltage is held fixed ( < 5 mV), and the gate 
voltage is swept; (2) the electrodes are con¬ 
nected together and the intensity of an illu¬ 
minating LED is scanned; (3) one gate bias 
is fixed while the other is swept. 



Figure 10 Sketch of the layout of a 4-terminal, 
double-barrier, planar-resonant-tunneling field-effect 
transistor (PRESTFET). 



Figure 11. Energy level diagram of a symmetric, 
double-barrier PRESTFET. 


In all three cases, clear resonances were 
observed in the drain-to-source current. In 
figure 12, a plot of the results of the third 
experiment is shown. The results of these 
three measurements provide clear evidence 
that resonant tunneling occurs through the 
quantized well states. 

We are currently fabricating PRESTFET 
devices with even finer electrodes and devel¬ 
oping a high-throughput PRESTFET fabri¬ 
cation process using x-ray nanolithography. 
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Figure 12. Drain-source current as a function of bias 
on gate 2 for three different bias conditions, below 
threshold, on gate 1. 

1.11 Submicrometer-Period 
Gold Transmission Gratings for 
X-Ray Spectroscopy and 
Atom-Beam Interferometry 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
X-Opt., Incorporated 

Project Staff 

Dr. Mark L. Schattenburg. Professor Henry I. Smith 

Gold transmission gratings with periods of 
0.2 to 1.0 /xm, and thicknesses ranging from 
0.1 to 1 jim are fabricated using x-ray litho¬ 
graphy and electroplating. The x-ray masks 
are made with holographic lithography. 
Transmission gratings are either supported on 
thin (1 /xm) membranes or are made self- 
supporting by the addition of crossing struts. 
They are supplied to external laboratories for 
spectroscopy of the x-ray emission from a 
variety of sources. Over 15 laboratories 
around the world are using MIT-supplied 
gratings, and this project constitutes the sole 
source for these diffractors. Self-supporting 
transmission gratings in both gold and thin - 
film chromium are being provided to Pro¬ 


fessor David Pritchard's group for use in 
experiments to study the diffraction of 
neutral sodium atoms by gratings. The 
sodium atoms have a de Broglie wavelength 
of ~ 17 pm. Clear demonstrations of atomic 
diffraction have been made with these 
gratings. The gratings can also be used to 
divide and recombine an atomic beam 
coherently, and may provide the easiest route 
to the realization of an atom wave interfer¬ 
ometer. 

1.12 High-Dispersion, 
High-Efficiency Transmission 
Gratings for Astrophysical 
X-Ray Spectroscopy 

Sponsor 

National Aeronautics and Space Administration 
Contract NAS8-36748 

Project Staff 

Professor Claude R. Canizares, Dr. Mark L. 
Schattenburg, Professor Henry I. Smith 

Gold gratings with spatial periods of 0.1 to 
1.0 /xm make excellent dispersers for high 
resolution x-ray spectroscopy of 
astrophysical sources in the 100 eV to 10 
keV band. These gratings are planned for 
use in the Advanced X-Ray Astrophysics 
Facility (AXAF) that will be launched in the 
mid-1990s. In the region above 3 keV, the 
requirements of high dispersion and high 
efficiency dictate the use of the finest period 
gratings with aspect ratios approaching 10 to 
1. To achieve this, we first expose a grating 
pattern in 1.0 /xm-thick PMMA over a gold 
plating base using x-ray nanolithography. 
Gold is then electroplated into the spaces of 
the PMMA to a thickness up to 1 /xm. Flight 
prototype gratings have been fabricated and 
undergone space-worthiness tests. In the 
initial s f age of this program we used the 
carbon K x-ray (A = 4.5 nm) which required 
that the mask and substrate be in contact to 
avoid diffraction. This, in turn, caused dis¬ 
tortion of the grating. To avoid this problem 
we are developing a new technology of 
microgap x-ray nanolithography using the 
copper L x-ray (A = 1.33 nm) and Si 3 N« 
membrane masks. 
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Electronic Materials 


Academic and Research Staff 

Professor Carl V. Thompson, Dr. Paul Evans, Professor Clifton G. Fonstad, Jr., Dr. En Ma, Dr. 
John Melngailis, Professor Henry I. Smith 

Visiting Scientists and Research Affiliates 
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Undergraduate Student 
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Technical and Support Staff 
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2.1 Modeling of 
Microstructural Evolution in 
Thin Films 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation 
U.S. Air Force - Office of Scientific Research 
Contract AFOSR 85-0154 

Project Staff 

Professor Carl V. Thompson, Harold J. Frost, 
Jerrold A. Floro, Yachin Liu 

We are developing analytic models for 
normal, secondary and epitaxial grain growth 


in continuous thin films as well as particle 
coarsening in discontinuous films. The 
effects of surface or interface energy 
anistropy play especially important roles in 
these processes. We have developed com¬ 
puter models for film formation by crystal 
nucleation and growth to impingement under 
a variety of conditions. The topology and 
geometry of grain structures have been 
shown to strongly depend on the conditions 
of film formation. We have also developed 
computer simulations for two-dimensional 
normal grain growth and secondary grain 
growth. 

We are developing analytic models and com¬ 
puter simulations in parallel with exper¬ 
imental studies in model systems. Computer 


1 MIT Lincoln Laboratory. 

2 Pontificia Universidade Catolica, Department of Materials Science and Metallurgy, Rio de Janiero, Brazil. 

3 Harvard-MIT Health Sciences Program. 

* Thayer School of Engineering, Dartmouth College, Hanover, New Hampshire. 

* IBM, Thomas J. Watson Research Center. 
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simulations have allowed us to explain 
microstructrual features which are generally 
characteristic of thin films, including: (1) 
stagnation of normal grain growth, (2) 
development of lognormal distributions of 
grain sizes and (3) abnormal grain growth 
which leads texture development or epitaxy. 
Grain growth models also provide a theore¬ 
tical context for research on microstructure 
engineering using impurities, ion 
bombardment, substrate-surface-topography 
and precipitates. 

Publications 

Thompson, C.V. "Coarsening of Particles on 
a Planar Substrate: Interface Energy 
Anistropy and Application to Grain 
Growth in Thin Films." Acta Metal/. 
36:2929 (1988). 

Thompson, C.V. 'Observations of Grain 
Growth in Thin Films.” In Microstructural 
Science for Thin Film Metallizations in 
Electronic Applications, 115. Eds. J. 
Sanchez, D.A. Smith, and N. Delanerolle. 
Warrendale, Pennsylvania: the Minerals, 
Metals, and Materials Society, 1988. 

2.2 Reliability and 
Microstructures of 
Interconnects 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
Semiconductor Research Corporation 

Project Staff 

Jaeshin Cho, Harold Kahn, Hai Longworth, Pro¬ 
fessor Carl V. Thompson 

We are investigating the use of precipitates 
to produce metallic films and interconnects 
with engineered microstructures. We deposit 
initially layered alloy films (e.g., films depos¬ 
ited in a sandwich structure with a pure W 
layer between pure Al layers). When heated, 
precipitates form in the center of these films. 
These precipitates lead to a drag force which 
impedes grain boundary motion and 


supresses normal grain growth. Eventually, 
at elevated temperatures, a minor fraction of 
the grains begin to break free of the precipi¬ 
tates and grow abnormally. This can ulti¬ 
mately lead to very large grains and also 
changes in the crystallographic texture of 
films. With this technique, it is possible to 
control the grain sizes, the distribution of 
grain sizes, and the distribution of grain ori¬ 
entations over very broad ranges. We have 
now observed these effects in a number of 
aluminum alloy systems including Al-Cu-Cr, 
Al-Ag-Cr, Al-Mn-Cr, Al-Mn, and Al-W. 

We are also developing new techniques 
which allow statistical characterization of 
failure of contacts and interconnects for 
silicon-based integrated cricuit technology. 
We are using these techniques to correlate 
failure rates and mechanisms with micro¬ 
structures of interconnect lines and contact 
diffusion barriers. We are investigating tech¬ 
niques for controlling microstructures in 
order to improve contact ana interconnect 
reliability, especially under conditions which 
can lead to electromigration. 

We have recently shown that interconnect 
lines with bimodally distributed grain sizes 
have drastically reduced reliabilities. Also, 
for lines with monomodally distributed 
(uniform) grain sizes, increasing the grain 
size (relative to line width) results in an 
increase of both the median time to 
electromigration-induced failure and the 
lognormal standard deviation in the time to 
failure, see figures 1 -3. The net result, in 
large populations of lines with monomodal 
grain size distributions, is little or no change 
in the time to the first failure. We have 
explained these results in terms of a "failure 
unit model" in which grain boundaries are 
taken to be the individual units which are 
responsible for the reliability of a line. The 
successful application of this model indicates 
the importance of the properties of individual 
grain boundaries in controlling interconnect 
reliability. We are now investigating the sta¬ 
tistical reliability of lines with single, con¬ 
trolled grain boundaries in order to 
understand in detail the type of microstruc¬ 
tural features which limit interconnect reli¬ 
ability. 
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Figure 1. Distribution of failure times for 
2.2-/rm-wide, AI-2%Cu-0.3%Cr lines: (a) as-deposited 
(MTTF = 9.8 h, DTTF = 0.197), (b) after annealing for 
30 min. at 450°C to develop a bimodal grain size dis¬ 
tribution (MTTF = 1.08 h, DTTF = 0.223), (c) after 
annealing for 20 min. a: 550°C to obtain mono-modally 
distributed large grains (MTTF = 26.3 h, DTTF = 
0.465), and (d) after annealing for 20 mm at 550°C 
before patterning, to obtain bamboo microstructures 
(MTTF = 35.0 h, DTTF = 0.580). 



Figure 3. Median time to failure and deviation in the 
time to failure vs. linewidth/grain size ratio for 
AI-2%Cu-0.3%Cr lines with monomodally distributed 
grain sizes. 

Publication 

Cho, J., and C.V. Thompson. “The Grain Size 
Dependence of Electromigration Induced 
Failures in Narrow Interconnects.” Appl. 
Phys. Lett. 54:2577 (1989). 


2.3 Epitaxial Grain Growth 


AT »t I 

U.S. Aimorce - Office cl Scientific R« 
ContrMt AFOSR 85-1154 A 


Project StafV I 

Jerrold A. Floro^fcofes»r CarL 
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Thompson, Pro- 


Figure 2. Sketches of the microstructures of the lines 
described in Figure 1, based on transmission electron 
microscopy. 
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Graduate Students 
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4.1 Chemical Reaction 
Dynamics on Semiconductor 
Surfaces 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 

Project Staff 

Professor Sylvia T. Ceyer, Dr. Kenneth B. Laughlin, 
David J. Gladstone, Marianne McGonigal, 
Michelle Schulberg, Gerald Cain 

The etching of semiconductor materials in 
halocarbon plasma environments is a 
complex chemical process. The plasma is 
used to produce species that are highly reac¬ 
tive with the semiconductor surface. 
Because many different reactive species, 
including radicals, ions, and highly excited 
neutral molecules are produced in a plasma, 
the mechanism of the etching reaction is dif¬ 
ficult to probe. Our program is aimed at 
using molecular beam reactive scattering 
techniques to systematically examine the role 
of each of these species in the etching 
process. 

The initial goal of this project was to 
measure the barrier to F 2 dissociative chem¬ 
isorption on Si(100). Contrary to popular 
belief, we found essentially no barrier to the 
dissociation of F 2 on the unfluorinated 
surface. The dissociative chemisorption of F 2 
on Si(100) proceeds with unit probability for 
translational energies of the incident F 2 mole¬ 
cules as low as 0.5 kcal/mol. However, there 
is a substantial barrier to dissociation above 
fluorine coverages of 0.5 monolayer. Higher 


fluorine coverages sufficient to produce the 
etch product SiF 4 require a higher transla¬ 
tional energy of the F 2 molecule. We have 
observed efficient etching of Si(100) at 300 
K for an incident kinetic energy of F 2 of 15 
kcal/mol. The angular and translational 
energy distributions of the SiF 4 — as meas¬ 
ured in our newly constructed scattering 
apparatus consisting of two beams incident 
on a surface and a rotatable, triply 
differentially pumped quadrupole mass spec¬ 
trometer — indicate that the last step of the 
reaction is the thermal desorption of the 
product SiF 4 molecule. 

We have established that silicon can be 
etched without the use of molecules incident 
with energies of hundreds of electron volts, 
such as those present in plasma etching 
environments. Molecular beam techniques 
utilize molecules with energies only slightly 
higher than thermal (<1 eV), and therefore 
do not introduce radiation damage or defects 
into the Si lattice, which is a typical result of 
plasma etching. We are continuing our work 
on this system to understand the lack of F 2 
energy dependence on the formation of a 
second reaction product, SiF 2 . 

We have also recently observed a new kind 
of mechanism for dissociative chemisorption. 
In the limit of zero fluorine coverage, we 
observed that the Si surface strips one of the 
fluorine atoms from the incident F 2 molecule, 
leaving the other fluorine atom to scatter into 
the gas phase. Although this stripping 
reaction is analogous to the well-known 
stripping reactions in gas phase chemical 
reaction dynamics, this is the first observa¬ 
tion of its kind from a surface. 
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Chapter 5. Measurement of Electron-phonon 
Interactions Through Large-amplitude Phonon 
Excitation 


Academic and Research Staff 

Professor Keith A. Nelson 
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Dan E. leaird, 1 Dr. Andrew M. Weiner 1 

Graduate Student 

Gary P. Wiederrecht 

5.1 Introduction 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 

The main goal of this project is to measure 
the strengths of electron-phonon interactions 
in semiconductors to understand the effect of 
these interactions on electron mobility at dif¬ 
ferent temperatures. A secondary goal — 
extremely important for ultrafast electronics — 
is to determine the responses of electronic 
materials to input signals which are both fast 
(i.e., pulses of short duration) and high in 
repetition rate (i.e., many short pulses in 
rapid succession). 

5.2 High Repetition-rate 
Signals and Resonant 
Responses of Crystalline Solids 

We have successfully characterized the elec¬ 
tronic and vibrational responses of a wide 
variety of solids to single femtosecond 
pulses. With earlier JSEP support, we dem¬ 
onstrated that ultrashort optical pulses 
passing through nearly any material initiate 
coherent vibrational oscillations that influ¬ 
ence the material’s transient optical and elec¬ 
trical responses. This finding is important for 
fast electronics because it indicates the 


response of materials to fast electrical and 
optical signals. The material’s response to 
sufficiently short signal duration (i.e., optical 
or electrical pulse duration) can include tran¬ 
sient vibrational oscillations, negatively 
affecting its usefulness in possible applica¬ 
tions such as information processing. 

At the same time, the initiation of coherent 
lattice vibrational oscillations with ultrashort 
pulses provides us with a unique spectro¬ 
scopic tool. We can record femtosecond 
time-resolved “stop-action" spectroscopic 
observations of coherently vibrating crystal 
lattices, permitting characterization of well 
defined, nonequilibirium crystal structures at 
various stages of vibrational distortion. 
Characterization of the effects of lattice 
vibrations (phonons) on electonic energy 
levels is of particular importance in semicon¬ 
ductors since temperature-dependent carrier 
mobilities are influenced by electron-phonon 
interactions. 

Figure 1 shows the typical effects of a single 
ultrashort pulse on a crystalline solid. A fem¬ 
tosecond pulse exerts an "impulse" force on 
lattice vibrational modes, initiating vibrational 
oscillations. These oscillations influence the 
material's optical and electrical properties. 
The Terahertz-frequency oscillations in the 
data (shown in figure 1) correspond to indi¬ 
vidual lattice vibrational oscillation cycles, 
observed in real-time. 


1 Bell Communications Research (Bellcore). 
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Figure 1. Impulsive stimulated Raman scattering 
(ISRS) data in which optic phonon oscillations and 
decay in lithium niobate are monitored. In this 
noncentrosymmetric crystal, the optic phonon mode is 
both Raman and IR active, and is strongly coupled to 
the polariton mode associated with infrared light propa¬ 
gation. The phonon oscillations induced through ISRS 
are therefore accompanied by IR radiation at the same 
frequency. Femotosecond pulses have therefore been 
used to excite phonons and also to produce tunable 
Terahertz-frequency (far-IR) radiation. The phonon 
frequency and damping rate depend on the ISRS scat¬ 
tering angle. 

We are monitoring the material responses to 
high-density signals instead of a single signal 
for two main reasons. First, in any real- 
world high-speed signal-processing applica¬ 
tion, it is essential to process signals which 
are not only short in duration but also fre¬ 
quent, i.e., signals with high information 
density. It is important, therefore, to under¬ 
stand the differences between a material's 
response to one isolated pulse and its 
response to many pulses in rapid succession, 
the latter being far more relevant to applica¬ 
tions. 

Second, the data in figure 1 illustrate our 
ability to use short light pulse to control 
solid-state structure and behavior in a very 
precise way. We can use pulses to literally 
move ions or molecules in crystals toward or 
away from each other. In addition, we can 
monitor the motions and their effects in real 
time. Unfortunately, the extent of our optica! 
control is limited because, with a single 
pulse, we can only move material very small 
distances — less than 10~ 3 A. These small dis¬ 


tances have dramatic effects in signal pro¬ 
cessing, as figure 1 illustrates. However, if 
we could drive lattice vibrations much harder 

— for example, to amplitudes exceeding 0.1 A 

— we could induce permanent structural 
changes in solids, fabricate new materials, 
etc. We could also characterize electron- 
phonon interactions over a wide range of 
phonon displacements. 

To measure material responses to high 
repetition-rate sequences of femtosecond 
pulses, we first needed to generate the pulse 
sequences. We could have approached this 
problem by splitting a single femtosecond 
pulse many times with partial reflectors and 
aiming all the resulting beams into a sample 
with desired timing. This extremely cumber¬ 
some approach would have involved dozens 
of partial and high reflectors, lenses, and pre¬ 
cision delay lines. Instead, we chose to col¬ 
laborate with scientists at Bellcore who have 


Femtosecond Pulse Sequence 



Figure 2. Cross-correlation trace of a sequence of 
femtosecond pulses at Terahertz repetition rates pro¬ 
duced through pulse-shaping techniques. A single 
femtosecond pulse was passed through two gratings, 
two lenses, and a spatially varying mask to produce the 
output shown. The output can be used to exert repeti¬ 
tive 'impulsive' driving forces on selected vibrational 
modes to achieve larger amplitudes than can be driven 
with a single pulse. 
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Figure 3. ISRS data from the a-perylene organic molecular crystal. Several phonon modes are excited by the fem- 
tosecona excitation pulses. This gives rise to the characteristic "beating" pattern which contains sum and difference 
frequencies. The "spike" at t = 0 is due to the essentially instantaneous electronic response of the crystal to the 
excitation pulses. The inset shows that the vibrational signal intensity is only about 4% that of the spike. In some 
signal processing applications (especially those involving threshhold devices), the comparatively slow response due 
to vibrations could be disregarded. 

recently developed femtosecond optical 
"pulse-shaping” techniques through which a 
single femtosecond pulse can be converted 
into a sequence with the desired timing. In 
this technique, a single pulse, passed 
through a simple optical network consisting 
of only two lenses, two gratings, and a spa¬ 
tially varying mask, yields a pulse sequence 
similar to the one shown in figure 2. 

Sequences with different temporal profiles 
are selected by changing the mask. 

Figures 3 and 4, respectively, show crystal¬ 
line responses to single- and multiple-pulse 
excitation. For the initial demonstration, we 
used an organic molecular crystal since its 


vibrational and electronic properties were 
well understood. Figure 3 shows the 
response to a single femtosecond pulse. The 
“impulse" force exerted by the pulse drives 
several lattice vibrational modes whose 
simultaneous oscillations yield the "beating" 
pattern in the data. The "spike" in the data 
at t = 0 is caused by the instantanous elec¬ 
tronic response of the crystals to the 
excitation pulses. The signal due to lattice 
vibrations is only about 4 percent as strong 
as that due to the electronic response. In 
some signal processing applications, this 
"slow" signal could be disregarded, using 
instead the ultrafast electronic response as 
the basis for a rapid switch. 
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Figure 4. ISRS data from the perylene crystal driven by the pulse sequence shown in figure 2, which is tuned to 
drive the 80 cm" 1 lattice vibrational mode. At first the signal is dominated by the instantaneous electronic response 
of the sample to each of the pulses in the sequence. After about 10 pulses, signal due to repetitively driven lattice 
vibrations is apparent. Signal from the vibrational mode grows stronger with each successive pulse, eventually 
reaching intensity levels comparable to the strongest electronic response. Selective amplification of the 80 cm -1 
mode has been demonstrated The strong response due to the amplified vibrational mode could not be ignored in 
any signal processing application. 


Figure 4 shows the crystal's response to the 
sequence of pulses shown in figure 2. We 
have crafted this sequence so that the time 
between pulses (419 fs) is precisely equal to 
the oscillation period of one of the lattice 
vibrational modes. The impulse forces exerted 
by the pulses resonantly drive this mode in a 
manner analogous to repetitively pushing a 
child on a swing so that larger oscillation 
amplitudes are reached in time. The data 
show that for the first several pulses, the 
intense "spike" due to the electronic 
response dominates the signal. After about 
10 pulses, however, the amplified vibrational 
response is apparent. By the end of the pulse 
sequence (t > 3ps), the signal due to 


vibrations is as intense as the strongest elec¬ 
tronic "spike." This illustrates one important 
difference between the response of the mate¬ 
rial to high-density input signals versus the 
response to a single pulse. We could not 
ignore the long-lived vibrational response in 
any real-life signal-processing application. 

The results also illustrate the increased 
optical control over the sample afforded by 
judicious multiple-pulse excitation. In this 
case, through repetitive excitation, we 
increased the vibrational amplitude by almost 
a factor of 10. Amplification is also highly 
mode-selective, since the pulse sequence is 
timed to drive only one of the lattice 
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vibrations. This is clear from the single¬ 
frequency oscillations in the data, as com¬ 
pared to the beating pattern in figure 3 that 
results from excitation of several vibrational 
modes. 

At present, we are beginning experiments in 
which these techniques are applied to semi¬ 
conductors. We expect to drive large- 
amplitude vibrational responses and to 
observe their effects on electronic energy 
levels. This will permit characterization of 
the electron-phonon interactions. Another 
very interesting application will be to mul¬ 
tiple quantum well structures. Repetitive 
excitation should lead to large-amplitude 
responses of acoustic-like phonons which 
traverse a single well in the time between 
pulses of the sequence. 

It is worth noting that our primary outside 
collaborator in this project, Dr. Andrew M. 
Weiner, performed research for his Ph.D. 
under the supervision of Professor Erich P. 
Ippen. Under earlier JSEP support, we col¬ 
laborated with Dr. Weiner and Professor 
Ippen to conduct some of the first femto¬ 
second time-resolved spectroscopy exper¬ 
iments on molecular crystals. 

Our work in this area is widely recognized 
within the ultrafast optics and electronics 
community. Optics News described our 
research in a summary of 1989 highlights, 
and an article in Science is forthcoming. At 
present, we are preparing more detailed 
reports. 


Publications 

Nelson, K.A. "Impulsive Stimulated Raman 
Scattering with Single-pulse and 
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the 12th International Conference on 
Raman Spectroscopy. Forthcoming. 

Nelson K.A., and E.P. Ippen. "Femtosecond 
coherent spectroscopy." Adv. Chem. Phys. 
75: 1-35 (1989). 
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6.1 Chemical Beam Epitaxy 
Facility 

Sponsors 

3M Company Faculty Development Grant 
AT&T Research Foundation 
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National Science Foundation 

Grants ECS 88-46919 and ECS 89-05909 
Purdue University 
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U.S. Navy - Office of Naval Research 
Contract N00014-88-K-0564 

We made substantial progress in 1989 
;?ward establishing a chemical beam epitaxy 
(CBE) facility for the epitaxial growth of 
compound semiconductors. The major part 
of our effort was in the development and 
assembly of the CBE hardware including the 
growth, transfer, introduction, and bake 
chambers. By the end of 1989, the CBE 
system was approximately 85 percent com¬ 
plete. We expect delivery from the manufac¬ 
turer by the middle of February 1990. 

Figure 1 is a photograph of the CBE system 
showing its growth reactor module and some 
associated electronics. Within the black 
metal frame outlining the exhausted enclo¬ 


sure are the gas manifold (upper left side), 
ultrahigh vacuum growth chamber (right 
side), and the pumping subsystem (bottom 
left). Figure 2 shows the transfer, introduc¬ 
tion, and bake chamber modules. The black 
metal frame outlines a fully exhausted, fil¬ 
tered laminar flow enclosure. 

The CBE growth reactor is equipped with 
four individual gas injectors for metalorganics 
of Zn, Ga and In, Se, and for the hydride 
arsine. The temperature of each gas injector 
can be varied between 0 to 1250 degrees 
Celsius so that the molecules can be 
"cracked" or thermally decomposed prior to 
impingement on the heated single crystalline 
substrate. Three typical effusion cells, such 
as those used in molecular beam epitaxy, can 
be used with dopant or alloying elements. In 
situ analysis facilities include reflection high 
energy electron diffraction, a residual gas 
analyzer, and a quartz crystal flux monitor. 
Four symmetrically placed heated sapphire 
viewports can be used to direct a laser beam 
at incident angles to the growing surface 
front and to collect the reflected laser beam. 

The design of the original CBE system was 
expanded to include an additional ultrahigh 
vacuum analytical/metalization chamber, 
which will added in mid-1990. Attached to 
the transfer chamber, this new chamber will 


1 Professor, Department of Electronics, Kanazawa Institute of Technology, Kanazawa, Japan. 
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Figure 1. Growth reactor module and partially complete electronics rack of the chemical beam epitaxy system (as 
of November 1989). 


allow samples to remain in the UHV environ¬ 
ment before and after analysis or 
metalization. The additional chamber is sub¬ 
divided to include both in situ analysis and 
metalization capabilities. Initially, the 
chamber will contain Auger electron spec¬ 
troscopy (AES) and will be configured to 
have the associated ports, sample heater and 
manipulator, shutters, etc., necessary to add 
the electron beam evaporators to perform 
metalizations. Additional viewports will be 
provided to illuminate samples during Auger 
electron spectroscopic analysis and 
metalization with a coherent source of light. 
The source of photons will be provided by 
various lasers or lamps having a multitude of 
powers and wavelengths. 

We envision that the laser-aided investi¬ 
gations carried out in the analytical/- 
metalization chamber will complement 
laser-assisted experiments planned during the 


epitaxial growth of semiconductor layers. 
We anticipate that both types of experiments 
will provide clues to understanding the effect 
of photons on the physisorption, chemisorp¬ 
tion and incorporation processes which 
occur during the chemical beam epitaxial 
growth process. 

Funds for obtaining the analytical/- 
metalization chamber were provided by the 
National Science Foundation, Defense 
Advanced Research Projects Agency, and by 
a donation from Emcore Corporation. 

The space which will house the new CBE 
system was completely remodeled. Reno¬ 
vation included: (1) major demolition of 
walls to create one large (1200 square foot) 
lab from three smaller labs; (2) an increase 
in the services available such as power, 
water, nitrogen, etc.; and (3) installation of a 
very sophisticated exhaust system, a design 
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Figure 2. Front view of the bake chamber (left), introduction chamber (center), and transfer chamber (right). 


requirement for the custom CBE system. 
(Both the growth reactor module and the 
introduction chamber module are fully 
enclosed and exhausted.) The laboratory 
also includes a small entry room, typical for a 
clean room environment, and a substrate 
preparation area. The substrates can be 
chemically cleaned in dedicated solvent 
hoods followed by etching in dedicated ^wd 
hoods in adjacent laboratory space. The lab¬ 
oratory will also house a surface profilometer 
for studying film thickness and a Normarski 
interference microscope for studying surface 
morphology. 


6.2 Controlled Substitutional 
Doping of ZnSe Grown by 
Chemical Beam Epitaxy 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
U S. Navy - Office of Naval Research 
Contract N00014-88-K-0564 

The potential of the 11-VI compound semi¬ 
conductor family for optical device applica¬ 
tions has been relatively underutilized. There 
are substantial difficulties in preparing them 
by any growth technique to be 
amphoterically doped in a controlled manner. 
In the past, because bulk, liquid phase, or 
vapor phase growth techniques required very 
high temperatures for fabrication of ZnSe, 
control of the carrier type could not be 
achieved. The resultant ZnSe exhibited 
n-type transport properties or the as-grown 
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material was so heavily compensated that the 
material was highly resistive. Newer growth 
technologies such as molecular beam epitaxy 
(MBE) and chemical beam epitaxy (CBE) 
provide promise for circumventing the 
doping problems of the 11-Vis and ZnSe as 
one example, due to the non-equilibrium 
nature of the growth, and by their ability to 
affect the growth, with photon illumination. 
Once ZnSe is controllably doped as p-type 
and n-type material, many device applica¬ 
tions could take advantage of the wide, 
direct bandgap (2.67 eV at room temper¬ 
ature) which provides for the emission of 
blue light. Blue light emitting diodes, blue 
semiconductor injection lasers, and fiat panel 
electroluminescent display devices are repre¬ 
sentative of the many possible optical 
devices. The blue wavelength range will be 
useful in (1) underwater communication, (2) 
visible holography, (3) laser disc recording, 
where a shorter wavelength would allow for 
a ten-fold increase in the density of stored 
information, and (4) for short range commu¬ 
nications where new polymer fibers can be 
utilized. 

Studies of dopant incorporation reported in 
the literature for ZnSe were complicated by 
the presence of unintentional impurities 
found in both elemental and compound 
source material which is used in the MBE 
growth process. In most cases, undoped 
ZnSe grown by MBE, using commercially 
available source material, has been reported 
to be n-type with low resistivity (~ Ifi-cm). 
The low resistivity of the "undoped" ZnSe 
material implied that the ZnSe was of correct 
stoichiometry because a relatively small devi¬ 
ation from a Zn-to-Se flux ratio of 
one-to-one toward either Zn-rich or Se-rich 
conditions was found to result in high 
resistivity material. (The defects generated 
during growth under non-stoichiometric con¬ 
ditions appeared to compensate for the non- 
intentionaily incorporated impurities giving 
rise to the high resistivity.) 

Although major progress was recently made 
in the molecular beam epitaxy of ZnSe, con¬ 
trolled substitutional n- and p-type doping is 
still significantly limited. The technique of 
chemical beam epitaxy (CBE), however, 
combines the many advantages of molecular 
beam epitaxy and metalorganic chemical 
vapor deposition (MOCVD) growth tech¬ 


niques while eliminating their many disad¬ 
vantages. 


An ultrahigh vacuum environment, similar to 
that used in MBE, is utilized in CBE. 
Because in CBE the solid source material of 
MBE is replaced by gaseous source material, 
CBE is similar to MOCVD in that respect. 
The gas sources are controlled by precision 
mass flow controllers, permitting constant 
flux ratios at reasonable growth rates, uy 
employing metalorganic gaseous sources for 
the Zn and Se species, we will obtain a 
greater degree of control of the flux intensi¬ 
ties. The MBE technique of employing low 
temperature, high vapor pressure sources of 
Zn and Se cannot achieve this control. 
Because the elements originate from a gas 
bottle, long-term operation is possible since 
very little gaseous material is consumed. We 
will elimin?**^ many hardware problems 
encounte ud dnriny the growth of 11 -VI com- 
r>oj:ids by using CBE. We also anticipate 
that the use of the CBE technique will allow 
us to solve crystal growth problems. Scien- 
t ; sts and device engineers will be able to 
capitalize on the many attractive material 
properties of the ll-VI compound semicon¬ 
ductors. 


With research support from the Office of 
Naval Research, we will explore the CBE of 
ZnSe. We will determine its optimum growth 
parameters following microstructural, optical, 
and electrical evaluation of the epitaxial 
layers. After we achieve high quality, ultra- 
pure, undoped ZnSe, our objective under the 
program supported by the Joint Services 
Electronics Program will be to incorporate 
the impurity species in ZnSe. 
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7.1 Introduction 


Sponsors 
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Charles S. Draper Laboratory 
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The goal of this project is to design, fabri¬ 
cate, test, and model submicron 
InAIAs/InGaAs Heterostructure Field-Effect 
Transistors (HFETs) on InP. These devices 
are of great interest for applications in long- 
wavelength optical communications and 
ultra-high frequency microwave telecommu¬ 
nications. We expect that extraordinary 
device performance will result from the 
unique combination of: (1) the large con¬ 
duction bandgap discontinuity between 
InAIAs and InGaAs, (2) the large T — L sep¬ 


aration in InGaAs, and (3) the high velocity 
of electrons in InGaAs. 

We have recently pioneered HFETs in which 
the InGaAs channel is heavily doped but the 
InAIAs insulator is undoped (MIDFETs for 
Metal-lnsulator-Doped semiconductor FET). 1 
We have found that these devices display a 
performance comparable to InAIAs/InGaAs 
Modulation-Doped FETs (MODFETs) of 
similar gate length. In addition, these 
devices offer a number of unique benefits not 
found in other device structures: reduced 
transconductance 1 and f T collapse, 2 higher 
breakdown voltage, 1 and freedom to optimize 
gate insulator parameters. 3 A new 
modulation-doped channel MIDFET also dis¬ 
plays unprecedented microwave character¬ 
istics. 4 

During 1989, we obtained the first working 
InAIAs/InGaAs MIDFETs on InP fabricated 


1 J.A. del Alamo and T. Mizutani, “An lno 52 Alo* 8 As/n + lno 53 Gao« 7 As MISFET with a Heavily-Doped Channel," IEEE 
Electron Device Lett. EDL-8 (11): 534-536 (1987). 

2 J.A. del Alamo and T. Mizutani, 'Bias Dependence of It and f m »* in an lno 62 Alo 48 As/n + -lno. 63 Gao 47 As MISFET," 
IEEE Electron Device Lett. EDL-9 (12): 654-656 (1988). 

3 J.A. del Alamo and T. Mizutani, 'A Recessed-Gate lnos 2 Alo* 8 As/n + -lno, 53 Gao 47 As MIS-type FET," IEEE Trans. 
Electron Devices ED-36 (4): 646-650 (1989). 

« J.A. del Alamo and T. Mizutani, "An ln 0 s 2 Alo 48 As/n*-ln 0 s 3 Gao 47 As MISFET with a Modulation-Doped Channel," 
IEEE Electron Device Lett. EDl-10 (8):394-396 (1989). 
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at MIT. 6 The performance and properties of 
these MIDFETs are similar to those fabricated 
elsewhere. 3 More relevant to the problem, we 
fabricated, for the first time, strained - 
insulator InAIAs/InGaAs MIDFETs. This 
report gives a detailed description of our 
research on these devices. 

We have exploited the inherent flexibility of 
gate insulator design of a MIDFET by 
increasing the AlAs content in the InAIAs 
insulator from the lattice-matching composi¬ 
tion to InP (lno 52 Alo 48 As), straining it towards 
larger bandgaps. This rapidly effects an 
enhancement in the conduction band discon¬ 
tinuity (AE C ) between the InAIAs and the 
InGaAs channel (figure 1). We expect, in 
this manner, to obtain the following benefits: 

• reduced gate leakage 

• improved gate breakdown voltage 

• reduced real-space transfer of hot elec¬ 
trons from the channel to the gate. 

We have demonstrated all of these benefits 
experimentally while largely unaffecting the 
rest of the device parameters (threshold 
voltage, transconductance, etc.). 


7.2 Experiments 

We grew four wafers by molecular-beam 
epitaxy in MIT's Riber 2300 system. Fol¬ 
lowing standard procedures, we prepared the 
surface of the InP wafer. Figure 2 shows the 
cross section of the grown heterostructures. 
The layer structure is similar to those utilized 
previously, 1 except for: 

1. The channel design is more aggressive: 
doping is higher, 4x10 18 cnrr 3 , and 
thickness is lower, 100 A. This channel 
is expected to result in superior perform¬ 
ance. 

2. The Al composition in three of the wafers 
is increased over lattice matching condi¬ 
tions, i.e., wafers were grown in which x 
in the ln ( i_ x| AI,As layer is 0.48, 0.52, 0.60, 



Figure 1. Cross-sectional equilibrium band diagram of 
device structure indicating the increase in AE g with 
increase in Al composition. 


and 0.70 (all insulator layers are 300 A 
thick). This is the main purpose of the 
present experiment. 

We continued fabrication after growth with 
mesa isolation utilizing a selective chemical 
etchant, AuGeNi ohmic contact formation, 
TiAu interconnect, and TiAu gate layer for¬ 
mation. Figure 2 shows the resulting device 
cross section. 

We have integrated a number of devices of 
varying dimensions and a large set of test 
structures into a 2.9 x 2.7 mm 2 chip. We 
present here DC l-V measurements on 
MIDFETs with a nominal gate length of 1.5 
/im and a gate width of 30 ^m. 


7.3 Device Results 

Figure 3 shows a plot of the transconduc¬ 
tance versus gate to source voltage for four 
representative devices with varying 
ln (1 _ x| Al„As insulator compositions. As 
expected from simple theory, the threshold 
voltage is, to the first order, unaffected by 
the introduction of strain in the insulator. 
The lattice-matched device has a maximum 
transconductance of 220 mS/mm. This value 
is higher than for any of the InAIAs/InGaAs 
MIDFETs previously reported in the literature 
with a similar gate length and overall 
epitaxial structure (the best was 160 


8 S. Bahl, W.J. Azzam, and J.A. del Alamo, "Fabrication of an lnAIAs/n + -lnGaAs MIDFET," unpublished paper pre¬ 
sented at the Fourth New England MBE Workshop, Cambridge, Massachusetts, May 4,1989. 


54 


RLE Progress Report Number 132 











Chapter 7. High Frequency Semiconductor Field-Effect Transistors 





50 A 

InGaAs 

300 & 

In ( i. x) Al x As *=o.a, o.m, o.6, o.r 

ioo A 

nMnGoAs (4xio“cm 3 ) 

ioo A 

InGaAs 

iooo A 

InAIAs 


SI-InP 


Figure 2. Cross section of device structure. 

mS/mm). 3 We believe that is due to the 
enhanced channel doping in this new device 
structure. 

An increase of AlAs composition in the 
ln ( i_ x( Al x As gate insulator from the lattice- 
matched condition x = 0.48 to x = 0.6 
brings about a slight reduction of peak g m 
due to an enhanced contact resistance. First 
order theory predicts this result. A further 
increase of x to 0.7 results in a catastrophic 
reduction in g m that cannot be explained by a 
larger source resistance. It appears that, from 
an electrical point of view, the critical layer 
thickness for the 300 A ln (1 _ x ,Al x As insulator 
has been exceeded in going from x = 0.6 to 
x = 0.7. 

While there is a small reduction of g m with 
AlAs composition, the g m curve broadens in 
V g „ as figuie 3 indicates, because gate 
current is being suppressed due to the 
enhanced conduction band discontinuity 
(discussed below). Thus, the maximum 
current driving capability of the devices is 
not degraded. In fact, this capability 
improves with AlAs composition up to 0.6. 
The maximum drain current obtained is 300 
mA/mm for the x = 0.6 device. This is better 
than that which had been reported previously 
in similar devices and is a very competitive 
value for any 1.5 nm gate length HFET. 



Figure 3. Transconductance versus gate to source 
voltage for Vd S = 4 V. 

Increasing the AlAs composition results in a 
drastic reduction of forward gate current. 
Figure 4 illustrates gate current as a function 
of drain-source voltage for V 8t = 1.2 V. As 
the AlAs composition increases, gate current 
at all values of V ds is drastically suppressed. 
This results in broader g m versus V g5 charac¬ 
teristics (shown in figure 3) and higher l Dmax , 
as mentioned above. 

A more significant result, shown in figure 5, 
is the impact of AlAs composition on the 
real-space transfer of hot electrons from the 
channel to the gate. Such a phenomenon is 
very clearly observed in the devices with 
AlAs composition of 0.48 and 0.52. Its 
unique signature is an increase in the gate 
current as V d5 increases for a given V a ,. This 
increase is due to electrons that are heated 
by the high electric field prevalent inside the 
channel and then acquire sufficient energy to 
overcome the conduction band discontinuity 
between the channel and the gate insulator. 
The drain current saturates when the device 
itself enters saturation. An AlAs composition 
of 0.6 results in a sufficient conduction band 
discontinuity to completely suppress this 
phenomenon. This may be a major advan¬ 
tage of submicron strained-insulator 
MIDFETs in which the electric field in the 
channel can reach a very high value. 

Figure 5 shows that real-space transfer of hot 
electrons to the gate causes a small negative 
differential output resistance. Because elec- 
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Figure 4. Gate current versus V<j$ at V gs = 1 2 V, as a 
function of AlAs composition. 


trons are lost to the gate, this phenomenon 
results in drain current reduction. In addi¬ 
tion, it may produce device instabilities and 
oscillations. For the x = 0.6 device, real- 
space transfer is completely suppressed and 
the drain current reaches higher values, as 
figure 5 indicates. Overall, the x = 0.6 
device has ten times less gate current than 
the x = 0.48 device at practical operating 
points. 

A final advantage of the strained-insulator 
MIDFETs is the enhanced reverse-bias gate 
breakdown voltage. Figure 6 shows the gate 
l-V characteristics measured with drain and 
source shorted. We obtain a rather hard 
breakdown for all devices. If we define 
breakdown at a reverse current of 1 mA, the 
breakdown voltage for the lattice matched 
device is 17 V. As the AlAs composition is 
increased in the insulator, the breakdown 
voltage increases very quickly, up to 28 V for 
x = 0.7. 

The combination of the unique features 
investigated above makes strained-insulator 
InAIAs/n^-InGaAs MIDFETs very suitable 


ID 


(mA) 


*> 



Figure 5. Drain current versus V ds at V 9S = 1.2 V, as a 
function of AlAs composition. 

candidates for high-power microwave appli¬ 
cations. Figure 7, for example, shows that 
ln 0 4oAlo6oAs/n + -lnos3Ga 0 47As MIDFETs dis¬ 
played excellent l-V characteristics up to 
V dt = 10 V. 


16 


(mA) 



Figure 6. Reverse gate current characteristics of FETs 
as a function of AlAs composition V dt = 0 V. 
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Figure 7. I-V characteristics of fabricated 
• n 0 4oAlo6oAs/n*-lno53Ga 0 47As MIDFET. Scales: Id: 
1mA/div,V ds :1 V-div,Vg,:0.2V/step, V g ,( max) = 0 V. 
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8.1 Project Description 
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In this program, we examine the behavior of 
the superconducting Josephson channel in 
electronically gated resonant tunneling struc¬ 
tures. These hybrid superconduc¬ 
ting/semiconducting structures represent the 
first attempt at Josephson coupling via reso¬ 
nant tunneling. A significant consequence of 
this approach is that the quantum confine¬ 
ment levels, not the semiconducting gap 
voltage, set the energy scale for modulating 
the Josephson coupling. Having made sig¬ 
nificant progress in the fabrication of such a 
device, we expect to have test devices avail¬ 
able for initial measurements in coming 
months. 

The envisioned device will use a thin silicon 
layer (<1000 A) as a superconducting weak 
link between the superconducting source and 
drain electrodes. The superconducting coun¬ 
terelectrodes are separated from the silicon 
layer by a thin oxide layer. The weak link 
also defines a quantum well in which the 
superconducting Cooper pairs can interfere 
with each other (either constructively or 
destructively). This creates a resonant tun¬ 
neling channel for Cooper pairs (i.e., the 
Josephson channel) between the source and 
drain in parallel with the normal quasiparticle 
(or single electron) channel. The applied 
gate voltage biases the levels within the 


semiconducting well, modulating the 
Josephson current by varying the wavelength 
of the tunneling electrodes. In many ways, 
this structure is analogous to a Fabry-Perot 
interferometer. Previous efforts on similar 
hybrid Josephson devices have not utilized 
either resonant tunneling or quantum con¬ 
finement. Thus, there is a mismatch between 
the semiconducting and superconducting 
gap energy scales. 1 The key motivation for 
our devices is that the energy scale is set by 
quantum confinement and can be tuned sep¬ 
arately from either gap energy. For silicon 
well widths in the range of 500-1000 A, the 
energy levels separation in the well can be 
set more than ten times smaller than the 
semiconducting gap. 

Current progress focuses on two areas. The 
first is the development of ultrathin uniform 
oxide tunneling barriers for the resonant 
quantum well. This tunnel barrier must be 
thin enough to enable a substantial number 
of Cooper pairs to tunnel through it, yet be 
thick enough to provide adequate isolation 
between the gate and the output. We ha»u 
successfully fabricated ultrathin oxide layers 
of very high quality, grown at 800°C in a 
dilute (0.8 percent) oxygen atmosphere. 
Oxide layers have been fabricated with thick¬ 
nesses of 15, 20, 25, 30 and 35 A, and all 
appear to have healthy barrier heights with 
no observable trap states (from C-V anal¬ 
ysis). 

The second area of progress is the examina¬ 
tion of two different device structures. The 
first device structure has a vertical geometry 


* A W. Kleinsasser, T.N. Jackson. G.D. Pettit, H. Schmid, J.M. Woodall, and D P. Kern, ‘Prospects for Proximity 
Affect Superconducting FETs,' IEEE Trans. Mag. MAG-25 (2): 1274-1277 (1989) and references within. 
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and uses a silicon membrane as the weak 
link, a geometry that has been fabricated 
before albeit for much thicker silicon mem¬ 
brane thicknesses. The second structure is a 
novel new structure, done in collaboration 
with Professor Henry I. Smith. This structure 
has a lateral geometry with a silicon wall as 
the weak link. Both are shown in figure 1. 
The operational physics of the two structures 
is identical. 



Figure 1. (a) Vertical device geometry, and (b) lateral 
device geometry Note: drawings are not to scale. 


In the vertical structure, the membrane is 
formed by an isotropic etch in ethylenedia- 
mine-pyrocatechol. This etch stops on a 
heavily boron-doped 1.75 nm thick layer 
deposited at the surface of the wafer. After 
this etch, the membrane is thinned to the 
desired thickness. Two disadvantages of this 
structure are its extreme fragility and its 
process limitation, being only heavily boron- 
doped. The advantage of the structure is 
that its fabrication requires only one level of 
conventional VLSI optical lithography. 
However, we have found the fabrication of 
membranes much thinner than 1 ^m to be 
exceedingly difficult, primarily due to 
fabrication-induced mechanical stresses. We 
have therefore stopped work on this partic¬ 
ular structure. 

The second device structure has a lateral 
topology, as shown in figure 1(b). The 


starting material for this device is (100) Si. 
The wall is formed by a strongly anisotropic 
potassium-hydroxide wet etch. As this weak 
link is defined by a mask layer, no etch stop 
is necessary. Therefore, this fabrication route 
places no constraints on the dopant type or 
level in the wafer, which is very desirable. 
Also, the mechanical integrity of this device 
is substantially more robust. In addition, the 
lateral structure allows it to be more easily 
integrated into circuits. 

The primary difficulty in fabricating this 
lateral device is the production of the mask 
and the alignment of the mask pattern to the 
crystalline axes of the underlying wafer. We 
believe that both of these problems are sur¬ 
mountable, and our work is progressing on 
both fronts. Indeed, although we need to 
make some refinements, we have successfully 
produced thin silicon walls on test chips 
(figure 2). After finalizing the fabrication 
route, we will then begin careful measure¬ 
ment and analysis of the device behavior. 



Figure 2. Nanow (500A) wall structure (on right) 
etched into Si<110> wafer. Note scale at bottom of 
figure. 
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9.1 Introduction 

The broad objective of our research effort is 
to develop lll-V quantum heterostructures for 
high performance electronic, optoelectronic, 
and photonic devices for high speed optical 
communications and signal processing. To 
this end, we are developing (1) new, higher 
performance materials systems including 
InP-based InGaAIAs heterostructures and 
<111 > oriented strained layer superlattices; 

(2) a new family of quantum-well-base, 
tunnel-barrier n-n-n transistors and near- 
and far-infrared optoelectrica! devices; and 

(3) new damage-free in situ processing tech¬ 
niques for fabricating advanced quantum 
structures and embedded heterostructures. 

The following sections describe our progress 
during the past year in the research areas 
listed above. Our group works closely with 
Professors Hermann Haus, Erich Ippen and 
James Fujimoto to develop the optical device 
application, characterization and modeling 
aspects of this program and with Professor 


Sylvia Ceyer to develop new in situ pro¬ 
cessing techniques. 

9.2 Computer Controlled 
Growth of Lattice-Matched 
InGaAIAs Heterostructures on 
InP 

Sponsors 

Charles S. Draper Laboratory 
Contract DL-H-315251 
Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation Grant 
Grant EET 87-03404 

Project Staff 

James Vlcek, Professor Clifton G. Fonstad 

The InGaAIAs material system, lattice- 
matched to InP substrates, includes alloy 
compositions with bandgaps that match the 
low-loss windows of silica-based optical 


1 Thomson CSF. 

2 Northeastern University, Boston. 
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fibers at 1.3 and 1.55 microns. Since the 
constituent "elements" of the pseudo-ternary 
alloy InGaAIAs (InAs, GaAs, and AlAs) each 
have the same anion, this compound is better 
suited to traditional molecular beam epitaxial 
(MBE) growth techniques than the more 
commonly used InGaAsP alloy system. 

In electronic and optical semiconductor 
devices, the need for both graded- 
composition and hyper-abrupt metallurgical 
junctions frequently arises, but it is difficult 
to build these structures with solid-source 
Molecular Beam Epitaxy. In the case of 
graded junctions, the fluxes of the constit¬ 
uent elements must be precisely controlled 
during the growth process to maintain lattice 
matching of the epitaxial layer to the 
substrate, while varying the ratios of two or 
more beam fluxes to each other. This control 
is accomplished by ramping the setpoint 
temperature of the effusion cells for the con¬ 
stituent elements of the growing crystal. 
Achieving hyper-abrupt junctions, on the 
other hand, is complicated by use of 
mechanical shutters to block the beam fluxes 
during the growth process. Opening and 
closing of these shutters changes the thermal 
characteristics of the effusion cell, leading to 
transients in the emitted flux which must be 
removed to achieve truly abrupt metallurgical 
junctions. 

Complications encountered in the growth of 
both graded layers and hyper-abrupt junc¬ 
tions can be addressed by a common sol¬ 
ution: construction of a model for the 
effusion cell which takes into account not 
only the nonlinear relationship between flux 
and temperature, but also the temporal 
response of the cell to setpoint changes. This 
model is then used to construct the proper 
setpoint temperature profiles to achieve 
either a time-varying flux profile for a graded 
layer, or a flat flux profile to remove a shutter 
transient. These temperature profiles are exe¬ 
cuted by a digital computer which communi¬ 
cates with the feedback-loop (or PID, for 
Proportional Integral Derivative controllers 
which stabilize the effusion cells). 

With this scheme, we have grown composi- 
tionally graded layers with maximum devi¬ 
ation from the desired flux profile maintained 
below one percent, reducing flux transients 
from values in the range of 10 to 25 percent 


to less than 1.5 percent. Current research on 
this topic includes extending models of the 
effusion cell and eliminating overlapping 
shutter transients similar to those that would 
occur in superlattice structures (multiple, 
closely-spaced abrupt heterojunctions). 

9.3 InGaAIAs Multiple 
Quantum Well Heterostructures 
for Guided Wave Optics 

Sponsors 

National Science Foundation Grant 
Grant EET 87-03404 

Project Staff 

James Vlcek, Professor Clifton G. Fonstad, Jr., in 
collaboration with Mary R Phillips and Professor 
Hermann A. Haus 

Through the use of quantum confinement at 
InGaAs/lnAIAs heterojunctions, emission and 
absorption edges of epitaxial grown struc¬ 
tures can be continuously varied over the 
wavelength range of 1.3 to 1.6 microns, 
making them compatible with silica fiber 
communications. In addition to use in semi¬ 
conductor lasers for optical emission and 
photodetectors for signal detection, these 
heterostructures can be employed for all- 
optical devices such as waveguides and 
modulators. The dynamics of optical nonlin¬ 
earities in the transparent region of these 
structures will be an important factor in the 
speed of refractive, all-optical switches. 

To measure subpicosecond optical 
nonlinearites just below the superlattice 
bandgap, we have fabricated InGaAs/lnAIAs 
multiple quantum well heterostructures. In a 
pump-probe experiment, we collaborated 
with Mary R. Phillips and Professor Hermann 
A. Haus of RLE to measure the refractive and 
absorptive nonlinearities of these samples. 

A typical heterostructure consisted of 72 
periods of an InGaAs/lnAIAs superlattice, 
each with periods of 70 A. This structure was 
grown using the computer-automated MBE 
system described in Section 9.2. Room tem¬ 
perature infrared absorption measurements of 
our best sample showed clear, strong 
excitonic absorption lines for both n = 1 and 
n = 2 transitions. Both the n = 1 light and 
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in the literature. Now, we are in the process 
of optimizing our layer structure, growth 
conditions, and device fabrication technology 
which should greatly enhance device per¬ 
formance. 

In addition to multiple quantum well laser 
diodes, we are investigating the use of 
graded-index optical confinement structures 
for (1) an optimal gain profile with minimal 
angular dispersion of the emitted beam and 
(2) structures utilizing controlled amounts of 
lattice mismatch strain in the active layer to 
further enhance the quantum size effects of 
the confined carriers. To achieve our long 
range goal of producing ultra-narrow 
linewidth devices, we are extensively ana¬ 
lyzing various theoretical techniques for 
reducing laser linewidth. Our ultimate goal is 
to produce room temperature laser diodes 
with emission linewidths below 100 Lrtz. 


9.5 InGaAIAs St rained-Layer 
Heterostructures on 111 GaAs 
and InP for Optical for 
Modulator Applications 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation Grant 
Grant EET 87-03404 

Project Staff 

Rhrtard Singer, Professor Clifton G. Fonstad, Jr., 
in collaboration with Stuart D. Brorson, Daryl 
Smith, 3 B. Laurich, 3 Christian Mailhoit, 4 Bruce D. 
McCombe, 5 and Bernie A. Weinstein 6 

When epitaxial layers with a slight lattice 
mismatch are deposited upon (111)-oriented 
GaAs or InP substrates, piezoelectric fields 
are generated. Whether a field points 
<111>Aor <111 > B is a function of whether 
a given epilayer is under tensile or 
compressive strain. Our group was the first 


to produce these structures and to demon¬ 
strate experimentally the existence of 
piezoelectric fields. We are now applying 
this effect to optoelectronic devices. 

If a quantum well structure is grown through 
appropriate choice of ternary composition 
and layering, the strain-induced field will 
red-shift the excitonic resonance of the well. 
This is the familiar Quantum Confined Stark 
Effect (QCSE). In the present structures, the 
internal strain-induced field provides a bias 
on the order of 105 V/cm. An applied 
voltage will either enhance or counteract the 
strain-induced field around a bias point in a 
highly sloped, nearly linear region of the field 
verses index and/or absorption characteristics 
of these structures. Thus, the external bias is 
expected to dramatically affect the transmit¬ 
tance or reflectance of the multiple quantum 
well structure at a given wavelength. Prac¬ 
tical, low voltage modulation of optical 
signals should be possible. 

Presently, we are seeking to grow, fabricate, 
and characterize window and waveguide 
p-i-n modulators in both the GaAs and InP 
material systems. Material of good surface 
morphology has been successfully grown by 
molecular beam epitaxy on (111) GaAs, and 
samples of (111 )-oriented strained layer het¬ 
erostructures have been distributed to a 
number of researchers interested in studying 
them. Recently, reflection modulators con¬ 
sisting of a mirror region of quarter wave 
plates and an active region of strained mul¬ 
tiple quantum wells have been grown and 
tested with no applied bias. Work is cur¬ 
rently underway to process these devices so 
that we can study their response to applied 
voltages. 

MBE growth on InP is more difficult because 
of the conflicting requirements of a low 
substrate temperature and high As overpres¬ 
sure to stabilize the volatile (111) surface, on 
the one hand, and the opposite growth con¬ 
ditions needed to enhance surface mobility 
of Al on the other. To overcome this diffi- 


3 Los Alamos National Laboratory. 

4 Lawrence Livermore National Laboratory. 
6 State University of New York at Buffalo. 
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culty, we will attempt a variation of migration 
enhanced epitaxy in which we will keep both 
the substrate temperature and As overpres¬ 
sure low. 

9.6 Molecular Beam Epitaxial 
Growth on (hi 1) Vicinal 
Surfaces 

Sponsor 
MIT Funds 

Project Staff 

Dr. Elias D. Towe, Professor Clifton G. Fonstad, 
Jr., in collaboration with H.Q. Le, 6 and J.V. 
Hryniewicz 6 

Traditionally, lll-V compound semicon¬ 
ductors have been synthesized on (100) 
substrate surfaces. Because it is desirable to 
fabricate devices on chemically etch- 
structured (100) surfaces, an understanding 
of molecular beam epitaxial growth on high- 
index planes vicinal to the (100) surface is 
important. Understanding the properties of 
material grown on such surfaces is also nec¬ 
essary. The particular planes of interest are 
those denoted by (hi 1). 

We have continued our work on the growth 
of epi-layers on (hi 1) vicinal surfaces of 
GaAs. In particular, the growth conditions 
and characteristics on the (111) B, (211) B 
and (511) B surfaces for acceptable quality 
layers are fairly well understood. GaAs and 
(Al, Ga)As layers with luminescence supe¬ 
rior to that of layers on the (100) surface 
have been grown on these surfaces. 

Our interest in the epi-lavers on the (hi 1) 
surfaces is due to the potential for improved- 
performance and for their use in novel optical 
devices. We are currently studying some 
laser devices with strained-layer active 
regions on (lll)-GaAs substrates. The 
devices have been fabricated, and we will 
test them soon. Some waveguide structures 
on (211) substrates have also been fabri¬ 
cated so that an undergraduate thesis student 


can evaluate them with a tunable Ti:sapphire 
laser. 

We expect that the threshold current densi¬ 
ties of the lasers on the (lll)-GaAs 
substrates will be lower than those on the 
conventional (100) substrates. This expecta¬ 
tion is based on the fact that the relative 
photoluminescence emission of the layers on 
these substrates is higher than that on the 
(100) substrates. 

9.7 Pseudomorphic 
AlAs/lnAs/InGaAs Resonant 
Tunneling Structures 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 

Project Staff 

Thomas P.E. Broekaert, Professor Clifton G. 
Fonstad, Jr. 

A fundamental objective in the design of res¬ 
onant tunneling diodes (RTDs) for micro- 
wave oscillators and digital switching is 
achievement of the highest possible peak 
current density at the lowest possible peak 
voltage, while maintaining a high peak-to- 
valley current ratio and wide valley region. 
The InGaAIAs material system, lattice 
matched to InP, has the design flexibility to 
optimize the RTD for microwave performance 
and switching applications in terms of peak 
current density and peak voltage of the RTD. 

Thus far, by using very thin, three monolayer 
AlAs barriers and InGaAs wells, within 
bounding regions of InGaAs lattice-matched 
to InP, we have achieved room temperature 
peak current densities in excess of 1 
MA/cm 2 , the highest reported value for RTDs 
to date, at peak voltages as low as 0.55 V. 
The peak-to-valley current ratio of these 
devices is approximately 2 to 4. 

For the diode structure, a turn-on threshold 
voltage is apparent in the l-V characteristics. 
This threshold characteristic can be removed 


* MIT Lincoln Laboratory. 
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by widening {he well or by including a low 
bandgap material, i.e., InAs, in the well. 
Diode structures of this type have been fabri¬ 
cated and display linear initial l-V character¬ 
istics with no turn-on threshold and a low 
peak voltage. 

The combination of thin AlAs barriers and 
InAs containing wells should allow us to 
further optimize the RTDs. 

9.8 Three-Terminal Quantum 
Well Base/Tunnel Barrier 
Devices 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 

Project Staff 

Thomas P.E. Broexaert, Lung-Han Peng, Professor 
Clifton G. Fonstad, Jr. 

We have gained our objective of making 
direct electrical contact to a populated con¬ 
duction band quantum well in a resonant 
tunneling configuration. Recently, in a sig¬ 
nificant breakthrough, we succeeded in 
etching away the layers above the InAs 
quantum well layer in an AlAs/lnAs/InGaAs 
resonant tunneling diode and making elec¬ 
trical contact directly to this 2 to 5 nm thick 
layer. This major advance in this technique 
opens the door to development of a new 
class of ultra-high performance electronic, 
optoelectronic, and photonic devices. We 
are currently beginning an aggressive 
program of research utilizing this advance to 
investigate quantum-well base, tunnel-barrier 
n-n-n transistors; tunnel injection, quantum- 
well-level laser diodes for the 2 to 5 micron 
region; sub-picosecond quantum-well-level, 
tunnel-barrier infrared photodiodes; and 
ultrafast quantum-well-level, tunnel-barrier 
optoelectronic modulators and non-linear 
optical components. 

The basis for all these devices is the 
AlAs/lnAs/InGaAs double barrier resonant 
tunneling structure discussed in the previous 
section. There are at least three unique fea¬ 
tures of this structure which make it ideally 
suited for these structures. First, the AlAs 


tunnel barriers are unusually high relative to 
the InGaAs injectors and the InAs well, 
which yields multiple confined well states 
and sharp quantum structure at room temper¬ 
ature. Second, the InAs well lies below the 
band edges of the InGaAs injectors so that 
when it is suitably wide and doped, the first 
well level will be populated. This provides 
all-important lateral conduction in the plane 
of the well and the necessary ground-state 
population for detector and optical 
modulator applications. Third, the very wide 
range of materials in these structures has 
allowed us to develop a selective etch which, 
in turn, makes it possible to make electrical 
contact to the very thin quantum well layer. 
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1.1 The Nonlinear Waveguide 
Interferometer 

Sponsors 

Charles S. Draper Laboratory 
Contract DL-H-404179 
Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Sciences Foundation 
Grant EET 87-00474 

Project Staff 

Professor James G. Fujimoto, Professor Hermann 
A. Haus, Professor Erich P. Ippen, Michael J. 
LaGasse, John D. Moores, Keren Bergman 


Section 1.2 describes a single-arm interfer¬ 
ometer invented in our laboratory that has 
been used for low power all-optical 
switching. We are currently investigating a 
switch very similar to the single arm interfer¬ 
ometer, which is based on a fiber (wave¬ 
guide) ring reflector, similar to, yet different 
from, the ring reflector demonstrated by the 
group at British Telecom. 2 In our research, 
we have shown that the fiber ring reflector is 
an excellent candidate for soliton switching 
and for producing squeezed light. 3 Conse¬ 
quently, since any experiments limited by 
quantum noise (i.e., detector shot noise) 
must take into account the quantum proper¬ 
ties of light, we have become interested in 
developing a self-consistent quantum theory 


1 NEC, Japan. 

2 N.J. Doran, K.J. Blow, and D. Wood, Proc. SPIE 836:238-243 (1987). 

3 J.D. Moores, K. Bergman, H.A. Haus, and E.P. Ippen, “Optical Switching Using Fiber Ring Reflectors,' submitted 
for publication; H.A. Haus, K. Watanabe, and Y. Yamamoto, 'Quantum-Nondemolition Measurement of Optical 
Solitons," J Opt Soc. Am. B 6:1138-1148 (1989); Y, Lai and H.A. Haus, 'Characteristic Functions and Quantum 
Measurements of Optical Observables," J. Quant. Opt., forthcoming; M. Shirasaki and H.A. Haus, "Squeezing of 
Pulses in a Nonlinear Interferometer,' J. Opt. Soc. Am. B 7(1): 30-34 (1990). 
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for optical solitons. 4 The promise of a simple 
squeezing mechanism using optical pulses in 
fibers has also kindled our interest in devel¬ 
oping a fiber gyro that would operate at a 
noise level below the standard quantum limit 
shot noise. 

The fiber ring is one realization of a fiber ring 
gyro utilizing the Sagnac effect. Unfortu¬ 
nately, the ring cannot be used directly as a 
gyro with squeezed noise, because the noise 
that is in phase with the Sagnac signal is not 
reduced below the vacuum level. However, 
our theory predicts that a cascade of two 
rings, one as a squeezer and one as a gyro, 
can be made to perform at a level below the 
standard quantum limit. Ongoing exper¬ 
iments should determine whether squeezing 
is achieved in the ring. The squeezing 
process is broadband so that pulses can be 
used in the experiments. Nonlinearity is 
enhanced because the duty cycle enhances 
peak power. Through our experiments, we 
need to identify the several classical sources 
of noise that could prevent squeezing. 

1.2 Picosecond Optical Signal 
Sampling 

Sponsor 

National Science Foundation 
Grant EET 87-00474 


Project Staff 

Professor James G. Fujimoto, Professor Hermann 
A. Haus, Wei-Ping Huang, Michael J. LaGasse 

The goal of our work is to develop optical 
waveguide devices for high rate all-optical 
switching. This objective involves develop¬ 
ment of (1) analytical and computer tools for 
the characterization of optical waveguide 
properties and (2) computer programs for 
describing optical pulse interactions and for 
experimentally verifying their predicted per¬ 
formance. 

Through our pioneering work, we have 
achieved a better understanding of coupling 
modes theory and prediction of crosstalk in 
waveguide couplers. 6 Our paper on this 
subject was selected for publication in an 
IEEE volume of reprints on optical couplers. 

Along with our work on coupled mode 
theory, we have clarified the theory of 
transfer of power in grating assisted 
couplers 6 and have developed a self- 
consistent theory for tapered couplers. 7 In 
addition, we have developed simple com¬ 
puter programs for the IBM PC that charac¬ 
terize optical waveguide modes. 8 In 
cooperation with Dr. Donnelly of MIT 
Lincoln Laboratory, we continued the work 
by Dr. L. Molter Orr, who was supported by 
the preceding National Science Foundation 
grant while working on her Ph.D. at MIT. 
This research has led to an improved under¬ 
standing of the crosstalk problem in wave¬ 
guide couplers. 9 


4 Y. Lai and H.A. Haus, ‘Quantum Theory of Solitons in Optical Fibers. I. Time-Dependent Hartree 
Approximations,' Phys. Rev. A 40:844-853 (1989); Y. Lai and H.A. Haus. ‘Quantum Theory of Solitons in 
Optical Fibers. II. Exact Solution," Phys. Rev. A 40:854-866 (1989); H.A. Haus and Y. Lai, 'Quantum Theory of 
Soliton Squeezing - A Linearized Approach," J. Opt. Soc. Am. B 7(3): 386-392 (1990). 

6 H.A. Haus, W.P. Huang, S. Kawakami, and N.A. Whitaker, 'Coupled-Mode Theory of Optical Waveguides," J. 
Lightwave Tech. LT-5:16-23 (1987). 

« W.P. Huang and H.A. H<.js, "Power Exchange in Grating-assisted Couplers," J. Lightwave Tech. 7:920-924 
(1989). 

i H.A. Haus and W.P. Huang, 'Mode Coupling in Tapered Structures," J. Lightwave Tech. 7:729-730 (1989); W.P. 
Huang and H.A. Haus, "Selfconsistent Vector Coupled-Mode Theory for Tapered Optical Waveguides," J. 
Lightwave Tech., forthcoming. 

« H.A. Haus, W.P. Huang, and N.M. Whitaker, "Optical Waveguide Dispersion Characteristics from the Scalar Wave 
Equation," J. Lightwave Tech. LT-5 1748-1754 (1987). 

» J.P. Donnelly, L A. Molter, and H.A. Haus, 'The Extinction Ratio in Opticsl Two-Guide Coupler A/? switches," J. 
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All-optical switching and all-optical logic 
gates require that the output of gates and 
switches are reasonably undistorted versions 
of input pulses. We have demonstrated an 
all-optical switch with a record low 
switching intensity of 0.5 W and with a very 
low distortion of the switched out pulse. 10 
We achieved the low switchout power by 
using a long propagation distance through 
an optical fiber. The switch, invented in our 
laboratory by a visiting Japanese researcher, 
is interferometric, using a single fiber and 
two orthogonally polarized versions of the 
pulse that interfere with each other. One of 
the two pulses is acted upon by a control 
pulse. Because only a single fiber is used, 
the interferometer is stable even when con¬ 
siderably long (i.e., 100 m). Modifications of 
the switch into integrable form will be fea¬ 
sible, when sufficient nonlinearities have 
been realized in semiconductor waveguides. 

1.3 Nonlinear Dynamics in 
Active Semiconductor 
Waveguides 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation 
Grant EET 88-15834 

U.S. Air Force - Office of Scientific Research 
Contract F49620-88-C-0089 

Project Staff 

Katherine L. Hall, Charles T. Hultgren, Professor 
Erich P. Ippen 

We have continued our study of carrier 
dynamics and nonlinear optical interactions 
in active semiconductor waveguide devices. 
Using 100-fs-duration pulses in the 
800-900 nm regime obtained by fiber com¬ 


pression of synch-pumped dye laser pulses 
and similar pulses in the 1.45-1,65-/im band 
from APM F-center lasers, we are investi¬ 
gating nonlinear dynamic behavior in both 
GaAIAs and InGaAsP devices under various 
excitation conditions. Varying the wave¬ 
length of the pump and probe beams, as well 
as injection current in our diode waveguide 
structures, we are studying interactions in the 
presence of gain, loss, or nonlinear transpar¬ 
ency. In all of these cases, an injected carrier 
density on the order of 10 18 /cm 3 makes the 
nonlinear optical behavior considerably dif¬ 
ferent from that observed in passive devices. 

In both GaAIAs and InGaAsP devices, we 
have discovered a strong nonlinearity due to 
nonequilibrium carrier heating. Heating of 
the carrier gas with respect to the lattice has 
a recovery time on the order of 1 ps in 
GaAIAs and 650 fs in InGaAsP. Since 
heating occurs via free electron absorption 
without a change in carrier number, recovery 
is complete. During the past year, we have 
documented gain and loss dynamics, which 
include two-photon absorption and hole¬ 
burning effects in addition to heating in a 
variety of different devices. In our present 
experiments, we are observing the index of 
refraction changes associated with these 
nonlinearities and evaluating their potential 
for use in all-optical switching. 

During 1989, we have also demonstrated a 
novel means for detecting nonlinear optical 
interactions in the active regions of diodes by 
monitoring changes in diode voltage. The 
forward bias voltage is measured as a func¬ 
tion of time delay between two pulses 
passing through the active region. The result 
is a nonlinear autocorrelation measurement 
that, when compared with optical pump- 
probe data, allows us to distinguish between 
those interactions involving carriers in the 
active region and other interactions. 


Quant. Electron. 25:924-932 (1989); J.P. Donnelly, H.A. Haus, and L.A. Molter, "Cross Power and Crosstalk in 
Waveguide Couplers," J. Lightwave Tech. 6:257-268 (1988). 

io M.J. LaGasse, D. Liu-Wong, J.G. Fujimoto, and H.A. Haus, 'Ultrafast Switching with a Single-Fiber Interfer¬ 
ometer," Opt. Lett. 14:311-313 (1989). 
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1.4 Saturation Characteristics 
of Semiconductor Optical 
Amplifiers 
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Although in the future optical amplifiers will 
be essential components of optical communi¬ 
cation networks and integrated photonic cir¬ 
cuits, very little is known about their 
nonlinear behavior. We need to improve our 
understanding of this behavior to properly 
evaluate the gain saturation and intermodu¬ 
lation distortion effects limiting device per¬ 
formance. In our laboratory, we have studied 
amplifier characteristics as a function of 
pulse duration and energy in the nonlinear 
regime. We have also developed an analyt¬ 
ical model to describe the observed behavior. 

In our experiments, we have been the first to 
observe pulsewidth-dependent gain in diode 
laser amplifiers. We have obtained data by 
measuring amplifier gain (output pulse 
energy/input pulse energy) as a function of 
input pulse energy. At low input energies, 
the amplification of ultrashort (100 fs) pulses 
is the same as that for longer (20 ps) pulses. 
As energy increases, however, the gain for 
the 100-fs pulses decreases more rapidly 
than that for the 20-ps pulses. This effect, 
observed at 1.5 n<r\ in both conventional and 
quantum-well quaternary devices, results in a 
3-db roll-off for the 100-fs pulses at an 
output energy that is about a factor of ten 
lower than that for the longer pulses. 

We have modeled this behavior with a 
system of coupled differential equations that 
includes changes in gain due to nonequilib¬ 
rium carrier distributions. The standard set of 
two equations describing the coupling 
between photon density and excited state 
carrier density is supplemented by a third, 
and in some cases a fourth, equation that 


represents the faster nonequilibrium dynamic. 
The time constant of this dynamic determines 
the pulsewidth at which gain roll-off occurs. 
The parameters of the model are consistent 
with results from time-resolved pump-probe 
experiments performed on the same devices. 
We have obtained excellent agreement with 
the observed saturation data. This model 
predicts that a more rapid sequence of fem¬ 
tosecond pulses should be ab'e to extract full 
energy. Although individual ultrashort pulses 
are restricted to lower gains, gain recovery 
between them is faster. 

1.5 Femtosecond Studies of 
Metallic and High Tc 
Superconductors 

Sponsors 
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Femtosecond pump-probe observations of 
reflectivity dynamics on metal surfaces 
provide a unique measure of nonequilibrium 
electron dynamics. Because of the high 
electron density, incident optical energy is 
absorbed directly by the electron gas. The 
resulting rise in electron temperature 
produces a dynamic change in reflectivity. 
Relaxation of this change occurs as the elec¬ 
trons lose energy to the lattice via phonon 
emission. The rate is governed by the 
electron-phonon coupling strength. 

We have begun to systematically study the 
strength of electron-phonon coupling, an 
important component in BCS theory of 
superconductivity, in collaboration with Pro¬ 
fessor M. Dresselhaus' group. In a series of 
experiments, we measured A{a> 2 } for ten dif¬ 
ferent metals (four of them supercon¬ 
ducting). The agreement between the values 
obtained with those described in the litera¬ 
ture was excellent. In addition, our method 
for measuring A{a> 2 }, the temperature relaxa¬ 
tion rate, has several advantages over other 
techniques such as tunneling or heat 
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capacity measurements because it: (1) is a 
direct measurement, (2) works at room tem¬ 
perature, and (3) can be applied to non¬ 
superconducting as well as superconducting 
samples. We have also demonstrated that 
thin overlayers of Cu, which has d-band 
transitions in the visible, can enhance the 
experimental reflectivity changes on metals 
which are otherwise too small to detect. This 
enhancement extends the method for use 
with any metal film. 

We have also performed a series of pump- 
probe and transmission experiments on three 
high-Tc thin films: YBa 2 Cu 3 0 7 .i, 

Bi 2 Sr 2 CaCu 2 0 8+ ,,, and Bi 2 Sr 2 Ca 2 Cu 30 io+y • In 
these materials, too, we have observed a 
strong correlation between observed relaxa¬ 
tion rates and T c . At present, there is no 
clear theoretical explanation for these results. 
We are performing additional experiments to 
improve our understanding of the observed 
correlation. 

1.6 Suppressed and Enhanced 
Spontaneous Emission from 
Microcavities 

Sponsor 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
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Optical microcavities could prove useful in 
the construction of efficient, high speed 
semiconductor lasers. One particularly inter¬ 
esting possibility is using a cavity to alter the 
spontaneous emission rate of the device. 
This alteration has been observed with atoms 
and molecules. However, altering the emis¬ 
sion rate in a semiconductor device is more 
difficult because the broad spontaneous 
emission bandwidth requires cavity dimen¬ 
sions on the order of a wavelength. To 
determine its potential application in semi¬ 
conductor devices, we have analyzed the 
radiation modes of oscillating dipoles in 
planar (one-dimensional confinement) and 


optical-wire (two-dimensional confinement) 
structures. When we used a mode counting 
method equivalent to, but less complex than, 
traditional classical and quantum electro¬ 
dynamical approaches, we found that an ide¬ 
alized planar metallic mirror cavity could 
suppress the spontaneous emission with 
respect to free space by no more than a 
factor of two. Emission was suppressed to a 
lesser degree with use of a real dielectric 
stack. Theory shows, however, that we 
could achieve greater suppression by 
restricting the dimensionality to that of the 
optical wire. Then, enhanced spontaneous 
emission should be easier to observe. With 
GaAs quantum-well cavities fabricated at 
NEC, Dr. Yokoyama observed reductions in 
the luminescence times by factors of two due 
to cavity effect. We have recently achieved 
femtosecond up-conversion detection of 
spontaneous emission from these structures 
and we are now attempting to observe more 
dramatic reductions with multiple-quantum- 
well structures grown at MIT by Professor 
Clifton G. Fonstad's group. 


1.6.1 JSEP Publications 

Ippen, E.P., H.A. Haus, and L.Y. Liu. “Addi¬ 
tive Pulse Mode Locking.” J. Opt. Soc. 
Amer. B 6(9): 1736-1745. 

Ippen, E.P., L.Y. Liu, and H.A. Haus. "Self¬ 
starting Condition for Additive-pulse 
Mode-locked Lasers.” Opt. Lett. 15(3): 
183-185 (1990). 

Kesler, M.P. and E.P. Ippen. "Femtosecond 
Time-Domain Measurements of Group 
Velocity Dispersion in AIGaAs Diode 
Lasers." Electron. Lett. 25(10): 640-641 
(1989). 

Nelson, K.A. and E.P. Ippen, "Femtosecond 
Coherent Spectroscopy." Advances in 
Chemical Physics. Vol. 75 (1989). 

Mark, J., L.Y. Liu, K.L. Hall, H.A. Haus, E.P. 
Ippen. "Femtosecond Pulse Generation in 
a Laser with a Nonlinear External 
Resonator.” Opt. Lett. 14(1): 48-50 
(1989). 
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1.7 New Ultrashort Pulse Laser 
Technology 
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1.7.1 Introduction 

Development of new ultrashort pulse laser 
technology is the cornerstone for studies of 
ultrafast phenomena and high speed optical 
measurement and communication tech¬ 
niques. In the area of ultrashort pulse gener¬ 
ation and measurement technology, we are 
developing new techniques to extend the 
range of applications and resolution of 
ultrafast studies. In particular, development 
of tunable wavelength femtosecond sources 
is important for a wide range of spectro¬ 
scopic as well as optoelectronic device 
studies. High repetition rate femtosecond 
amplifier technology permits the generation 
of high peak intensities for the study of non¬ 
linear optical effects. Finally, our research 
places special emphasis on developing new 
solid state laser materials and diode pumping 
techniques. These developments could rep¬ 
resent next generation technology, ultimately 
replacing existing dye laser systems to yield 
more compact and viable short pulse sources 
suitable for a wider variety of engineering 
applications. 


1.7.2 Tunable Femtosecond Ti:Al203 
Laser 

Solid state lasers provide an attractive alter¬ 
native to existing ultrashort pulse dye laser 
technology. Working in collaboration with 
researchers at MIT Lincoln Laboratory, we 
have developed a novel ultrashort pulse gen¬ 
eration modelocking technique in the 
Ti:Al 2 0 3 laser. 11 Ti:AI 2 0 3 is a promising new 
solid state laser material which features room 
temperature operation, high energy storage, 
and a broad tuning range of 670 nm to 1000 
nm. The broad gain bandwidth of this mate¬ 
rial make it ideal for femtosecond pulse gen¬ 
eration and amplification. In addition, the 
tuning range is especially useful for the study 
of GaAs-based electronic and optoelectronic 
materials and devices. 

Our ultrashort pulse generation technique is 
an extension of a pulse shortening technique 
developed in color center lasers known as 
additive pulse modelocking (APM). 12 
However, in the Ti:AI 2 0 3 laser, modelocking 
is achieved without the need for active gain 
or loss modulation. This self-starting addi¬ 
tive pulse modelocking results in a significant 
reduction in the system cost and complexity 
since high frequency active modulation is net 
required for the laser system. Extending this 
technique to other laser materials holds the 
promise of developing a compact and low 
cost ultrashort pulse laser technology suit¬ 
able for commercial application. 

The self-starting additive pulse modelocking 
technique uses an external cavity which con¬ 
tains an optically nonlinear medium (a length 
of optical fiber). The nonlinear external 
cavity functions resemble a nonlinear Fabry 
Perot. The reflectivity of the Fabry Perot is a 
function of the relative cavity lengths or 
phase of the optical field. Since the optical 
fiber has a nonlinear index of refraction, the 
reflectivity of the Fabry Perot is intensity 


11 J. Goodberlet, J. Wang, J.G. Fujimoto. and P.A. Schulz, 'Femtosecond Passively Mode-locked Ti:AI 2 0 3 Laser 
with a Nonlinear External Cavity," Opt. Lett. 14:1125-1127 (1989). 

12 P.N. Kean, X. Zhu, D.W. Crust, R.S. Grant, N. Langford, and W. Sibbett, "Enhanced Mode Locking of Color-Center 
Lasers," Opt. Lett. 14:39-41 (1989); K.J. Blow and D. Wood, ‘Mode-Locked Lasers with Nonlinear External Cav¬ 
ities," J. Opt. Soc. Am. B 5:629-632 (1988); J. Mark, L.Y. Liu, K.L. Hall, H.A. Haus, and E.P. Ippen, "Femto¬ 
second Pulse Generation in a Laser with a Nonlinear External Resonator," Opt Lett 14:48-50 (1989); E.P. Ippen, 
H.A. Haus, and L.Y. Liu, 'Additive Pulse Mode Locking," J. Opt. Soc. Am. B 6:1736-1745 (1989). 
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dependent. If the relative cavity lengths are 
controlled appropriately, the nonlinear cavity 
will function like a fast saturable absorber to 
enhance intensity fluctuations in the laser. 
This produces the pulse formation process. 

The additive pulse modelocked laser incorpo¬ 
rates an all-optical switch in the form of a 
nonlinear Fabry Perot which produces pulse 
shortening. Alternatively, this system can be 
viewed as the analog of an RF or microwave 
feedback system, except that the feedback is 
nonlinear, operating at optical frequencies. 
Figure 1 shows a schematic diagram of the 
Ti:AI 2 03 laser. In Ti:AI 2 0 3 , pulses as short as 
1.4 ps have been generated directly from the 
laser with pulse durations of 200 fs obtained 
by externally compensating for chirp. 

During the past year, we have studied the 
starting dynamics of the ultrashort pulse for¬ 
mation process and the scaling behavior of 
the laser system and performance. This 
research is relevant to extending the short 
pulse generation technique to other solid 
state laser materials at other wavelengths as 
well as to higher and lower power solid state 
lasers. 


1.7.3 Diode Pumped Modelocked 
Nd:YAG Lasers 

During preliminary studies in collaboration 
with investigators at MIT Lincoln Laboratory, 
we demonstrated the extension of this short 
pulse generation technique to a diode 
pumped Nd:YAG laser. 13 Other studies at RLE 
by Professors Ippen and Haus have achieved 
short pulse generation in a lamp pumped 
Nd:YLF laser. Results suggest that self¬ 
starting additive pulse modelocking can be 
applied to a range of solid state laser mate¬ 
rials. 

The diode pumped Nd:YAG laser combines 
the advantages of diode pumping which 


TUNABLE FEMTOSECOND Ti.AljOj LASER 



Figure 1. 

offers high inherent stability, compactness, 
and low cost with the technique of additive 
pulse modelocking which has a simple 
design and especially short pulse generation. 
We have achieved a pulse width of 1.7 ps, 
which is the shortest pulse generated to date 
in a modelocked Nd laser. Preliminary meas¬ 
urements show that this system has excellent 
interpulse arrival time or phase stability, 
making it a suitable source for optoelectronic 
diagnostic studies that require synchroniza¬ 
tion of optical and electronic signals. 

One of the key features of this technology is 
that the laser power and repetition rate can 
be scaled up or down. The extension of 
diode pumped self-starting additive pulse 
modelocking to lower output powers would 
permit development of extremely compact 
short pulse sources. Since the cavity length 
can be varied, this approach provides a 
method for obtaining variable pulse repetition 
rates. The development of compact short 
pulse sources could prove useful for a 
number of commercial applications including 
high speed optoelectronic sampling measure¬ 
ments and device diagnostics as well as pos¬ 
sible applications for high speed signal 
processing using optical techniques. 


’3 J. Goodberlet, J. Jacobson. J.G. Fujimoto, P.A. Schulz, and T.Y. Fan, "Self Starting Additive Pulse Modelocked 
Diode Pumped Nd:YAG Laser," Opt. Lett., forthcoming. 
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1.7.4 Multistage Femtosecond Dye 
Amplifier 

While solid state laser systems are an impor¬ 
tant area of investigation for the next gener¬ 
ation of ultrashort pulse laser technology, 
current measurement technology is based on 
dye laser systems and amplifiers. We are 
investigating new laser generation and high 
repetition rate amplification techniques for 
studying ultrafast phenomena, placing 
special emphasis on development of a multi¬ 
stage, high repetition rate, dye amplifier. The 
amplifier has a modular construction and may 
be easily configured for a variety of tasks. 

The ultrashort pulse source for our system is 
a colliding pulse modelocked ring dye laser 
(CPM). 14 This laser produces 50 fs pulses at 
620 nm wavelength. The wavelength is not 
tunable because the laser is passively 
modelocked. The output of the CPM is 
amplified using a copper vapor laser pumped 
dye amplifier. 16 The copper vapor laser oper¬ 
ates at a repetition rate of 8 kHz which facili¬ 
tates the use of signal averaging and lock-in 
detection techniques for high sensitivity 
measurements. 

We have developed a new multistage copper 
vapor laser pumped dye cell amplifier 
system. 16 The compact, multistage design 
allows the amplifier to be easily configured 
for advanced short pulse generation tech¬ 
niques. In comparison with earlier designs, 
we have improved its efficiency and exper¬ 
imental convenience while keeping the 
overall system complexity within manageable 
limits. 

The high pulse intensity available from just 
one stage of amplification (2—3 microjoules) 
permits use of nonlinear optical techniques 
to generate wavelength tunable pulses in the 
near infrared. After the first amplification 


stage, the femtosecond beam is focused in a 
cell of flowing water to generate a 
broadband continuum. A single wavelength 
may be selected from this continuum and 
amplified using an infrared dye in the second 
stage. 

We are also currently investigating a third 
configuration which will produce com¬ 
pressed amplified pulses with pulse durations 
shorter than 20 fs. In this arrangement, a 
short single mode optical fiber is placed after 
the first stage of amplification. This 
increases bandwidth of the pulses because of 
self-phase modulation. After a second stage 
of amplification, the pulses can be com¬ 
pressed with a dispersive delay line or phase 
filter. 

Thus, the high repetition rate amplifier forms 
the basis for a wide range of studies of 
ultrafast phenomena. Using nonlinear con¬ 
tinuum generation techniques, the system 
provides a wavelength tunable source. For 
very high time resolutions, we can apply 
pulse compression techniques, achieving 20 
fs measurement resolutions. 

1.8 Ultrafast Processes in 
Waveguide Devices 
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1.8.1 Optical Switching and 
Nonlinear Index in GaAs and MQW 
Waveguides 

Femtosecond investigations of bandedge 
nonlinear dynamics in guided wave devices 
are directly relevant to the development of 
high-speed all-optical switches. Our group 
has recently developed a novel technique for 
nonlinear index measurements. This tech¬ 
nique uses time division multiplexing to 
perform high resolution transient interfero¬ 
metric measurements on a nonlinear index. 
This approach reduces contributions from 
thermal and acoustic artifacts which tradi¬ 
tionally limit interferometric measurements. 
The system can detect phase shifts as small 
as /i/1000. 17 

For studies in semiconductor devices, we 
took measurements using a synchronously- 
pumped, tunable femtosecond dye laser 
system. Using the dye Styryl 9, we obtained 
430 fs pulses in the wavelength range 780 
nm to 880 nm. We performed systematic 
studies to characterize the nonlinear response 
in AIGaAs ridge waveguide devices. 18 We 
took the first direct measurement of the 
wavelength dependence of the instantaneous 
nonlinear index in this semiconductor. The 
results show values of n 2 in the range of 
10" 12 cm/W. Femtosecond time-resolved 
studies permit us to separate the real from 
the virtual contributions to the nonlinear 
index. 

In one area of our investigation, we study 
new materials for all-optical switching. 
Other materials we are currently investigating 
include multiple quantum well (MGW) GaAs 
waveguides. After performing preliminary 
measurements of femtosecond nonlinear 
effects, we are now in the process of system¬ 
atically characterizating the wavelength 
dependence of the nonlinear index as well as 
its linear optical properties. 

A second material system of interest to us is 
nonlinear optical polymers. We have recently 


begun studies in collaboration with investi¬ 
gators at Bell Communications Research to 
characterize the nonlinear index and absorp¬ 
tion behavior in polydiacetelene waveguides. 
Nonlinear optical polymers are a rapidly 
emerging research area, and the combination 
of Bellcore's unique materials growth and 
fabrication capabilities with our capabilities 
in femtosecond optics will permit state-of- 
the-art studies to be performed. 


1.8.2 Time Domain Optoelectronic 
Diagnostics 

Femtosecond dye laser systems as well as 
future generation Ti:AI 2 0 3 solid state laser 
systems can provide a tunable source of fem¬ 
tosecond pulses in the near IR wavelength 
range. This technology can form the basis 
for a wide range of diagnostics in optoelec¬ 
tronic devices. Diagnostics of optical wave¬ 
guide devices can be divided roughly into 
two catagories, measurement of nonlinear 
and linear properties. 

As discussed in the previous section, high 
peak intensities associated with short pulses 
can be applied to study nonlinear index or n 2 
effects in waveguide devices. In addition, 
other nonlinear effects such as two photon 
absorption can also be important in wave¬ 
guide devices. We have measured the two- 
photon absorption coefficient p in GaAs 
waveguide devices. Since two-photon 
absorption produces excited carrier popu¬ 
lations, this is an important limiting process 
for almost any application requiring high 
intensities in a waveguide. 

In addition to nonlinear properties, we can 
also develop time domain techniques to 
measure linear absorption, group velocity, 
and dispersion. Since femtosecond pulse 
durations have a spatial extent which is less 
than the cavity round-trip of most waveguide 
devices, the transient linear response can be 
measured. For example, if a short pulse is 


17 M.J LaGasse, K.K. Anderson, H.A. Haus, and J.G. Fujimoto, "Femtosecond All-Optical Switching in AIGaAs 
Waveguides Using a Single Arm Interferometer," Opt. Lett. 14:314 (1989). 

19 M.J. LaGasse. K.K. Anderson, C.A. Wang, H.A. Haus, and J.G. Fujimoto, 'Femtosecond Measurements of the 
Nonresonant Nonlinear Index in AIGaAs," Appt. Phys. Lett. 56:417 (1990). 
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injected into the waveguide, the pulse will be 
reflected internally within the device and 
several echoes will be emitted corresponding 
to multiple round-trip paths through the 
device. Using a cross-correlation technique, 
the time resolved behavior of the echoes can 
be measured. The decay of the echoes con¬ 
tains information about the facet reflectivity 
and absorption of the device. The delay time 
or time of flight between the echoes gives 
the group velocity dispersion. Finally, meas¬ 
urements performed as a function of wave¬ 
length permit us to characterize absorption 
versus wavelength and group velocity 
dispersion. We have used GaAs passive 
waveguide devices to demonstrate these 
techniques. 

By extending techniques further, we should 
be able to measure a wide range of optoelec¬ 
tronic device properties. With development 
of more compact and cost effective solid 
state tunable wavelength femtosecond lasers, 
a time domain optoelectronic diagnostic 
technology can become a viable alternative 
to continuous wave and other diagnostic 
approaches currently used for device charac¬ 
terization. 

1.9 Image Potential and 
Electron Dynamics in Metals 
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An improved understanding of electron 
dynamics in metals and at metal- 
semiconductor interfaces is necessary for the 


development of high speed electronic 
devices. Semiconductors and metals are the 
major constituents of devices. As device 
speeds increase with the advent of new nigh 
electron mobility and quantum transport 
effects in semiconductors, dynamic processes 
in metah will also become an important 
factor in affecting device performances. 
Since the density of electrons in metals is 
high, nonequilibrium processes and excited 
state dynamics occur extremely rapidly. 
Femtosecond optical measurement, with its 
high sensitivity and ultrafast temporal resol¬ 
ution, provide one of the few techniques for 
directly measuring dynamical processes in 
metals. 

We have continued our investigation of the 
transient dynamics of image potential states 
in metals. The image potential state has 
been a topic of recent active investigation in 
surface science. Image potential occurs 
when an electron, excited from the surface of 
a metal, is bound to its image charge in the 
bulk. 19 The electronic states form a Rydberg- 
like series somewhat analogous to a two- 
dimensional electron gas in semiconductor 
quantum wells. Because the electronic 
wavefunction is maximized outside the metal 
surface, scattering processes are reduced, 
and the electron relaxes by tunneling from 
the image potential state to bulk states. 
Thus, the image potential state provides an 
important model system for investigating 
electron dynamics in metals. 

Working in collaboration with researchers 
from General Motors Research Laboratories, 
our previous investigations combined femto¬ 
second pump-probe spectroscopy with 
photoemission techniques to directly 
measure, for the first time, the lifetime of an 
image potential state in a metal. 20 In our 
studies, ultraviolet femtosecond pulses are 
used to populate the image potential states. 
The image potential dynamics are probed by 
photoemitting electrons out of these levels 
with a visible femtosecond pulse delayed in 


19 D. Straub and F.J. Himpsel, "Identification of Image-Potential Surface States on Metals," Phys. Rev. Lett. 
52:1922-1924 (1984). 

20 R.W. Schoenlein, J.G Fujimoto, G.L. Eesley, and T.W. Capehart, 'Femtosecond Studies of Image-Potential 
Dynamics in Metals," Phys. Rev. Lett. 61:2596 (1988). 
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time. The energy spectrum of the 
photoemitted electrons is measured with a 
cylindrical mirror analyzer. This approach 
has allowed us to map the population of 
image potential states as a function of both 
time and electron energy. Studies performed 
on the the surface of Ag(100) revealed an 
n = 1 image potential lifetime of 15-35 fs, 
representing one of the highest time resol¬ 
ution transient photoemission measurements 
ever made. 

To test theoretical models for image potential 
electron dynamics, it is important to system¬ 
atically study different states in the Rydberg 
series using different metal surfaces. For the 
image potential state to be long-lived, the 
image potential energy must occur in a range 
where there is a gap of bulk states. The 
image potential electron must be confined in 
the potential barrier outside the metal from 
the crystal lattice by Bragg reflection. In this 
case, the image potential state relaxes by ine¬ 
lastic scattering. Conversely, if a different 
metal surface is examined where the image 
potential is resonant with allowed bulk 
states, then the dominant relaxation channel 
is elastic scattering with electronic states in 
the bulk. 

We have continued our earlier measuring 
efforts to examine the n = 1 and n = 2 image 
potential dyanmics in Ag(100) and 
Ag(111). 21 These experimental results can 
help explain differences in theoretical models 
based on simple single partical wavefunction 
analysis versus many-body effects. The 
image potential state has well defined 
quantum properties, providing an attractive 
model system for studying electronic 
quantum effects in metals. 


k 1.10 Laser Medicine 
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1.10.1 Heterodyne Ranging ft 
Opthalmic diagnostics ft 


The ability o 
corneal tissue i 
led to interest 
surgery. In this 
refractive power 
surgically to achieve refractive 
person's vision witlrout glassj 
the potential for m\ron-pr; 


excimer lasers t* remove 
submicron increments has 
lasers in kerttorefractive 
chnique, the^urvature or 
if the corn« 


is altered 
irrection of a 
To exploit 
:ision ablation 


control available with%jser sflrgery requires a 
precise method for nlbnitamng the incision 
depth. Previously, we McUraeveloped a fem¬ 
tosecond optical rangin^jjfchnique for meas¬ 
urement of corneal inalion depth. 22 This 
technique is analogous tjlthose used in radar 
or ultrasound except th|J s^ort pulses of light 
are used to perform precisameasurements of 
distance. Using thisr teclroique in corneal 
and intra-ocular surgjery, distances in the eye 
can be measured ncminvasivtjy, i.e., without 
contact with the Jye. Althaigh this tech¬ 
nique has achievjp high senswvity measure¬ 
ment perform^&e, it required use of 
ultrashort pulswfemtosecond lasers, and thus 
has limited pojpntial for clinical applications. 

Our recent rjkearch has focused fcn develop¬ 
ment of anralternate technique ernploying a 
short c^ierent length lignk source 
(AR-co^ed laser diode) and optitel hetero¬ 
dyne ejection to perform optica ranging. 
This npthod achieves measurement wvith lon- 
gitucftial and transverse resolutionsl(10 pm 


21 R.W. Schoenlein, J.G. Fujimoto, G.L. Eesley, and T.W. Capehart, "Femtosecond Image-Potential Dynamics in 
Metals: n = 2." Phys Rev. B 41:5436-5439 (1990). 

22 D. Stern, W Z. Lin, C.A. Puliafito, and J.G. Fujimoto, "Femtosecond Optical Ranging of Corneal Incision Depth," 
tnv. Ophthal. Vis. Sci. 30:99-104 (1989). 
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1.1 Introduction 
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Correlated fluctuations play an important role 
in systems with electronic and structural 
degrees of freedom. This role is most ubiqui¬ 
tous at phase transitions, but is also consid¬ 
erable away from phase transitions, 
extending down to the lowest temperatures 
due to quantum mechanics. The renormal¬ 
ization-group method is a new calculational 
method that can systematically deal with cor¬ 
related fluctuations at each successive scale 
length. Since we can include even the con¬ 
sequences of defects, for the first time we 
can obtain predictive microscopic theories for 
realistic systems. 

1.2 Finite-Temperature Phase 
Diagram of Vicinal Si(100) 
Surfaces 

Project Staff 

Professor A. Nihat Berker 

With the collaboration of Professor John D. 
Joannopoulos and Dr. Oscar L. Alerhand, we 
have combined electronic energy calculations 
and such statistical mechanics to obtain ab 
initio descriptions of finite-temperature semi¬ 


conductor surfaces and interfaces. The 
entropy, free energy, and other properties 
have been evaluated for the silicon (100) 
surface. The single-step/double-step phase 
diagram in the variables of crystal cut angle 
and temperature, as well as other observable 
properties such as step profiles, are predicted 
in very good agreement with ongoing exper¬ 
iments. Contrary to previous suggestions 
that only double-layer steps should appear 
on the equilibrium surface, it is predicted that 
the single-layer stepped surface is at equilib¬ 
rium for small misorientation a’ gles. This 
structure is stabilized by strain relaxation and 
by the thermal roughening of the steps. For 
annealed surfaces, the critical angle at which 
the transition between the single- and 
double-layered stepped surface occurs is cal¬ 
culated to be 9 c x2°. 

1.3 Absence of First-Order 
Phase Transitions in Physical 
Surface Systems 

Project Staff 

Professor A. Nihat Berker, Kenneth Hui 

Most recently, we made a theoretical predic¬ 
tion using the renormalization-group method 
that appears to have general and far-reaching 
consequences: we discovered that even an 
infinitesimal amount of randomness in inter¬ 
actions (e.g., distribution of defects) in 
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surface systems, converts first-order phase 
transitions, characterized by discontinuities, 
to second-order phase transitions, character¬ 
ized by infinite reponse functions. In bulk 
systems, as (calculable) threshold random¬ 
ness is needed for this conversion to occur. 
This general prediction appears to be sup¬ 
ported by experiments on doped KMnF 3 . 

1.4 New Orderings in Systems 
with Competing Interactions 

Project Staff 

Professor A. Nihat Berker, William Hoston, Roland 
Netz 

Our studies of realistic, complex systems with 
competing interactions have led to several 
new results. We have recently developed a 
new method that blends Monte Carlo simu¬ 
lation and mean-field theory. We are able to 
distinguish, for the first time, the effect of 
dimensionality on frustrated magnetic 
systems. We have obtained two ordered 
phases that nevertheless have considerable 
entropy. Also, we have recently obtained 
novel phases and multicritical points in 
systems with competing dipolar and 
quadrupolar interactions. 
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2.1 Introduction 
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In this research program, we use modern 
x-ray scattering techniques to study struc¬ 
tures and phase transitions in thin films and 
on surfaces. We have two principal exper¬ 
imental facilities, one at MIT and the other at 
the National Synchrotron Light Source at 
Brookhaven National Laboratory. At MIT, 
we have four high-resolution computer- 
controlled x-ray spectrometers using high 
intensity rotating anode x-ray generators. 
The angular resolution can be made as fine 
as 1.8 seconds of arc, which enables us to 
probe the development of order from dis¬ 
tances of the order of the x-ray wavelength, 
~lA, up to 30,000 A. The sample temper¬ 
ature can be varied between 2 K and 500 K 
with a relative accuracy of 2 x 10~ 3 K. At the 
National Synchrotron Light Source at 
Brookhaven National Laboratory, we have, in 
collaboration with IBM, three fully instru¬ 
mented beam lines. Two of these beam lines 
allow studies with photons varying in energy 
between 3 and 12 keV; the third has a fixed 
energy of 18 keV. These facilities make pos¬ 
sible high resolution scattering experiments 
with a flux more than three orders of magni¬ 
tude larger than that from a rotating anode 
x-ray generator, opening up a new gener¬ 
ation of experiments. 

As part of this JSEP program, we have built 
an x-ray compatible high vacuum single 
crystal apparatus. This enables us to use 
synchrotron radiation to study the structures 


and transitions occuring at a single surface, 
and, indeed, such experiments are now 
becoming routine. Our current experiments 
in this program are concentrated on the 
phases and phase transitions of metal and 
semiconductor surfaces and surface over¬ 
layers. 


2.2 Metal Surface Studies 

We have carried out detailed studies of the 
reconstruction of the (110) surface of gold: 
Au(110). Many noble metal facets are 
known to favor reconstructed structures at 
low temperatures. Further, at some temper¬ 
ature T c a reconstructed surface will typically 
undergo a reversible “deconstruction" to a 
high temperature structure which is no 
longer reconstructed. One such transition 
which has been studied theoretically as well 
as with electron diffraction is the Au(110) 
1x2 to 1x1 deconstruction. While the 
"missing row" model of the reconstructed 
1x2 surface is well established, the nature 
of the deconstruction transition itself remains 
controversial. 

Previous work by our group has demon¬ 
strated that x-ray scattering is an effective 
means of studying surface phase transitions. 
Specifically, the lineshapes of bulk-forbidden 
surface peaks provide information about the 
surface height-height (or step-step) corre¬ 
lation functions and are thus a sensitive 
probe of the surface roughening transition. 
This transition is characterized by a prolifer¬ 
ation of atomic steps which results in loga¬ 
rithmically divergent height fluctuations of 
the crystalline surface. We have carried out 


in 
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glancing angle synchrotron x-ray scattering 
experiments as a means of studying the 
Au(110) structures and deconstruction tran¬ 
sition. The results have been quite remark¬ 
able. The so-called 2x1 "missing row" 
model discussed above is a special case of a 
more general model in which one has (111) 
microfacets on the Au(110) surface. More 
extended microfacets can lead to recon¬ 
structions with longer periods: 1x3, 1x4, 
1x5, etc., and these are very close in energy 
to the 1x2. 

In our first set of experiments, the stable 
phase of the specific Au(110) surface turned 
out to be 1 x 3 rather than 1x2. The basic 
structure and the possible step excitations are 
shown in figure 1. For this structure super¬ 
lattice peaks occur at (0 0 if 3) positions. 
We find that at 485°C the surface undergoes 
a continuous, reversible transition. At this 
transition the 1x3 long range order is lost 
and the local period becomes incommensu¬ 
rate. This behavior results from the prolifer¬ 
ation of atomic steps (such as those shown 
in figure 1b and 1c), which act as domain 
walls between coherent regions of the recon¬ 
structed surface. In fact, by further studies 
we can determine that both types of steps 
are present. This means that the 1x3 
surface both deconstructs and roughens 
simultaneously. 

After further treatment of this same Au 
crystal, we were able to obtain a stable 1 x 2 
phase on the (110) surface. This latter struc¬ 
ture is the one most commonly observed in 
previous experiments. The deconstructive 
transition of this surface appears to be fun¬ 
damentally different from that of the 1x3 
phase. First, the structure remains 
commensurate even in the disordered phase. 
Second, there is no evidence of roughening. 
Third, we have some indication of 2D Ising 
critical behavior. Further analysis and meas¬ 
urements are required to elucidate the 1x2 
deconstruction fully. 


2.3 Semiconductor $ur*^e 
Studies 

We have investigated the surface disordering 
transition of Ge (111) near the bu'k melting 
point. The experiment '' c.s performed at the 
National Synchrotron Light Source beamline 
X20A during the summer of 1989. Our 
observation was made on Ge wafers 
resistively heated inside a surface chamber 
which has a base pressure of 1.3 x 10~ 10 torr. 
The sample surface was prepared by sput¬ 
tering for 30 minutes and subsequent 
annealing at the same temperature for one 
hour. Auger-electron spectroscopy was done 
in situ to check the surface cleanliness. The 
sample and surface chamber were oriented 
such that the in-plane (i.e., grazing inci¬ 
dence) (10) and (20) surface peaks, as well 
as the (11) bulk peak, were observable. 

We have studied the evolution of the (10) 
and (20) surface peaks between 700 K and 
1170 K, finding radical changes in intensity. 
The surface peaks remain resolution limited 



Figure 1. (a) Schematic representation of the 

Au(110) 1 x 3 missing-row reconstruction. The arrows 
indicate the southward atomic distrotions observed in 
the second and third rows from the surface, (b) A 
single atomic step thermal excitation of the 1 x 3 
surface, (c) A paried-step thermal excitation of the 1 x 
3 surface. 
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in the radial direction and the transverse 
widths stay the same throughout the entire 
temperature range. As the temperature goes 
up from 700 K, the intensity of the (10) 
surface peak decreases exponentially with a 
Debye-Waller factor of 0.001 K _1 . At about 
1000 K the rate of decrease of the (10) peak 
intensity per degree K intensifies. The peak 
intensity goes through a minimum at about 
1100 K and then recovers partially. The tran¬ 
sition shown here is reversible with temper¬ 
ature, indicating that it is indeed an 
equilibrium phase transition. Similar 
behavior is observed for the (20) surface 
peak. We have also looked for signs of a 
surface fluid layer which should manifest 
itself as a broad ring in reciprocal space, but 
found no indication of such a layer. Subse¬ 
quent to the x-ray scattering experiment, we 
have performed mass spectroscopy measure¬ 
ment on some Ge(111) samples. Each 
sample was placed one inch away from and 
with the normal pointing towards the probe. 
Our data show that there is an exodus of Ge 
atoms from the surface starting at about 
1050 K. 

None of the existing modeis of surface disor¬ 
dering — surface roughening and surface 
melting — can explain the data. The partial 
recovery of the surface peak intensity at 
about 100 K below bulk melting is partic¬ 
ularly intriguing. We propose a model which 
would give a surface peak intensity that is 
consistent with our observation. A Ge crystal 
with an ideally terminated (111) surface con¬ 
sists of bilayers of atoms, stacked in an 
..ABCABC.. sequence. The structure factor 
is identically zero at the surface peak position 
in reciprocal space because the scattering 
from every three bilayers of atoms exactly 
cancel. However, the penetration depth of 
x-rays in the direction normal to the surface 
is drastically reduced at grazing incidence, 
due to absorption and total internal 
reflection. The magnitude of the wave 
amplitude caused by a layer is maximum at 
the surface, dropping quickly as one moves 
into the bulk. As a result the sum of wave 
amplitudes from every three bilayers is 
nonzero: the amplitude from the top layers 
dominates that from below. Suppose that at 
1000 K atoms on the sample surface pop-out 
randomly, leaving vacant sites behind. This 
reduces the amplitude from the top layer of 


atoms. The net scattering intensity from the 
top three bilayers of Ge atoms reduces due to 
more effective cancellation. At 1100 K, the 
cancellation is optimized. At higher temper¬ 
ature, the random vacancy at the surface 
layer may be so large that the amplitude from 
the bottom layers now dominates. Since 
intensity is proportional to the modulus 
squared of the total wave amplitude from 
every layer, the intensity starts to increase 
once the random vacancy exceeds some crit¬ 
ical value. In this model there is no broad¬ 
ening of the surface peaks, which is 
consistent with our observation. The mass 
spectroscopy measurement also supports this 
picture. Of course, this model is by no 
means unique. Further studies are required. 


2.4 Rare Gases in Graphite 

We have studied the phases and the phase 
transitions for a xenon solid layer on graphite 
in the monolayer coverage regime. The 
incommensurate solid undergoes the 
sequence of transitions: aligned -» rotated -» 
reentrant aligned, before the solid transforms 
into the commensurate solid. The current 
theory of orientational epitaxy correctly pre¬ 
dicts only the point at which the rotation to 
the reentrant aligned phases occurs. None of 
the current theories can explain the whole 
picture of orientational epitaxy. Neither the 
high temperature aligned solid nor the reduc¬ 
tion of the rotation angle from that of the 
static prediction are explained. 

The reentrant aligned phase at low temper¬ 
atures is an incommensurate solid with a 
network of superlight domain walls. The 
observed first-order phase transition can be 
explained within the domain-wall theory. 
The incommensurability at constant coverage 
decreases approximately as [(T — T 0 )/T 0 ] 1/3 
until the first order C-IC transition point is 
hit. This one-third power law is close to the 
behavior seen for Kr on graphite following a 
similar thermodynamic pathway. 

We estimated the potential corrugation of the 
surface adsorption from the intensity ratios 
between the main and satellite peaks in the 
domain-wall incommensurate phase. The 
corrugation of the potential estimated was 
near the one predicted by Steele. Along with 
the estimation from the onset of the rotated 
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phase, we set the limits for the dimensionless 
wall width, 36 < < 42. 

At low temperatures the xenon layer on 
vermicular graphite did not transform into the 
commensurate phase as on the single-crystal 
substrate. Instead, the xenon layer seems to 
have a pinned domain-wall network due to 
defects. 
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Understanding the properties of surfaces of 
solids and the interactions of atoms and mol¬ 
ecules with surfaces is extremely important 
from both a technological and academic 
point of view. The recent advent of ultrahigh 
vacuum technology has made microscopic 
studies of well-characterized surface systems 
possible. The way atoms move to reduce 
surface energy, the number of layers of atoms 
involved in this reduction, the electronic and 
vibrational states that result from this move¬ 
ment, and the final symmetry of the surface 
layer are all important in achieving a funda¬ 
mental and microscopic understanding of the 
nature of clean surfaces, chemisorption proc¬ 
esses, and the initial stages of interface for¬ 
mation. 

The theoretical problems associated with 
surface systems are quite complex. We are 
currently, however, beginning to solve for the 
properties of these real surfaces systems 
(rather than for simple mathematical models). 
In particular, we continue to pursue our goal 
of calculating the total ground-state energy 
of a surface system from "first principles," to 
provide accurate theoretical predictions of 
surface geometries. In this program, we are 
concentrating our efforts on the areas of 
surface growth, surface reconstruction geom¬ 
etries, structural phase transitions, and chem¬ 
isorption. 


3.2 Thermal Amplitudes of 
Surface Atoms 

The vibrational excitations of semiconductor 
surfaces provide a direct probe of the struc¬ 
tural and bonding properties that result from 
their reconstruction. Moreover, studies of the 
vibrational properties of surfaces are neces¬ 
sary to understand the dynamical processes 
that occur on these surfaces (such as kinetics 
of adsorbates, diffusion and desorption, 
epitaxial growth, and possibly surface 
melting). The thermal fluctuations of the 
surface atoms are directly related to these 
phenomena. Furthermore, the analysis of a 
number of experimental probes requires a 
detailed knowledge of the vibrational ampli¬ 
tudes and correlations of the surface atoms. 
The Debye-Waller factors, which appear in 
the expressions for the intensity of radiation 
scattered from a surface, and the displace¬ 
ment-displacement correlations of the surface 
atoms enter into the interpretation of ion 
scattering and other experiments. 

New developments in a number of exper¬ 
imental techniques for studying vibrational 
properties of surfaces have generated further 
interest, the need for realistic calculations, 
and better theoretical understanding. These 
include He-atom scattering, high-resolution 
electron-energy-loss spectroscopy, surface 
infrared spectroscopy, medium- and high- 
energy ion scattering, and surface extended 
x-ray-absorption fine-structure experiments. 
From a theoretical point of view, an under¬ 
lying formalism for surfaces exists within the 
general framework of the theory of lattice 
dynamics, but we have made a few realistic 
calculations for specific surfaces. 

In this work, we have made calculations of 
the atomic displacement-displacement corre¬ 
lation functions for atoms on the Si(111) 
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2x1 and Si(001) 2x1 surfaces. Both of 
these surfaces have been the subject of 
intense research for many years, and even 
though we have made considerable progress, 
several questions regarding their structure 
and basic excitations remain unresolved. 
Using a semi-empirical tight-binding theory 
for structural energies in Si, we investigate 
the thermal fluctuations of the atoms on 
these surfaces. Specifically, we determine 
the anisotropy of the vibrational amplitudes 
of the surface atoms between vibrations par¬ 
allel and normal to the surface plane, their 
temperature dependence, and the interatomic 
vibrational correlations. 

The Si(111) 2x1 surface is currenlty under¬ 
stood in terms of a rr-bonding chain model. 
According to this model, the surface recon¬ 
structs by forming zigzag chains of surface 
atoms that run parallel to each other. The 
bonding along the chain is provided by the n 
interaction between the surface electrons. 
Figure 1 shows a side view of the Si(111) 
2x1 surface. The chain model also helps us 
understand the unique vibrational properties 
observed on this surface. Its lattice dynamics 
are dominated by the chain geometry and the 
strong coupling between the electronic 
surface states and the structural degrees of 
freedom of the zigzag chains. In fact, the 
most salient features in the surface phonon 
spectrum are related to modes that have the 
same form and elementary vibrational 
excitations of the surface zigzag chains. 

The appearance of a small-wavelength 
surface phonon with an unusually low 
energy is particularly relevant to the thermal 
fluctuations of the surface atoms. Indeed, 
one of the main reasons for this study is to 
investigate the thermal effects of this mode. 
This intrinsic vibration of the surface is 
dispersionless across the surface Brillouin 
zone, having an energy of fecu = 10meV. 
Near the zone edge, it is located well below 
the acoustic continuum. Because of the low 
energy of this phonon and its lack of 
dispersion (leading to a large contribution to 
the phonon density of states at that energy), 
this mode is thermally populated with a large 
probability. Thus, one might expect that the 
vibrational amplitudes of the surface atoms 
on Si(111) 2x1 would be greatly enhanced. 



Figure 1. Side view of the n-bonded chain model of 
Si(111) 2x1. The two surface atoms, numbered 1 
and 2, form zigzag chains that run parallel to the [112] 
direction (into the page); these chains are tilted with up 
(2) and down (1) atoms. 

In contrast with the result for Si(111) 2x1, 
this low-energy small-wavelength surface 
phonon does not appear in the phonon spec¬ 
trum of Si(001) 2x1. The basic building 
block in the reconstruction of Si (001) 
involves the formation of surface dimers 
between two neighboring surface atoms, 
halving the number of dangling bonds on the 
surface. Figure 2 shows a side view of the 
Si (001) 2x1 surface. The surface dimers do 
not remain symmetrical, tilting so that each 
dimer is formed by an up and down atom. 
With this model, we can understand the dif¬ 
ferent periodicities observed on the surface 
[2x1, p(2 x 2) and c(4 x 2)] in terms of dif¬ 
ferent ordered arrangements of the tilted 
surface dimers. We find that the normal 
modes of vibration of the higher-order recon¬ 
structions p(2 x 2) and c(4 x 2) are essen¬ 
tially the same as for the 2x1 
reconstruction, since the surface phonons are 
essentially characterized by the vibrations of 
the nearly isolated surface dimers. Conse¬ 
quently, we concentrate only on the 2x1 
reconstruction. 
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Figure 2. Side view of the tilted dimer model of 
Si(001)2 x 1. 


In the harmonic approximation, the corre¬ 
lation function between the displacement u iv 
of atom i along the direction n (from its equi¬ 
librium position), and the displacement of 
atom u ( * of atom j along the direction v, can 
be expressed in terms of the phonons or 
normal modes of the system: 




2N(M,M ) ) 1/2 


I 

q. * 


e iu (q, A)e jv (q, A) 


^(q) 


x e iq,R ' - R i'{i + 2n[cu,(q)]}, (1) 

where the sum is over normal modes x and 
wave vectors q in the Brillouin zone, cu^q) 
e(q, A) are the phonon frequencies and 
eigenvectors, respectively, Ri denotes the 
two-dimensional lattice vector of the surface 
unit cell where atom i is located, is its 
mass, and N is the number of unit cells 
within the periodic boundary conditions used 
in the calculations. The temperature T enters 
through the Bose-Einstein distribution func¬ 
tion n[<y]. 


The main difficulty in evaluating eq. (1) is 
that it requires the knowledge of the phonon 
energies and eigenvectors of the system over 
the entire Brillouin zone. Therefore, the cal¬ 
culation of vibrational correlation functions is 
usually limited to simple systems where the 
normal modes are determined by symmetry. 
Also, eq. (1) is usually approximated by low- 
and high-temperature expansions using 
models for the dispersion of phonon energies 
and the density of states, like the Debye 
model. 

To make reliable predictions and to extract 
microscopic information from the thermal 
vibrational amplitudes of the system, eq. (1) 
must accurately represent the system’s 
phonon energies and eigenvectors. For 
semiconductors surfaces, this representation 
is further complicated because their recon¬ 
struction leads to bonding configurations 
that can be quite different from the bulk; 
thus, dynamical properties cannot be simply 
described using force constants calculated 
for the bulk and transferred to the surface. 

A realistic theory must be used to calculate 
phonon energies and eigenvectors for each 
particular reconstructed surface. In this work, 
we use a formalism based on a semi- 
empirical tight-binding theory for structural 
energies in semiconductors. To explain the 
model briefly, the total energy of the system 
is expressed in terms of the band-structure 
energy, an empirical potential energy of 
interionic interactions, and the Coulomb 
energy of on-site electronic repulsions. The 
band-structure energy is calculated using an 
sp 3 tight binding representation of the 
valence electrons, and the interionic inter¬ 
actions are approximated with a sum of 
pairwise potentials. The on-site Coulomb 
repulsion term for the electrons is added to 
incorporate the effects of charge transfer. 
With this model, the dynamical matrix of a 
system can be calculated using perturbation 
theory, in which the bare ion-ion forces are 
renormalized by the electron polarizability; 
exponentially decaying interactomic forces 
are generated by the exchange of virtual 
electron-hole pairs between lattice distor¬ 
tions at different sites on the lattice. 
Because this approach incorporates the 
effects of reconstruction into the lattice 
dynamics of the surface, it successfully pro- 
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vides a description of the vibrational proper¬ 
ties of Si surfaces. 

The results of our calculations are presented 
in figures 3 and 4. Here we notice that for 
both Si(111) 2x1 and Si(001) 2x1, the 
surface atoms have thermal vibrational ampli¬ 
tudes that are larger along the direction per¬ 
pendicular to the surface plane than parallel 
to it. This is usually expected, since the 
surface atoms have reduced coordination 
normal to the surface. However, we find that 
the degree of anisotropy between normal and 
in-plane vibrations is veiy different for the 
two surfaces. Indeed, for Si(111) 2 x 1 we 
find that this anisotropy is quite large (1.8 in 
the mean-squared amplitudes of vibration), 
while on Si(001) 2x1 the fluctuations of 
the surface atoms are nearly isotropic. Fur¬ 
thermore, in the case of Si(111) 2x1, the 
anisotropy increases rapidly with increasing 
temperature, a feature we attribute to the 
presence of a low-energy dispersionless 
surface phonon. This large temperature 
enhancement is not observed for the surface 
atoms of Si(001) 2x1, which does not 
support a similar low-energy surface phonon. 
In comparing the atomic vibrational ampli¬ 
tudes on the two surfaces, the most salient 
feature is that, in the case of Si(111) 2x1, 
the surface atoms have very large amplitudes 
of vibrations normal to the surface. The dif¬ 
ference in the vibrational amplitudes of 
surface atoms on Si(111) 2x1 and Si(001) 
2x1 shows that the specific microscopic 
properties of the surface must be considered 
in the calculation of the atomic vibrational 
correlation functions. 

Specifically, we predict the degree of 
anisotropy between normal and in-plane 
vibrations and the difference in vibrational 
amplitudes between inequivalent surface 
atoms. Experimental techniques exist that 
can directly measure these variables. 
Perhaps the most important contribution of 
the calculations we have presented here is in 
their use in the analysis of experiments. Our 
calculated atomic vibrational amplitude and 
correlation function results might greatly 
reduce the uncertainties in the interpretation 
of scattering experiments that probe the 
structure of the surface; Debye-Waller factors 
might also be easily calculated with the 
information we have provided. 



Figure 3. Atomic rms vibrational amplitudes at T= 
270K as a function of penetration from the surface for 
(a) Si(111) 2x1 and (b) Si(100) 2x1. The 
Cartesian directions are defined in figures 1 and 2, 
respectively. 

3.3 Heteroepitaxial Growth 

Heterostructures of semiconductors, ele¬ 
mental or compound, have recently been the 
subject of intense study as technologically 
promising artificial materials. GaAs on Si is a 
system of particular interest, not only from a 
technological but also from a theoretical 
point of view. Usually GaAs growth is initi¬ 
ated on Si substrates oriented along the 
(100) direction. The quality of the GaAs 
overlayers is far from ideal with a very high 
density of defects and dislocations present, 
presumably due to lattice mismatch and 
polarity differences. However, it is not known 
to what degree other interface factors, such 
as surface reactions or substrate morphology, 
play an important role in growth inhibition or 
initiation. These issues are at the core of 
current growth models of GaAs on Si (100). 

In this work, we investigate, from an ab initio 
theoretical point of view, the possibility of 
layered growth of GaAs on atomically flat 
regions (terraces) of the Si surface. We shall 
focus on the initial stages of growth by con¬ 
sidering very low coverages (one to two 
monolayers) where a number of important 
issues can be isolated and addressed in 
detail. For example, what is the structure of 
the initial monolayers? Are the Ga and As 
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Figure 4. Atomic rms vibrational amplitudes of surface 
atoms on (a) Si(111) 2x1 (surface atom 2) and (b) 
Si(001) 2 x 1 (surface atom 1) as a function of tem¬ 
perature. The bold line is for a bulk atom. 

atoms locally mixed or do they separate into 
domains of each species? Is the Si surface 
covered by Ga and As atoms to the largest 
possible extent at monolayers thickness or do 
these atoms form GaAs bilayer islands on 
part of the surface? What is the structure of 
the initial bilayer of GaAs? Does it resemble 
that of bulk GaAs with alternating places of 
Ga and As or does it have a novel structure? 

To address the possibility of two-dimensional 
growth on Si terraces, we consider a slab of 
Si bounded by (100) surfaces on which we 
place Ga and As atoms. The issue of interest 
is the lowest-energy local structure and its 
effects on prospects for growth. There are 
two possibilities: the lowest-energy local 
structure may be the same as the structure of 
bulk GaAs, so that two-dimensional growth 
can proceed, or the structure may be dif¬ 
ferent, in which case layered growth of 
zincblende GaAs will be inhibited. The 
easiest way to envision formation of the 
lowest-energy structure is to consider a Si 
surface covered by a monolayer of one 
species (corresponding to pre-deposition, for 
example, of Ga or As) and to introduce some 
amount of the opposite species. One must 
then consider the possible local atomic con¬ 


figurations and identify the lowest-energy 
structure. In this sense, the exact relative 
abundance (chemical potential) of Ga and 
As atoms is not important, but the presence 
of both species is crucial. 

Each surface layer of the slab is taken to be 
reconstructed, for simplicity, in the standard 
(2x1) dimer pattern (this choice will be jus¬ 
tified below). The Ga and As atoms are 
placed on the slab in five different ways as 
shown in figures 5(a)-5(e). The total ener¬ 
gies resulting from ab initio calculations 
allowing these geometries to fully relax, are 
presented in table 1. 



trrr^=-H!awBH!!3 am _!■ smn 



Figure 5. Schematic representation of overlayer 
arrangements. (a) Pure monolayer coverage, 
(Ga)[Si)(As). (b) Mixed monolayer coverage, (GaAs) 
[Si](AsGa). (c) Bilayer coverage with pure layers and 
As-exposed surface, (As)(Ga)[Si]. (d) Same as (c) 
with Ga-exposed surface, (Ga)(As)[Si]. (e) Bilayer 
coverage with mixed layers, (GaAs)(AsGa) [Si]. 
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Figure 

Configuration 

Relative 

Energy 

Difference 

5(a) 

(Ga)[Si](As) 

0.8 

5(b) 

(GaAs)[Si](AsGa) 

0.0 

5(c) 

(As) (Ga) [Si J 

2.2 

5(d) 

(Ga) (As)[Si] 

2.8 

5(e) 

(GaAs)(AsGa)[Si] 

1.5 


Table 1 . Relative energy difference of the various over¬ 
layer configurations as depicted in figure 5 in eV/(2X1) 

sue. 


Based on the energies of the different con¬ 
figurations, a number of important conclu¬ 
sions can be drawn: 

1. When the available amount of Ga and As 
atoms is less than a full monolayer, under 
equilibrium stoichiometric conditions, it 
wets the Si surface and formation of 
thicker GaAs islands is inhibited. 

2. At monolayer coverage, the two species 
of atoms do not separate into pure Ga or 
pure As domains, but are locally mixed 
on the Si surface. 

3. When GaAs is available in stoichiometric 
amounts sufficient for bilayer coverage, 
the equilibrium structure of the bilayers is 
not a succession of pure Ga and As 
layers as in bulk GaAs along the (100) 
direction. Rather, chemical and rehybrid¬ 


ization reactions are energetically very 
favorable, leading to a passivating mixed 
overlayer configuration. Once this mixed 
layer is formed, it could persist in a disor¬ 
dered phase even under non- 
stoichiometric conditions because of its 
considerable stability. The mixed layer 
resembles locally a wurtzite structure 
under large uniaxial stress. A structural 
barrier (requiring creation of extensive 
defects at large energy cost) is encount¬ 
ered in reverting from the mixed layer 
structure to the zincblende structure of 
bulk GaAs. 

This conclusion has important implications: 
Since the mixed layers have very different 
structure from the thermodynamically stable 
zincblende phase of bulk GaAs, the growth 
of mixed layers on terraces will be under¬ 
mined by any surface topology which could 
lead to nucleation of zincblende GaAs. Thus, 
nucleation centers of zincblende GaAs at 
specific surface topologies will prevail in the 
growth process. It is possible to grow bulk 
GaAs experimentally on vicinal Si(100) sur¬ 
faces, that is, surfaces with large density of 
steps. It is likely that surface steps are pre¬ 
cisely the topology that nucleates zincblende 
GaAs, and growth occurs directly on the 
steps. Indeed, experimental evidence 
appears to support this view. Search for a 
possible microscopic mechanism for the initi¬ 
ation of this kind of growth, which involves 
formation of zincblende GaAs seeds at 
double-layer steps on the Si surface, is cur¬ 
rently in progress. 
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4.1 Project Description 
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In the late 1970s, the rapidly improving the¬ 
oretical description of disordered metals 
combined with the rapidly evolving tech¬ 
nology of submicron lithography to create a 
new subfield of condensed-matter physics. 
The electronic properties of artificial struc¬ 
tures — which are large compared with 
atomic size scales but comparable to 
quantum mechanical coherence lengths — 
have been shown to be fundamentally dif¬ 
ferent from both those of larger and smaller 
structures. For this reason, the term 
“mesoscopic" was coined to differentiate 
them. 

Almost from the start, the hallmark of 
mesoscopic structures has been the random, 
noise-like, but reproducible fluctuations in 


1 IBM Thomas J. Watson Research Laboratories. 


their conductance. These fluctuations were 
first seen in the insulating regime of one¬ 
dimensional field-effect transistors when 
their Fermi energy was varied. The most dra¬ 
matic discovery was that, in metallic samples, 
independent of their size and geometry, the 
fluctuations had magnitude of order e 2 /h in 
the limit of zero temperature. These Uni¬ 
versal Conductance Fluctuations were seen 
when the magnetic field, Fermi energy, or 
impurity configuration was varied, consistent 
with the theory of Altshuler, Lee and Stone. 2 

Recently, we reported the discovery of a new 
phenomenon in mesoscopic field effect 
structures. In the same regime where expo¬ 
nentially large random fluctuations had been 
seen in earlier devices, 2 we found fluctu¬ 
ations which, while still exponentially large, 
were accurately periodic in the density of the 
electron gas. We were the first to observe 
these oscillations in Si MOSFETs 3 fabricated 
at MIT and in novel GaAs inverted 
semiconductor-insulator-semiconductor field 


* For a review see M. A. Kastner, R. F. Kwasnick, J. C. Licini, and D. J. Bishop, 'Conductance Fluctuations near the 
Localized-to-Extended Transition in Narrow Si MOSFETs," Phys. Rev. B 36: 8015 (1987). 

3 J.F.H. Scott-Thomas, M.A. Kastner, D A. Antoniadis, H.l. Smith, and S B. Field, 'Si Metal-Oxide Semiconductor 
Field Effect Transistor with 70-nm Slotted Gates for Study of Quasi-One-Dimensional Quantum Transport," J. 
Vac. Sci. Technology B 6: 1841 (1988); J.F.H. Scott-Thomas, S.B. Field, M.A. Kastner, H.l. Smith, and D.A. 
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effect transistors ISIS structures fabricated at 
IBM. 4 An example of the latter is shown in 
figure 1. 

A rich phenomenology is associated with this 
oscillatory behavior. The conductance is 
thermally activated between 1 K and about 
0.3 K, and the activation energy oscillates 
with the number of electrons per unit length. 
At a lower temperature there is a tunneling 
component which also oscillates. The con¬ 
ductance is highly non-ohmic, showing dra¬ 
matic threshold behavior with threshold 
voltage that oscillates out of phase with the 
conductance. Most surprising of all, the 
oscillations are independent of magnetic 
field. This is probably why the phenomenon 
was not observed earlier by other groups that 
typically study MODFETs in which the carrier 
density can not be varied over a wide range. 
For the quantum Hall effect, the variation of 
carrier density is equivalent to the variation 
of magnetic field. This is not the case for our 
new phenomenon. 

The phenomenology suggests that the oscil¬ 
lations result from the opening of an energy 
gap whenever an integral number of elec¬ 
trons resides in a subsection of the one¬ 
dimensional channel. Several models have 
been proposed to explain the observations. 
It is probable that the electron-electron inter¬ 
actions are crucial. If so, this will be the first 
time that many-body effects will have been 
shown to dominate the behavior of 
mesoscopic systems. 

The silicon device fabrication and part of the 
GaAs fabrication was sponsored by the 
National Science Foundation under Grant 


ECS 88-13250. A portion of the ultralow 
temperature measurements was sponsored by 
the Joint Services Electronics Program under 
Contract DAAL03-89-C-0001. 



Figure 1. (a) Conductance vs. gate voltage of a 

narrow channel, 2 ftw long, measured at T = 50 mK. 
The inset shows an expansion of the first few oscil¬ 
lations. (b) Fourier power spectrum of the above. The 
peak corresponds to a period of AV»1.7 mV. 


Antoniadis, "Conductance Oscillations Periodic in the Density of a One-Dimensional Electron Gas," Phys. Rev. 
Lett 62: 583 (1989). 

4 U. Meirav, M.A. Kastner, M. Heiblum, and S.J. Wind, ‘A One-Dimensional Electron Gas in GaAs: Periodic 
Conductance-Oscillations as a Function of Density," Phys Rev. (Rapid Comm.) B 40: 5871 (1989); U. Meirav, 
M. Heiblum, and F. Stern, Appl. Phys. Lett. b2: 1268 (1988). 
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5.1 Project Description 
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We are investigating the conductance of 
narrow MOSFETs that are subject to a strong 
magnetic field. This investigation was moti¬ 
vated by an experiment conducted by Pro¬ 
fessor Marc A. Kastner's group two years 
ago, 2 in which the group studied the trans¬ 
port properties of narrow MOSFETs 
(width ~ 500A) subject to a large magnetic 
field of 5T to 10T. In these fields, only the 
lowest few Landau levels — in the quantum 
Hall regime — are occupied. At the same 
time, the cyclotron orbit is comparable to the 
sample width, making the finite size effect of 
the sample important. The experimentalists 
observed thresholds for large increases in 
conductance as a function of magnetic field, 
as well as ill-defined steps in the conduc¬ 
tance of approximately 2e 2 /h. We tackled this 
problem using both computer simulation and 
analytic treatment of many-body effects. 3 


We use computer simulation to model disor¬ 
dered conductors subject to a large magnetic 
field. The model is defined on a tight 
binding lattice of finite width and length 
which is extended to infinity in both 
directions by connecting it to ideal lattices 
without randomness to simulate the source- 
drain leads. The novel feature of the model 
is that the leads are also subject to the same 
quantizing magnetic fields so that (1) the 
edge states in the leads consist of Landau 
level subject to a lateral boundary condition, 
and (2) the states at the Fermi level are edge 
states which are bound to the edge and 
propagate to the right along the upper edge 
and left along the lower edge. The conduc¬ 
tance is then related to the transmission 
probability of these edge states using the 
Landauer formula. These transmission prob¬ 
abilities are computed numerically. 

When the width of the wire is large com¬ 
pared with the cyclotron radius, the edge 
states are basically forward scattered in the 
disorder region so that the transmission 
probability is unity, leading to quantized Hall 
conductance of e 2 /h. It is only when the 
width becomes comparable to the cyclotron 


1 Department of Physics, University of Munich, Germany. 

2 J. Scott-Thomas, S B. Field, M.A. Kastner, H.l Smith, and D.A. Antoniades, Phys. Rev. Lett. 62:583 (1989). 

* J. Kinerct and P.A. Lee, to be published. 
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radius that the edge state can scatter across 
the width of the sample and is then reflected. 
The reduction in the transmission probability 
leads to ordinary metallic conduction. As the 
field is increased, the simulation shows a 
threshold for the appearance of quantized 
steps in the conductance in agreement with 
the experiment. The simulation also shows 
that in the step between plateaus, the system 
behaves as a metal in the sense that it 
exhibits conductance fluctuations in the 
order of e 2 /h. The period of the fluctuation as 
given by the condition that the flux through 
the sample area is a flux quantum just as in 
universal conductance fluctuation. 

The next puzzle to be resolved is under¬ 
standing why the quantum Hall step is in 
units of 2e 2 /h. We believe that the 2e 2 /h step 
is due to the presence of two degenerate 
valleys in the silicon interface, which is not 


resolved. In ordinary two dimensional 
samples, the valley degeneracy is resolved 
due to the exchange effect. It is more ener¬ 
getically favorable to occupy one valley first, 
rather than occupying both valleys equally, 
because the exchange energy is nonlinear in 
the density. However, this mechanism will 
be weakened in narrow wires, because the 
Landau level is no longer completely flat and 
it costs kinetic energy to preferentially 
occupy one valley. We calculated the 
exchange energy for a model system con¬ 
fined by a parabolic potential demonstrating 
this competition between kinetic and 
exchange terms in quasi-one-dimensional 
MOSFETs. We conclude that the combina¬ 
tion of numerical simulation and analytic 
work to account for the many-body effect 
gives a quite satisfactory account of the 
transport properties of a narrow MOSFET 
wire in a quantizing magnetic field. 
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6.1 Project Description 
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Under the Joint Services Electronics 
Program, we have studied the Pt(001) 
surface as a model for the behavior of an 
incommensurate overlayer on a substrate, 1 
This research will lead to a more complete 
understanding of what determines the orien¬ 
tation of a thin film on a substrate. 

The (001) surface of platinum, like the (001) 
surfaces of gold and iridium, exhibits a 
radical reconstruction in which the topmost 
layer of atoms forms a hexagonal layer, in 
spite of the layers of square symmetry 
beneath. At T = 300 K, the hexagonal over¬ 
layer is rotated by 0.75 degrees away from a 
high symmetry direction of the cubic lattice 
(the substrate). Our recent synchrotron x-ray 
measurements have revealed, however, that 
the overlayer experiences a completely new 
class of rotational transition at high temper¬ 
atures. For T > 1670 K, the hexagonal over¬ 
layer is aligned with the cubic substrate. The 
temperature T c = 1670 K marks the onset of 
rotation, and as the temperature decreases 
further the rotation angle increases contin¬ 
uously, reaching 0.75° at T = 1300K. Figure 


1 shows typical experimental profiles. The 
two peaks are arranged symmetrically about 
the high symmetry [1, 1, 0] direction and 
correspond to positive and negative rotation 
angles of the overlayer. The accompanying 
change in the overlayer lattice constant is on 
the order of one part in 1000. 

Current theories of orientational epitaxy 
predict a definite relationship between the 
overlayer lattice constant and the rotation 
angle, 2 but our observations contradict this 
understanding. Nevertheless, the evolution 
of rotation angle with temperature follows a 
power-law behavior (see figure 2). In fact, 
the observed exponent is one-half, indicative 
of so-called mean-field behavior. An expo¬ 
nent of one-half is unusual for a two- 
dimensional phase transition because the 
effect of fluctuations are expected to be par¬ 
ticularly severe in two-dimensions. However, 
an elegant interpretation is possible: If one 
supposes that the rotational transition occurs 
as a result of a rotational elastic constant 
becoming zero, then by considering the 
effect of such behavior on the low-energy 
modes (lattice vibrations) of the Pt overlayer, 
it emerges that there is no mode that 
becomes soft as the transition is approached. 
In other words, there simply are no fluctu¬ 
ations, and character of the transition must 
be mean-field. 3 With this observation, we can 


1 D.L Abernathy, S.G.J. Mochrie, D M. Zehner, G. Griibel, and D. Gibbs, ‘Rotational Transitions of the Pt(001) 
Surface," unpublished.;D. Gibbs, G. Grubel, D.M. Zehner, D.L. Abernathy, and S.G.J. Mochrie, ‘Phase Behavior of 
the Pt(0G*i) Surface. X-Ray Measurements," unpublished. 

2 J.P. McTague and A.D. Novaco, Phys. Rev. B 19: 5299 (1979); J. Villain, Phys. Rev. Lett. 41:36 (1978); H. 
Shiba, J. Phys. Soc Jap. 46:1852 (1979), 48:211 (1980). 

3 R.A. Cowley, Phys. Rev. 8 13: 4877 (1976). 
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also understand the mean-field character of 
other high temperature rotational transitions. 4 


Pt(QOl) 



Figure 2. Evolution of the rotation angle vs. temper¬ 
ature for the Pt(001) surface. 


Figure 1. Transverse scans through a near hexagonal 
peak of the reconstructed Pt(001) surface for temper¬ 
atures decreasing from the aligned phase (single peak) 
into the rotated phase (two peaks). 


4 K.L D'Amico. D.E Moncton, E D. Specht, R.J. Birgeneau, S.E. Nagler and P.M. Horn, Phys. Rev. Lett. 53: 2250 
(1984). 
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Cavity quantum electrodynamics — the study 
of interactions of individual atoms with 
quantum fields in cavities — has grown into 
an active area of quantum optics. 2 A seminal 
experiment in this area was the inhibition of 


spontaneous emission of an excited atom, 
first demonstrated in our laboratory some 
years ago. 3 During this past year, we have 
completed a study of the evolution of small 
ensembles of atoms in a microwave cavity. 
Our goal is to observe the dynamics of an 
atom-vacuum system in the quantum regime. 

The states of our collective N-atom system 
are Dicke states. 4 * The total atom-field Hamil¬ 
tonian for the system is 6 

H t0 t = H A + H F + H af 

= fcct> 0 D° + ha>(a^ a + 'A) 

+ 'Ahco a (aO + a D ), (1) 


1 Department of Physics, Wellesley College, Wellesley, Massachusetts. 


2 S Haroche and D Kleppner, Phys. Today B 42(1): 24 (1989). 


3 R.G. Hulet, E S. Hilfer, and D. Kleppner, Phys. Rev Lett. 55: 2137 (1985). 

* R.H. Dicke, Phys. Rev B 93: 99 (1954). 


s S Haroche, in New Trends in Atomic Physics, eds G Grynber t , nd R Stora (Amsterdam North Holland, 1984) 
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where H A and H f are the atomic and field 
Hamiltonians, respectively, and H AF describes 
the atom-field interaction, a and a' are the 
conventional destruction and creation opera¬ 
tors for the electromagnetic field. D* and D 
are collective atomic raising and lowering 
operators, respectively, that act on the Dicke 
states like standard angular momentum 
raising and lowering operators. D° acts on 
the Dicke states like the z - component of 
the angular momentum operator. 

The evolution of the system can be described 
by the density matrix equation 

^ ~ C^tot' Pa+f] + ^F^A+F' (2) 


The last term in equation (2) accounts for 
dissipation. The probability that a member of 
the system is excited at time t is related to 
the expectation value of D° by 

P e (t,N) = < D° > "I ~ (3) 

For more than fifty atoms, an alternative 
method based on the Bloch vector model 
becomes more practical than the above 
approach. The time derivatives of the 
quantum mechanical operators D°, D\ and a' 
are found in the Heisenberg picture. The 
atomic operators are related to the compo¬ 
nents (r/, p) of the Bloch vector, a classical 
angular momentum vector with length 
J = N/2 

D° -► t] = %N cos 0 and D* -» p = sin 0e"\ 
where 0 and <f> are the standard angles in 
spherical coordinates. The field operator a' is 
identified with c, a classical electric field in 
the cavity. The evolution of the system is 
described by® 


d rj_ 
dt 





(4a) 


d/> 

dt 


= - '<"a 


(4b) 


—- = i<5c - ye + p + F^t) (4c) 
dt 2 

where 6 = o> — a> 0 Equation (4c) includes a 
damping term and a random force F'(t) that 
account for the coupling of the radiation 
mode to the thermal reservoir of the cavity 
walls. It can be shown that the probability 
of excitation at time t is given by 

P a (t.N) = -jL < »|(t,N, 0(0)) > 0(0) + -1- (5) 

where the average over 0(0) is carried out 
over an appropriate distribution of initial 
angles arising from thermal fluctuations. 

There is no general way to parametrize the 
behavior of the atom-cavity system, but for 
underdamped motion one feature is partic¬ 
ularly easy to observe experimentally: this is 
the time T m , r , at which P„(t,N) achieves its first 
minimum. We define the "collapse 
frequency" v col by 

v coi = (2t min )“ 1 (6) 


Although P„(t,N) is a complicated function of 
the system parameters, Haroche 5 has pointed 
out that the variation of v col with the number 
of atoms N can be accurately parameterized 
by 


v C0 .(N) = 


a x /N~ 

b + InN 


(7) 


where a and b are adjustable parameters that 
depend on a>„ y, and /?. The accuracy of this 
expression is illustrated by figure 1 which 
displays values of v co , obtained from the exact 
solution [equation (3)] for N < 50 and from 
the Bloch method [equation (5)] for 
N > 200, along with the best fit of equation 
(7) to the calculations. 

The experiments employ an atomic beam of 
calcium Rydberg atoms and a split supercon¬ 
ducting cavity operated at 35 GHz. At the 
ambient temperature of 2 K, the mean 
blackbody photon number is 0.8. Selective 


8 B J Hughey, Cavity Modified Atom-Photon Interaction, Ph D. diss., Dept of Physics, MIT, 1989 
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Figure 1. Dependence of the "collapse" frequency v co! 
on number of atoms N for Q = 1.2x10 7 and 
T = 2.17K The filled circles are obtained from the exact 
method and the open squares from the Bloch method, 
as discussed in the text. The solid curve is a fit of the 
calculations to equation (7), which yields 
a = 0.3187MHz and b = 2.245. 

field ionization allows us to monitor simul¬ 
taneously the populations of the initial and 
final states. The time evolution of the atomic 
system is probed by "Stark switching." The 
collective oscillations of energy between 
ensembles of atoms and a cavity with 
Q > 10 7 are being studied for one to several 
hundred atoms. 

The experiment is carried out with the appa¬ 
ratus shown in figure 2. A calcium atomic 
beam is prepared in the 46p state (abbrevi¬ 
ation for the 4s46p 1 P-, state) inside the 
cavity using a three-step pulsed dye laser 
system. The final laser beam is polarized 
along the cavity electric field so that the 
46p(m=0) state is excited relative to this 
quantization axis. The two states involved in 
the atom-cavity oscillations are thus the 
46p(m=0) (upper state) and 46s (lower 
state). The frequency of the 46p -* 46s tran¬ 
sition in calcium is 35.332 GHz; the atom- 
cavity coupling frequency a>J2n calculated 
for this transition and our cavity dimensions 
is 105 kHz. Observing underdamped 
behavior for a single atom in the cavity 
requires Q>2x10 5 . To avoid undesirable 
effects of blackbody radiation, the one-atom 
experiments require temperatures <2K, for 
which b<0.8. 

The atoms are excited inside the cavity. The 
mean transit time to the cavity exit is 18 



l£ on 


Figure 2. Diagram of the apparatus, side view. A: 
atomic beam, H: liquid helium temperature shield, N: 
liquid nitrogen temperature shield, M: single layer mu¬ 
metal shield, P: collimating plate of 2 channel field 
ionization detector, C: cavity, R: ramped plate, S: 
slotted plate, Ch: channel electron multipliers, W: wave¬ 
guide, L: laser beams, B: waveguide holder, Co: 
coupler. The tuning mechanism is not shown for 
clarity. 

/isec, which corresponds to two cycles of the 
single-atom-cavity oscillation. After the 
atoms exit the cavity, they enter the detector 
region where selective field ionization is used 
to detect and differentiate the two states. 
We measure P e (t, N) for each value of N by 
scanning the time t when a 200 mVjsulse is 
applied between the cavity halves. (N is the 
mean number of atoms in the cavity.) In this 
way, we count the number of atoms detected 
in the 46p and 46s states as a function of the 
atom-cavity interaction time. Sample data 
are shown in figure 3. 

An important feature of the experiment is the 
use of a split superconducting cavity. By 
applying an electric field between the two 
halves, the atomic resonance can be shifted 
far from resonance with the cavity, thereby 
terminating the atom-cavity interaction. This 
effectively freezes the atomic population at 
the moment at which the voltage pulse is 
applied. By using a carefully designed choke 
groove structure, we suppress leakage from 
the cavity and achieve a quality factor 
Q>10 7 in niobium and lead-plated copper 
cavities. 

A series of_time evolution curves for various 
values of N was studied. Figures 4 and 5 
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Figure 3. Observed atom-cavity oscillations for 
Q = 1.18x10 7 and T = 2.17K for two values of Nd- 
The solid curves are fits of the data to damped 
sinusoids, from which v CO i is extracted. The open 
circles are the atoms detected in the 46p state, and the 
filled squares are the atoms detected in the 46s state. 
Ti nt is the atom-cavity interaction time. 

display some of the data along with the cal¬ 
culated time evolution curves. 

We also compared our results with theory for 
a single atom in the cavity. Single atom 
experiments are of special significance 
because they represent the limit of small 
systems. Furthermore, they have the exper¬ 
imental advantage of being free of much of 
the uncertainty in detection efficiency, 
assuming that the mean number of atoms 
detected per pulse is small compared to 
unity. For these studies, we generally oper¬ 
ated with an average of 0.04 atoms detected 
per laser pulse. The observed and calculated 
evolution of a single atom in the cavity with 
Q= 1.5x10® and T=1.95K are shown in 
figure 6. We believe that the damping which 
is evident in the data is due to stray electric 
and magnetic fields, and possibly thermal 
effects. 

In support of these studies, we also carried 
out a high resolution study of the microwave 
spectrum of calcium Rydberg states. The 
results are described in a thesis, 7 and in the 
publications listed below. 



Figure 4. Comparison of calculated and observed 
evolution of the system for large ensembles of atoms in 
the cavity. The filled circles are the number of atoms 
detected in the 46p and 46s states as a function of 
atom-cavity interaction time Tj„ t . The amplitudes of the 
theory curves are scaled to compare with data but there 
are no other adjustable parameters. (a) and (b): 
N = 380. (c) and (d): N = 300. 

Publications 

Gentile, T.R., B.J. Hughey, T.W. Ducas, and 
D. Kleppner. "Experimental Study of 
Two-Photon Rabi Oscillations." Paper 
presented at the Sixth Conference on 
Coherence & Quantum Optics, Rochester, 
New York, June 1989. 

Gentile, T.R. Microwave Spectroscopy and 
Atom-Photon Interactions in Rydberg 
States of Calcium. Ph.D. diss. Dept, of 
Physics, MIT, 1989. 

Gentile, T.R., B.J. Hughey, T.W. Ducas, and 
D. Kleppner. "Microwave Spectroscopy 
of Calcium Rydberg States." Phys. Rev. A. 
Accepted for publication. 

Hughey, B.J., T.R. Gentile, T.W. Ducas, and 
D. Kleppner. "A Split High Q Supercon¬ 
ducting Cavity." Submitted to Rev. Sci. 
Instrum. 


7 T.R. Gentile, Microwave Spectroscopy and Atom-Photon Interactions in Rydberg States of Calcium. Ph D diss. 
Dept, of Physics, MIT, 1989. 
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polarizers. A scanned antenna-cyfpled 
pyro^ectric detector system has beejfcievel- 
oped \r mapping the millimeter beam 

profiles^ map of the radiation pj^mle at the 
atomic barn is shown in figure 

We have b«un designing the Controlled field 
region in wroh the actual measurement will 
take place Vid hope toVcomplete con¬ 
struction of th^apparatus Crthe coming year. 
Efforts are alsoLcontinu^j to optimize the 
hydrogen beam Ipurce^md the efficiency of 
the circular state otodu^ion. 


Fig^e 11. Intensity profile of the 256 GH^|diation 
u \jBa to drive the n = 29 -* n = 30 circular tranWion. 

2.4 Precision Mass 
Spectroscopy of Ions 

Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
National Science Foundation 
Contract PHY 86-05893 

Project Staff 

Kevin R. Boyce, Eric A. Cornell, Deborah Kuchnir, 
Professor David E. Pritchard 

In 1989, we made the first mass comparison 
of single trapped ions, raising the state of the 


art for fractional precision in mass measure¬ 
ment to 4 x 10 10 . This is a first step toward 
our ultimate goal of determining the mass of 
individual atomic and molecular ions with 
precisions of 10 " or better. This precision 
will give us the capability of making exper¬ 
iments that address issues in both funda¬ 
mental and applied physics: 

• The 3 H + — 3 He + mass difference is impor¬ 
tant in ongoing experiments to measure 
the electron neutrino rest mass. 

• Excitation and binding energies of typical 
atomic and molecular ions can be deter¬ 
mined by weighing the small decrease in 
energy: Am = E b md/c 2 . 

• Experiments that weigh y-rays can be 
used in a new method to determine 
Avogadro's number, N A , a key funda¬ 
mental constant whose accurate deter¬ 
mination would permit the replacement 
of the artifact mass standard by an 
atomic mass standard. 

• Traditional applications of mass spectros¬ 
copy should benefit from the several 
orders of magnitude improvement in both 
accuracy and sensitivity which our 
approach offers over conventional tech¬ 
niques. 

In our experimental approach, we measure 
ion cyclotron resonance on a single mole¬ 
cular or atomic ion in a Penning trap, a 
highly uniform magnetic field with axial con¬ 
finement provided by weaker electric fields. 
We monitor ions, oscillating along magnetic 
field lines, by the currents induced in the trap 
electrodes. Working with only a single ion is 
essential because space charge from other 
ions leads to undesired frequency shifts. 
This work in trapping and precision reso¬ 
nance draws on techniques developed by 
Hans Dehmelt at the University of 
Washington and Norman Ramsey at Harvard 
University, for which they shared the 1989 
Nobel Prize. 

Our most notable accomplishment this year 
was determining the carbon monoxide- 
molecular nitrogen mass ratio to an accuracy 
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of 4 x 10~ 10 . 13 The ratio of the cyclotron fre¬ 
quencies of the two ions is the inverse of the 
ratio of the masses as long as the magnetic 
field remains constant. By trapping first a 
single Nj ion and measuring its frequency, 
and then swapping to a single CO" ion, and 
then back again to Nj, we are able to correct 
for drifts in the magnetic field to about 0.5 
parts per billion (see figure 12). 

We have developed techniques for driving, 
cooling, and measuring the frequencies of all 
three normal modes of Penning trap motion. 
Thus, we can manipulate the ion position 
reproducibly to within 30 microns of the 
center of the trap, correcting for electrostatic 
shifts in the cyclotron frequency to great 
accuracy. We use a rr-pulse method to 
coherently swap the phase and action of the 
cyclotron with the axial modes.' 4 Therefore, 
although we detect only the axial motion 
directly, we can determine cyciotron fre¬ 
quency by measuring how much phase accu¬ 
mulates in the cyclotron motion in a known 
time interval (see figure 13). 

Currently, precision is limited by magnetic 
field imperfections and temporal instabilities. 
Achieving our long range goal of 10’ 11 preci¬ 
sion requires either dramatic improvements in 
field stability or simultaneous comparison of 
two different ions. With two ions of equal 
charge, the ion-ion perturbations are very 
similar for each ion and, hence, do not intro¬ 
duce significant uncertainty in the mass ratio. 
Our group has succeeded in trapping a single 
ion and a single Nj ion simultaneously. 15 
,(i the coming year, we plan to develop tech¬ 
niques for making precision resonance meas¬ 
urements on two ions simultaneously. 


Publications 

Cornell, E.A., R.M. Weisskoff, K.R. Boyce, 
R.W. Flanagan, G.P. Lafyatis, and D.E. 
Pritchard. "Single-Ion Cyclotron Reso¬ 
nance Measurements of M(CO + )/M(Nj)." 
Phys. Rev. Lett. 63:1674 (1989). 

Cornell, E.A., R.M. Weisskoff, K.R. Boyce, 
and D.E. Pritchard. "Mode Coupling in a 
Penning Trap: n Pulses and a Classical 
Avoided Crossing." Phys. Rev. A 41:312 
(1990). 

Kuchnir, D.L. Trapping and Detecting Two 
Different Single Ions at Once . A Step 
Towards Ultra-High-Precision Mass Com¬ 
parison Measurements. B.S. thesis. Dept, 
of Physics, MIT, 1989. 


2.4.1 Atom Interferometry 
Sponsors 

Joint Services Electronics Program 
Contract DAAL03-89-C-0001 
U.S. Army Research Office 

Contract DAAL03-89-K-0082 
U.S. Navy - Office of Naval Research 
Contract N00014-89-J-1207 

Project Staff 

Chris R. Ekstrom, David W. Keith, Bruce G. 
Oldaker, Quentin Turchette, Professor David E. 
Pritchard 

Using fabricated transmission gratings as 
optical elements for the matter waves, we are 
constructing an atom interferometer to phys¬ 
ically separate atom waves before recom¬ 
bining them. This interferometer will be 
useful in studies of atomic properties, tests of 
basic quantum physics, for metrology, as 
rotation sensors, and, perhaps, ultimately as 
devices to make ultra-small structures using 
atom holograms. 


13 E.A Cornell, R.M. Weisskoff et al., Phys. Rev. Lett. 63:1674 (1989). 

14 E.A. Cornell, R.M. Weisskoff et al., Phys. Rev. A 41:312 (1990). 

15 D.L Kuchnir, Trapping and Detecting Two Different Single Ions at Once: A Step Towards Ultra-High-Precision 
Mass Comparison Measurements, B.S. thesis, Dept, of Physics, MIT, 1989. 
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Figure 12. The ojta from a mass comparison run are shown. The solid squares (open squares) are the cyclotron 
frequency of Nj(C f*). A total of three ions were loaded in the order NJ - CO* - NJ. The solid lines are a fit to the 
two frequencies assuming a field drift that is linear in time. The dotted-line assumes a quadratic field drift. The 
indicated value for the difference frequency results from the latter assumption, and corresponds to M(CO + )/M(N£) 
=0.9995988876(' . 


During the last year, our atom interferometer 
evolved from a rough plan to an essentially 
complete devic At present, we have tested 
all its major components at least once. We 
will test the complete system during the 
coming year. 

Our interferometer consists of three 0.2 
/tm-period diffraction gratings equally spaced 
~ 0.65 m apart in our atomic beam machine. 
The maximum separation of the beams will 
be ~ 60jum. The first two gratings separate 
and redirect the atomic beam forming a 
standing wave interference pattern in the 
atomic flux at the third grating, which acts 
like a mask to sample this pattern. 

We can estimate our anticipated final signal 
strength from the properties of the individual 


gratings. Attenuation caused by the primary 
grating and the grating support structure 
gives an intensity in the 0 ,h order of ~ one- 
eighth of the incident intensity and one- 
sixteenth in each of the ± 1 orders. The final 
intensity detected at the maximum of a fringe 
after transmission through all three gratings 
can be calculated by summing the ampli¬ 
tudes for the two sides of the interferometer 
and will be 0.005 of the incident intensity. 
The fringe constant will be 4 to 1, so the 
interference signal will be 0.004 of the inci¬ 
dent intensity. We are anticipating that the 
final interference signal through the interfer¬ 
ometer will be ~ 4 x 10 3 counts per second. 
This should generously exceed the noise of 
the detector ( < 100 sec- 1 ). 
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Figure 13. For each plotted point, we perform the following experiment: The initially cold ion is pulsed into a 
cyclotron orbit of known initial phase and then allowed to evolve "in the dark’ for an indicated amount of time, t. 
Then a pulse is applied which exchanges cyclotron and axial motions, bringing the ion's cyclotron action and phase 
into the axial mode. As the ion's axial motion rings down, its phase is detected. The appropriate multiple of 360° is 
added, and a line is fitted to the points. The slope of the line is the frequency difference between the frequency 
generator and the trap cyclotron frequency. 


The mechanical vibrations of our machine are 
a principal technical obstacle because they 
could blur the interference pattern. There are 
two types of required limits on vibrations. 
First, the three gratings must move relative to 
each other by less than ~ 1/4 period (50 nm) 
during the time the final grating samples the 
intensity at a given position. Thus, the rms 
amplitude of relative vibrations integrated 
over all frequencies greater than the recip¬ 
rocal of the integration time must be less 
than ~ 50 nm. The second requirement is 
related to the motion of the gratings due to 
acceleration of, or rotation about, the center 
of mass of the grating system during the 1.3 


ms time it takes for the atoms to traverse the 
interferometer. This means that below ~ 900 
Hz the rms acceleration must be less than 
10 2 ms -2 , and the rms angular velocity must 
be less than 10 _s radians per second. 

We have solved our vibration problem by 
using a combination of passive isolation and 
active feedback. The passive isolation system 
consists of small pneumatic feet which act 
like damped springs to support the machine 
and give it a 3-Hz resonant frequency. This 
isolates the machine from building noise at 
higher frequencies. We have used the active 
feedback system to stabilize the relative posi- 
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Figure 14. Our current atom interferometer with laser interferometer vibration isolation system is shown. (Not to 


scale.) 


tions of the three gratings at frequencies 
below ~ 150 Hz. This system works best at 
low frequencies < 10 Hz) where the passive 
system is least effective. The reduction of 
relative motion provided by the active system 
will allow us to use much longer integration 
times when we are looking for the interfer¬ 
ence signal. The active feedback system uses 
a laser interferometer which has the same 
transmission grating geometry as the atom 
interferometer. We mounted the gratings for 
the optical interferometer on the same three 
translation stages as the matter wave gratings 
to record the exact relative orientation of the 
matter wave interferometer. We apply the 
signal from the optical interferometer to 
provide a measure of the relative alignment 
of the three grating platforms to a Piezo¬ 
electric translator (PZT) through a feedback 
network to stabilize the platforms (figure 14). 
With this system, we have adequately 
reduced the relative rms motion (at frequen¬ 
cies less than 0.3 Hz) of the gratings from ~ 
1500 to 40 nm. We also reduced the rms 
acceleration in a frequency range to which 
the interferometer is sensitive from 1.1 x 10~ 2 
to 2.3 x ICHms 2 , which is safe by a factor 


~ 5, ensuring sufficiently low angular veloci¬ 
ties of the apparatus. 

We have also made significant progress in 
overcoming another technical obstacle, the 
relative alignment of the atom gratings. For 
all points along the height (3 mm) of our 
ribbon shaped beam to have the same phase 
of interference signal, the gratings must be 
aligned to an angle of ~10~ 6 rads, with 
respect to rotations about the beam axis. We 
accomplished this by using a technique 
based on the optical polarizing properties of 
the gratings. 

In addition to the work on vibrations and 
alignment mentioned above, our main 
progress during 1989 was construction of 
the various mechanical components that 
position the gratings inside the vacuum 
envelope. Also, we have written computer 
software and built electronic hardware to 
control the position of the three grating plat¬ 
forms and the detector, and the height, angle 
and position of the second collimating slit. 
Instead of directly varying the voltage of the 
PZT that controls the position of the last 
grating, now we will use the computer to 
vary the null point of the active feedback 
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system when we are searching for interfer¬ 
ence fringes. 

When we have successfully demonstrated 
this interferometer, our first experimental 
objective will be a demonstration of Berry's 
phase with bosons. Another experiment 
could be an improved measurement of the 
Aharonov-Casher effect. 

Publications 

Keith, D.W., and D.E. Pritchard. "Atom 
Optics." In New Frontiers in QED and 
Quantumoptics. New York: Plenum Press. 
Forthcoming. 

Keith, D.W., M.L. Schattenburg, H.l. Smith, 
and D.E. Pritchard. "Diffraction Gratings 
from Atoms." QELS Conference, Balti¬ 
more, Maryland, April 1989. 

Oldaker, B.G., P.J. Martin, P.L. Gould, M. 
Xiao, and D.E. Pritchard. "Experimental 
Study of Sub-Poissonian Statistics in the 
Transfer of Momentum from Light to 
Atoms." Submitted to Phys. Rev. Lett. 

Pritchard, D.E., and B.G. Oldaker. "Light 
Forces and Atom Diffraction: An Illus¬ 
trated Summary." Paper presented at the 
Sixth Conference on Coherence, 
Rochester, New York, 1989. 

Pritchard, D.E. Experimental Studies of Atom 
Diffraction and the Mechanical Forces of 
Light on Atoms, NICOLS Conference. 
Bretton Woods, New Hampshire: Aca¬ 
demic Press, 1989. 


2.5 Neutral Atom Trap 

Sponsors 

U.S. Navy - Office of Naval Research 
Contracts N00014-83-K-0695 and 
N00014-89-J-1207 


Project Staff 

Kristian Helmerson, Michael A. Joffe, Dr. Min 
Xiao, Ke-Xun Sun, Professor David E. Pritchard 

We have trapped large numbers of neutral 
atoms, and cooled them to millikelvin tem¬ 
peratures. Our next objective is to cool them 
to microkelvin temperatures. Dense samples 
of atoms cooled to microkelvin temperatures 
promise to open up new and exciting areas 
of physics. The lack of interaction of the low 
velocity atoms due to their reduced thermal 
motion, together with the possibility of 
indefinitely long interaction times, make 
samples of trapped atoms ideal for high 
resolution spectroscopy and for use as 
atomic frequency standards. High density 
samples of ultra-cold atoms will also make 
possible new studies of interatomic collisions 
and collective effects, such as Bose 
condensation. We have made progress using 
our existing magnetic trap. In addition, we 
started a new project to develop a contin¬ 
uous source of slow atoms to load into 
future magnetic traps. 


i.1 Magnetic Trap for Neutral 
Atoms 


tral j 

gr rieut 


Now fltet techniques for trapping'neutral 
atoms a^ well established, 16 jbfi field of 
neutral albm trapping has /moved from 
infancy to adolescence andyfhe emphasis is 
now on doinVexperiment^with the trapped 
atoms. 

Currently, our flain effort in neutral atom 
trapping is coolmk/trapped atoms to low 
temperatures. Wh^Lthis remains a difficult 
and elusive goal/(t^^iate, we have only 
achieved micyefKelvin ^temperatures with 
untrapped aj*fms), the ^M/ards for super¬ 
cooling trajrfped atoms ar^eignificant. The 
long confinement times, together with the 
reduce*/thermal motion of co^atoms, could 
resujrin a new era of ultra-higb resolution 
spectroscopy and precise frequency stand- 
fas. Potentially more exciting is trie possi¬ 
bility of combining the high alnsities 
achievable in traps and the long deBrdalie 


19 D.E. Pritchard, K. Helmerson, and A.G. Martin, 'Atom Traps,' in Atomic Physics,. 11, eds. S. Haroche, J.C. Gay, 
and G. Grynberg (Singapore: World Scientific, 1989), pp. 179-97. 
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We are developing a generalized formulation 
for the electromagnetic fields in multilayered 
uniaxially anisotropic media containing arbi¬ 
trary distribution of current sources. A spec¬ 
tral domain dyadic Green's function for 
multilayered uniaxially anisotropic media 
containing three-dimensional sources is 
derived. Tractable forms are shown to be 
easily deduced from the physical picture of 
the waves radiated by the primary sources 
and the multiple reflections from the stratified 
medium. The formulation decomposes the 
dyadic Green's function into TE and TM 
waves. The dyadic Green's function in the 
source region is properly represented by 
extracting the delta function singularity. A 
simple procedure to obtain the fields in any 
arbitrary layer is described. Recursion 
relations for appropriately defined reflection 
and upward and downward propagating 
transmission coefficients are presented. 
Forms suitable for transmission line applica¬ 
tions in multilayered media are derived. 

Finite difference time domain (FDTD) tech¬ 
niques show great promise in their ability to 
solve three-dimensional problems with arbi¬ 
trary geometry. Advantages of this method 
include the ability to model spatially or 
temporally varying media. These advantages 
are due to the complete discretization of both 
space and time. Considering the volume of 
information being calculated, these tech¬ 
niques are very efficient and well suited to 
calculation on future parallel processing 
computers. This method was first formulated 
by Yee in 1966, and his basic algorithm is 
still in use. 

Recent work has demonstrated the applica¬ 
bility of the FDTD technique to microstrip 
problems. The centered finite difference 
approximations used are second order accu¬ 
rate in both space and time yielding good 
results for reasonable mesh sizes. Numerical 
techniques used to solve electromagnetic 
problems must limit the domain over which 
the fields are calculated. This mandates the 
use of an absorbing boundary condition to 
simulate the outward propagation of waves 
incident on the walls of the mesh. An 
absorbing boundary condition, developed by 
Mur and based on the work of Enquist and 
Majda, is used in this work. 


Our work in this area includes development 
of the algorithms mentioned above into a 
general purpose computer code which can 
be used to solve for the transient response of 
electromagnetic problems having arbitrary 
geometries. In addition to the transient 
response, frequency domain parameters can 
be obtained by Fourier transform of the time 
domain results. Since the fields are calcu¬ 
lated throughout space and time all other 
desired parameters can be calculated from 
the field quantities. 

Specifically, we are analyzing rectangular 
microstrip structures with as many as two or 
more ports. Such structures can be used in 
MMIC filters or antennas. This problem is of 
interest for several reasons. First, there are 
existing frequency domain solutions to the 
resonance problem of a rectangular micro¬ 
strip patch, which we compare with the 
FDTD solution. Secondly, the FDTD tech¬ 
nique can be used to analyze coupling of 
microstrip lines to the rectangular structure. 
This coupling can be either a direct con¬ 
nection or a gap coupled connection. Advan¬ 
tages of the FDTD solution to this problem 
are that (1) the full wave solution allows for 
radiation or surface wave loss, (2) no empir¬ 
ical values such as “effective” dimensions are 
needed for the analysis, and (3) the geom¬ 
etry can be altered easily to allow for various 
connections or coupling to the patch. This 
method is a significant improvement over 
those that rely on a planar circuit approach in 
which the substrate thickness must be small 
compared to wavelength and inherently 
three-dimensional coupling problems are not 
easily handled. We will compare our results 
with various planar circuit approaches. 

The leakage phenomenon is important in the 
area of millimeter-wave integrated circuits 
and integrated optics. We have performed 
theoretical analyses and experiments to 
investigate this phenomenon. The leakage is 
due to the TE-TM coupling occurring at the 
geometrical discontinuities, and the leaky 
power in the form of surface wave propa¬ 
gates in the background medium. 

We are using various methods to analyze the 
dielectric strip waveguides including the 
approximate field matching, effective dielec¬ 
tric constant (EDC), and mode matching. 
Because the first two methods are approxi- 
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mate, they can not be used to predict the 
imaginary part of the propagation constant. 
When using the third method, we must place 
ground planes at some distance from the 
guiding structure and omit the effect of radi¬ 
ation loss. 

To solve for the dispersion relation of single 
and coupled dielectric strip waveguides, we 
derive an integral equation formulation using 
dyadic Green's function. We also present a 
method for predicting the leakage. Three dif¬ 
ferent dielectric strip waveguides are investi¬ 
gated. optica! rib waveguide, strip dielectric 
guide, and insulated image guide. Both 
single and coupled strip waveguides are 
studied. The cross section of the dielectric 
strips are assumed to have rectangular shape. 
Applying the Galerkin's method, the field dis¬ 
tribution on the cross section is represented 
by a set of unit pulse basis functions. By 
substituting these basis functions into the 
integral equations and choosing the same set 
of basis functions as the testing functions, 
we can obtain a determinant equation from 
which the propagation constant can be 
solved. 

For single dielectric strip waveguide, we have 
observed that the leakage occurs when the 
effective refractive index is smaller than that 
of a surface wave mode in the background 
medium. We also observed that if the lowest 
TE-like (TM-like) mode is leaky, the lowest 
TM-like (TE-like) mode is non-leaky. When 
the lowest order mode leaks, the surface 
wave mode of opposite polarization is 
excited. When the higher order mode leaks, 
the surface wave modes of both polarizations 
can be excited. 

For two symmetrical dielectric strip wave¬ 
guides, we investigated both the even and 
odd modes. In the leaky mode, the total 
leakage is due to the leakage from each indi¬ 
vidual strip waveguide. At the separation 
where the even mode has a maximum 
leakage, it is implied that the surface wave 
modes excited by each waveguide add in 
phase. For the odd mode at about the same 
separation, these surface wave modes add 
out of phase; hence, a null in the leakage is 
observed 

We are analyzing the propagation properties 
of single and coupled inhomogeneous slab 


waveguides and propose an integral equation 
formulation using the dyadic Green's func¬ 
tion which covers both the TE and TM 
modes. The dispersion relations are obtained 
by applying Galerkin's method to solve the 
integral equation. We also investigate the 
coupling between two symmetrical inho¬ 
mogeneous slab waveguides. This method is 
applicable to arbitrary dielectric constant pro¬ 
files. 

Full modal analysis is used to study the 
dispersion characteristics of microstrip lines 
periodically loaded with crossing strips in a 
stratified uniaxially anisotropic medium. 
Dyadic Green's functions in the spectral 
domain for the multilayered medium in con¬ 
junction with the vector Fourier transform 
(VFT) are used to formulate a coupled set of 
vector integral equations for the current dis¬ 
tribution on the signal line and the crossing 
strips. Galerkin's procedure is applied to 
derive the eigenvalue equation for the propa¬ 
gation constant. We investigate the effect of 
anisotropy for both open and shielded struc¬ 
tures on the stopband properties. 

We are analyzing the quasistatic fields gener¬ 
ated by an electrode mounted on a perfectly 
conducting pad of finite extent buried in a 
planar stratified medium. An integral 
equation in the spectral domain is derived for 
the outflowing current density distribution on 
the pad-electrode surface. The method of 
moments is then applied to solve the integral 
equation. We are investigating the effects of 
the electric properties of the stratified 
medium and the standoff thickness on the 
total electrode current. Several conductivity 
profiles modelling different practical meas¬ 
urement environments are also considered. 
Numerical results reveal that the total 
electrode current is insensitive to the standoff 
thickness, and can be used to prospect for 
the conductivity of rock formation. 

Based on a hybrid transmission lines-lumped 
element circuit model, we analyze the tran¬ 
sient propagation characteristics of VLSI 
interconnects with discrete capacitive loads 
at various locations. Exact expressions of the 
Laplace transform of unit step responses are 
first obtained through the ABCD matrix for¬ 
mulation. We then apply the equivalent 
dominant pole approximation to the transfer 
function with the propagation delays fac- 
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tored out. The approximated transfer func¬ 
tion can be inverted in closed form and 
quickly evaluated. These results provide effi¬ 
cient ways of finding approximately the 
effects on delays and rise time brought by 
VLSI off-chip interconnects. 

Because of the dramatic increase in device 
densities on microelectronic chips, the prop¬ 
agation delay for off-chip interconnects has 
become the factor limiting the speed of VLSI 
packages. Typical scales of these intercon¬ 
nects will be comparable or larger than the 
characteristic wavelength of the high fre¬ 
quency components of the signal. Therefore, 
to calculate the delays caused by these inter¬ 
connects properly, a hybrid circuit model 
containing transmission line sections as well 
as lumped elements must be used in place of 
the all-lumped elements. Nevertheless, most 
circuit simulation packages are based on the 
latter and have to resort to subsection 
approximation when dealing with trans¬ 
mission lines. This scheme will undoubtedly 
lead to lengthy computation — which is not 
desirable when a quick, heuristic estimate of 
bounds is needed for the initial phase of the 
design cycle. 

We have developed two approaches for 
obtaining the approximate transient response 
without lengthy simulation. The first empha¬ 
sizes the calculation of bounds to voltage 
responses from the differential equations 
either by direct integration or by using the 
optimal control theory. The second analyzes 
the properties of Laplace transform domain 
solution. Thus far, their applications are 
limited to all lumped-element and distributed 
RC networks, which can only handle on-chip 
interconnects. For off-chip delay estimation, 
we took the second approach by incorpo¬ 
rating transmission line elements. 

Extensive work has been published on radi¬ 
ation by a microstrip patch excited by a 
probe. Most of the published work did not 
rigorously model the correct current distribu¬ 
tion on the probe. In the case of electrically 
thin substrates, the probe current is usually 
approximated as uniform along the probe. 
This leads to acceptable numerical results for 
the computation of the radiation pattern by 
the patch and the mutual impedance 
between the probe and the patch. However, 
this approximation is not sufficiently accurate 


to solve for the current distribution on the 
patch or for the computation of the self 
impedance across the terminals of the probe. 
Furthermore, this approximation is not appro¬ 
priate in the case of electrically thick 
substrates. 

We rigorously analyzed the radiation problem 
of a circular patch which is center fed by a 
coaxial-line driven probe over a ground plane 
and situated in an arbitrary layered medium. 
The current distribution on both the patch 
and the probe is rigorously formulated using 
a planar stratified medium approach. A set of 
three coupled integral equations is derived 
which governs the axial current distribution 
on the probe, the radial current distribution 
on the patch and the azimuthal magnetic 
current sheet across the aperture of the 
driving coaxial line. This set of equations is 
then solved using the method of moments. 
The resulting matrix equation is obtained in 
terms of Sommerfeld-type integrals that take 
into account the effect of the layered 
medium. These integrals are efficiently com¬ 
puted by a simple deformation in the 
complex wavenumber domain. The probe 
current distribution, input impedance and 
radiation pattern are presented and compared 
to the case of an uniform probe current dis¬ 
tribution. 

The analysis of microstrip antennas with 
electrically thick substrates has many appli¬ 
cations pertaining to millimeter wave systems 
and to achieving wide bandwidth. Most of 
the published work on the input impedance 
and other parameters of a probe-fed micro¬ 
strip antenna employ an approximation to the 
probe feed where it is replaced by a uniform 
current ribbon of equivalent dimensions. 
This approximation fails to give sufficiently 
accurate results for the input impedance of 
the probe especially for thick substrates. 
Recently, we have made some progress in 
this area — we obtained the current distribu¬ 
tion on the metallic probe in a closed mag¬ 
netic wall cavity. Then, we accounted for 
the radiation losses by lumping them arti¬ 
ficially into an effective loss tangent. 

We have formulated the problem of a center 
fed circular microstrip antenna in terms of a 
Weber transform, which allowed us to 
develop the Green’s function of the layered 
medium with the probe and the microstrip 
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patch as part of the medium. Using the 
Weber transform automatically enforces the 
boundary conditions on the probe and the 
patch. This allows us to cast the problem as 
the solution of a set of two coupled integral 
equations governing the tangential compo¬ 
nent of the electric field across the aperture 
of the coaxial line feed and that aperture 
across the interface where the patch lies. 
Next, we solve this set of equations using the 
method of moments. We then compute the 
current distribution on both the probe and 
the patch from the component of the mag¬ 
netic field tangential to their surfaces. Fur¬ 
thermore, from the computed electric field 
across the aperture of the coaxial line feed, 
we obtain the reflection coefficient for the 
TEM mode which allows us to compute the 
input impedance across the terminals of the 
probe. The probe current distribution, input 
impedance and radiation pattern are pre¬ 
sented, and the obtained results compared 
with those using the stratified medium for¬ 
mulation. 

The transient fields of a current source on a 
layered medium are calculated using the 
double deformation technique, in which 
complex integrals are deformed in the trans¬ 
verse wavenumber and frequency planes. 
Singularities from these complex planes cor¬ 
respond to physical modes of the structure, 
such as guided and leaky waves, and the rel¬ 
ative importance of each to the overall 
response can be discovered. Unlike the 
Cagniard-de Hoop method, double deforma¬ 
tion can be applied to dispersive and 
dissipative media. Also, the causality of the 
electromagnetic signal can be shown analyt¬ 
ically. 

We developed a modification to the double 
deformation technique, which splits the 
Fourier transform of the source current into 
two halves, one for times before the arrival at 
the observation point, and one after. This 
greatly increases the range of sources to 
which the double deformation technique can 
be appliel Another advantage of the mod¬ 
ification is the individual causality and conti¬ 
nuity of each mode. 

We have computed results for both line and 
strip currents on the surface of a coated 
perfect conductor for cases where the dielec¬ 
tric coating is both lossless and dissipative. 


In most cases, only a small number of modes 
suffices to reproduce the important features 
of the response, including the arrivals of 
reflected and lateral rays. The importance of 
each type of arrival depends on certain fea¬ 
tures of the time function, especially the 
initial slope. The response due to a strip 
current resembles that of a line current 
although there is some smoothing of the 
sharper features. 

The complex resonant frequencies of the 
open structure of a microstrip antenna con¬ 
sisting of two circular microstrip disks in a 
stacked configuration have been rigorously 
calculated as a function of the layered 
parameters and the ratio of the radii of the 
two disks. Using a dyadic Green's function 
formulation for horizontally stratified media 
and the vector Hankel transform, the mixed 
boundary value problem is reduced to a set 
of coupled vector integral equations. 
Employing Galerkin's method in the spectral 
domain, we calculated the complex resonant 
frequencies and demonstrated convergence 
of the results. We show that for each mode, 
the stacked circular microstrip structure has 
dual resonant frequencies which are associ¬ 
ated with the two coupled constitutive 
resonators of the structure and a function of 
the mutual coupling between the two disk 
resonators. This mutual coupling is a func¬ 
tion of the geometrical configuration of the 
stacked structure; the layered parameters, 
permittivities, permeabilities, and heights; and 
the ratio of the radii of the two disks. The 
dual frequency behavior of the stacked 
microstrip structure can be used to broaden 
the bandwidth or to provide for dual fre¬ 
quency use of the antenna. 
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We are developing mdEete of various kinds of 
earth terrain using J^two^ayer configuration 
to investigate tly polambetric scattering 
properties of renptely sensed media. The 
scattering layer/s an anisOTppic random 
medium characterized by a three-dimensional 
correlation fun/ion with lateralpid vertical 
correlation lengths and variances.^Based on 
the wave th^ry under Born appropriations, 
we applied/his model to derive t^e fully 
polarimetriopackscattering coefficients^ the 
Mueller anS covariance matrices. We cowd- 
ered a si«jle scattering process, taking iMo 
account#all the multiple reflections at tip 
boundafes. For an anisotropic randonr 
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Through our VLSI CAD research, we seek the 
means to produce custom integrated circuits 
correctly, quickly, and economically. Tradi¬ 
tionally, correctness has been checked at 
several representational levels such as layout 
(via design m!; checking), circuit, and logic 
(both via simulation). The techniques for 
checking correctness are local to the partic¬ 
ular representational level involved. While 
these techniques are important components 
of design testing, they do not attempt to 
provide alignment and consistency checks 
between the different levels of representation 
and an input behavioral specification. In 
addition, they do not characterize the set of 
possible designs at each representational 


level corresponding to the initial functional 
specification, while ranging over a variety of 
performance levels. For this reason, there is 
an increasing need to provide CAD tools that 
can serve as the framework for design explo¬ 
ration, providing the desired performance 
together with consistently aligned represent¬ 
ations at all levels. 

In this research group, we are studying a 
variety of research topics with an emphasis 
on performance directed synthesis of custom 
VLSI designs. Professor Allen has recently 
provided an overview of these issues, 1 in 
which he emphasizes the need for coordi¬ 
nating the design optimization process over 
the several levels of representation. 

Since there is so much to be gained by 
aggressive circuit design, our research group 
has placed considerable emphasis on this 
level. In addition, the group has focused on 
optimized architectures for digital signal pro¬ 
cessing (1) because of the highly parallel 
nature of the algorithms involved, and (2) 
because of the need for a very high level of 
real time performance. We have also 
increased our emphasis on developing means 
for formally specifying designs that facilitate 
performance directed design exploration, as 
well as efficient verification of correctness of 
the coordinated set of design level represent- 


1 J. Allen, 'Performance Directed Syn»hesis of VLSI Systems," Proc IEEE 78(2): 366-355 (1990). 
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ations. All of these studies are taking on an 
increasingly theoretical approach, signifying 
the transition of digital system design from 
an art to a science. 

The circuit level of design representation is 
the most abstract level at which performance 
can be characterized explicitly in terms of 
delay, and also where performance can be 
substantially manipulated through changes in 
circuit style and device sizing. It is an impor¬ 
tant level of design because optimization 
techniques which manipulate circuit repre¬ 
sentation can have great impact on resulting 
performance. More abstract representations 
can be transformed into a circuit represen¬ 
tation, but the need to transform a circuit 
into layout in a general yet efficient way is 
needed for the rapid design of custom cir¬ 
cuits. Baltus 2 has developed a very general 
yet efficient technique for compiling circuits 
into layout. The algorithm converts both 
NMOS and CMOS circuits in any style, and 
provides special means for the efficient 
layout of circuits with highly varying tran¬ 
sistor width-to-length ratios. While initially 
focused on implementation technologies 
with only one level of metal, it is now being 
generalized to multiple levels of metal. 
Recent comparisons of this program with 
other techniques at a benchmark competition 
have demonstrated the effectiveness of these 
procedures. 

Now that there is assurance of the generation 
of a high quality layout from any circuit 
specification, we are turning our attention to 
design performance exploration at the archi¬ 
tectural level. For this research, our goal is to 
study the performance of quasi-regular archi¬ 
tectures through manipulations of an initial 
behavioral specification. A new language has 
been developed for algorithms which can be 
compiled into either regular systolic architec¬ 
tures, or structures with a limited degree of 
irregularity. From the input language specifi¬ 
cation, a prototype descriptor is generated 
which compactly and efficiently represents 
both the structural hierarchy and the degree 


of replication of architectural elements. This 
descriptor can be readily expanded into a 
fully parallel data flow graph if necessary, but 
has the advantage that its complexity does 
not appreciably grow with circuit element 
replication. 

After obtaining the prototype descriptor, we 
apply a series of three types of transforma¬ 
tions in order to explore the architectural 
design space. First, various amounts of 
serialization can be instantiated, cutting 
down the amount of inherent parallelism in 
the input language specification, but also 
reducing the area of the resultant imple¬ 
mentation. Next, we explore a variety of 
architectural transformations that map the 
index spaces associated with each array vari¬ 
able into a common global index space. 
Finally, we explore pipelining to localize 
communication between architectural ele¬ 
ments. This approach is expected to lead to a 
flexible but powerful framework for architec¬ 
tural design exploration in the context of 
known capability for the conversion of the 
circuit form of the architectural elements into 
efficient layout. One important research issue 
is the degree of irregularity that is desirable 
in the resultant arrays due to non-local com¬ 
munication. It may be possible to make the 
distance of communication a designer- 
controlled parameter, thus providing a means 
to depart from regularity when the communi¬ 
cation needs of the algorithm favor it, 
thereby expanding the class of designs that 
can be efficiently converted to layout. 

Another example of the mapping of a high 
level specification to the circuit level repre¬ 
sentation has focused on the technology 
mapping process. Much research has focused 
on the conversion of a high-level functional 
specification to a highly optimized logic 
specification. This logic representation is 
usually then converted to the circuit form by 
means of a technology mapping in terms of a 
cell library. Heretofore the cell library was 
made up of typical restoring gates, such as 


2 D.G. Baltus and J. Allen, "SOLO: A Generator of Efficient Layouts from Optimized MOS Circuit Schematics,' 
Proceedings of the 25th Design Automation Conference. June 1988. 
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NAND and NOR cells, but Kaplan 3 has 
recently expanded this process to include 
pass transistor gates as well as restoring 
forms. By broadening the class of circuit 
forms, and hence basic library gate cells, new 
flexibility is introduced into the technology 
mapping process. New logic-level represent¬ 
ations have been introduced to facilitate the 
optimization of these circuits. By coupling to 
the existing SOLO layout generator, savings 
in the area of approximately 20 percent have 
been demonstrated for a number of designs 
that take advantage of the use of pass tran¬ 
sistor gates as well as restoring gates. This 
research exemplifies the need to coordinate 
design at several levels. Logic synthesis 
must contemplate the class of library cells 
that will be used in implementations, while 
increased circuit flexibility leads to new 
layout savings by means of a more powerful 
technology mapping. 

Driven by the requirements of performance 
optimization, we have devoted a considerable 
amount of research to delay characterization 
associated with circuit waveform transitions. 
McCormick 4 has developed an efficient 
scheme for computing these transitions accu¬ 
rately using a small number (four or less) of 
circuit transform moments. Interactions 
between interconnection lines, even as 
general as RLC lines, can be readily com¬ 
puted with these moments. In addition to 
passive interconnect, we have expanded 
these techniques to also include the efficient 
characterization of active circuits. Using 
these techniques, efficiency of delay calcula¬ 
tion, together with accuracy, is obtained, 
thus avoiding the substantial numerical load 
of comprehensive circuit simulation. Moment 
calculations are thus part of an overall 
strategy of performance estimation that must 
be done quickly and accurately in order to 


facilitate rapid exploration of design alterna¬ 
tives. It is also important to note that delay 
due to interconnect is now dominating the 
delay caused by active circuits in large 
systems. The ability to accurately charac¬ 
terize circuit delay in complex branching 
interconnect structures, including the effect 
of interconnect material together with driving 
and terminating impedances, is essential to 
the accurate exploration of circuit perform¬ 
ance. 

With the advent of very small devices oper¬ 
ating at very high speeds, we need to revise 
device models for circuit simulation. Rather 
than use average value parameters, we must 
now solve the device equations directly to 
provide the accuracy required. Unfortu¬ 
nately, direct use of these equations implies 
very long simulation run times, making it 
important to find ways of obtaining required 
accuracy with reduced computational 
expenditure. Reichelt 5 has applied waveform 
relaxation techniques to this problem in a 
way that readily allows for the use of parallel 
waveform relaxation simulators. Initial results 
have proved the effectiveness of this 
approach, which guarantees convergence of 
the simulation. Investigations of possible 
parallel architectures for this purpose are 
now under way. 

Recently, our emphasis on verification tech¬ 
niques has grown as the complexity of 
designs, and the need for overall correctness, 
has increased. There are a variety of 
approaches to verification, including logical 
proofs that implementation structures corre¬ 
spond to behavior specifications, and sym¬ 
bolic simulation. Because the need for 
formal characterizations of design has grown, 
we are searching for representations that can 
compactly characterize the class of correct 


3 J T. Kaplan, Pass and Restoring Logic Technology Mapping for the Efficient Automated Generation of Combina¬ 
tional Leaf Cells. S.M. thesis. Dept, of Electr. Eng. and Comput. Sci., MIT, 1989. 

4 S.P. McCormick and J. Allen, "Waveform Moment Methods for Improved Interconnection Analysis.' Proceedings 
of the 27th Design Automation Conference, Orlando, Florida, June 1990. 

8 M. Reichelt, J. White, and J. Allen, 'Waveform Relaxation for Transient Simulation of Two-Dimensional MOS 
Devices.' Proceedings of the International Conference on Computer Aided Design. November 1989. 
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designs. Bamji 6 has applied context free 
formal grammars to the specifications for the 
class of MOS digital integrated circuits. 
Compact grammars can represent classes of 
infinite cardinality, so that, for example, the 
class of all CMOS static combinatorial cir¬ 
cuits can he gei.erated by only five grammat¬ 
ical rules. Many circuit styles can be 
characterized in this way, and the corre¬ 
sponding parsing procedures can deal with 
mixed styles in an efficient way. Bamji has 
represented the class of all layouts generated 
by a regular grammar by use of the Regular 
Structure Generator, and he has shown how 
to verify the correctness of the resulting 
layout by parsing techniques. Finally, both 
of these grammatical representations, one at 
the circuit level and the other at the layout 
level, have been combined into a coupled 
grammar that can be used to verify the 
proper correspondence between layout and 
circuit. Thus, not only can both layout and 
circuit structures be verified independently by 
similar grammatical techniques, but we can 
show that their alignment is correct so that 
they are both coordinated representations of 
the same underlying abstract circuit in a con¬ 
sistent fashion. Successful verification dem¬ 
onstrates that all circuit and layout structures 
have correspondents at related layout and 
circuit representations, respectively, so that 
no loose fragments at either level are isolated 
from those at other design levels. Since these 
grammatical techniques have been successful 
at the circuit and layout levels, we are now 
focusing on their application at higher levels 
of design. 

Grammatical techniques are appropriate for 
characterizing correct structure of circuits, 
but they are not easily adapted to the verifi¬ 
cation of nonstructural circuit properties such 
as device sizing and allowable capacitance 
sizes that could lead to undesirable charge 
sharing. To deal with these attributes. Van 


Aelten 7 has developed an efficient technique 
for verification of circuit properties. This 
approach builds on the grammatical tech¬ 
niques of Bamji, parsing the circuit first to 
reveal its constituent structure. Semantic 
checks are then applied to the circuit to 
verify its specified performance. The circuit is 
effectively mapped into logic level represent¬ 
ations by a process of recursive interpreta¬ 
tion. Checks for electrical constraints are 
made at the level of the individual subcir¬ 
cuits, and the behavioral description of these 
elements is composed in order to provide a 
match with the user supplied specification. 
Extensions to this approach require more 
comprehensive circuit grammars that may 
have to be interleaved with the semantic 
interpretation process, thus breaking down 
the traditional separation of syntax and 
semantics. 

To continually maintain a consistent set of 
design representations that are properly 
aligned, it is important to build a unified 
framework for all levels of representation. 
Armstrong 8 has developed a technique for 
specifying design at all levels of abstraction 
with a uniform formalism. The way in which 
each of these abstractions is a consistent 
projection of one underlying abstract design 
entity is explicitly described, thus facilitating 
the automatic alignment of these represent¬ 
ations. In this way, when one design 
projection at a specific level of representation 
(e.g., circuit level) is manipulated, all other 
levels are changed automatically to provide 
the correct alignment. In general, this is a 
one-many process, since there is usually an 
infinite number of less abstract design 
projections corresponding to a more abstract 
projection of the same overall design. All that 
is required, however, is that all levels are 
consistent, so that performance optimization 
algorithms can select the best design 
projection at a given level, subject to some 


« C.S. Bamji and J. Allen, 'GRASP: A Grammar-based Schematic Parser." Proceedings of the 26th Design Auto¬ 
mation Conference, June 1989. 

7 F. Van Aelten and J. Allen, "Efficient Verification of VLSI Circuits Based on Syntax and Denotational Semantics," 
Proceedings. IMEC/IFIP Workshop on Applied Formal Methods for Correct VLSI Design, Houthalen, Belgium, 
1989. 

8 R.C. Armstrong, A Formal Approach to Incremental Consistency Maintainence in Multirepresentation VLSI Data¬ 
bases, Ph.D. diss. proposal, Dept, of Electr. Eng. and Comput. Sci., MIT, 1987. 
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overall optimization criterion. Upon com¬ 
pletion, this new formal framework will 
provide a basis for interconnecting an entire 
suite of design tools and will easily contem¬ 
plate the addition of new levels of represen¬ 
tation without the need to provide an entirely 
new framework. 

The above descriptions show how several 
aspects of performance-driven synthesis are 
being studied and adapted to new CAD 
tools. There is a clear trend toward adopting 
more formal notions of behavior, structure, 
and attributes, together with techniques for 
coordinating these several levels of 
description in one overall consistently main¬ 
tained structure. In the future, we will 
emphasize completion of the theoretical 
framework and unification of optimization 
procedures across all levels of design. 


The Vision Chip Project 

kors 

DuPonm^orporation 
NationalEcience Foundation 
GrantVlP 88-14612 

Project St* 

Professor Joffb L. Wyatt, Jr., Professor Berthold 
K.P Horn, Professor Hae-Seuy Lee, Professor 
Tomaso Poggi<\ Professor Charles G. Sodini, 
Steven J. Deckit Ibrahim th. Elfadel, Mikko 
Hakkarainen. CraicjVeast. Igdacio McQuirk, Mark 
N. Seidel, David \ Starkey, Christopher B. 
Umminger, Woodwar*Yanc 


6.2.1 Overall Proifeatt Description 
Problem Definition anti Methods 


Computational Cfemands of'Vision and the 
Smart Sensor Paradigm: A ^ajor problem 
in machine vis/n is the extraordinary volume 
of input data/at must be acquir^L managed 
and processes. This large quantity^ data in 
real-time cu<6y-level images leads to\ommu- 
nication 7 bottlenecks between \nager, 
memor^^nd processors, while the computa- 
tionaWrlemands remain high once data have 
reached the processors. The result is that 


conventional machines are incapable of any 
Jt the most rudimentary forms of real-time 
v*ion processing. m 

TheL smart rensor paradigm shov^E great 
proirose as a new approach to thisproblem. 
Alth<\gh potentially applicable to/a variety 
of sensor types, e.g., tactile an/ acoustic 
arrays, \n this project we are restricting its 
application to vision. The key/idea is to 
incorporate signal processing as early as pos¬ 
sible into\a system's signal flpw path to 
reduce deniands for transmission bandwidth 
and subsequent computation. / 

t 

In this paradigm, the major thepies are: 

processing/circuitry inte- 
^ingle chip, j 

• parallel comcititation, / 

p 'ocessing circuits distrjbuted throughout 


sensors ar 
irated on a^ 



array for s 
with minim 


jtially peal operations in 
wiring overhead, 

4 

^of analog ciriWsyfor early processing 
to avoid the speed and area costs 
ialog-to-digiteU/conversion on high- 
L /idth data straams, 

r 

• select^pn of tafekfc and algorithms 
requiring low to moderate precision, and 

omputations that 
cal processes in 
processes or to 


special;emphas& on 
map naturally bo phy^ 
silicon, e.g., tq^relaxatit 
resistive grids.! 


The Gi 
matic 
Requiri 
aithou 
achie 



Analog Intei 

continuous-tim 
integrated circ 
research. An 
advantages in 
silicon area, a 
ently analog 


ted Circuits: Both 

and sampled-data analog 
ts play a keyi role in our 
og circuits offer substantial 
peed, functionality per unit 
ease of interfacing to inher- 
put data, while digital circuits 
unlimited 
jly longer 


are easier to Resign and can obtair 
precision through use of progressiv 
wordlengths 


2-quadrant multiplier i\ a dra- 
of analog capability, 
only four transistors, its similar, 
slightly more complex circuits, 
).02 percent precision multiplication 


in unoer 200 nanoseconds with a dynamic 
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SUBJECT: Hold and Denial Records for FDAB FL-88 Orders(e.g. H-0475, D-0476) 


PROBLEM: FDAB regularly generates FL-88 orders, resulting from user 

requests, work unit summaries, bibliographies, etc. On occasion, FDAC 
must hold one of our orders or return the document to the origi’«ntor, 
resulting, in a denial. Reasons for holds include lack of distribution 
statement, missing pages, improper classification markings, reproducibility, 
illegibility. If answers to the above reasons are not received, denials 
are the result. There is a lack of consistent communication between 
FDAB and FDAC regarding FDAB orders that are held or denied. For example, 
FDAB is iml always informed if a hold or denial is taking place. Also, 
the AQ inputter in FDAC inputs her own record instead of using the 
AO record established by FDAB. 

PROPOSED SOLUTION: FDAB is requesting more control over holds and denials 
on FL-88 orders. The proposed solution is for FDAB to input the hold/denial 
information into the AQ database for FL-88 orders. With denials, FDAB 
will be able to close out the files and inform users, when appropriate. 

With holds, FDAB records will be more up-to-date. It is imperative that 
iDAB be inlormed of the outcomes on holds. 


Sharon Serzan 
16 Oct 89 


Discussed with C. Gibson, Acting Chief of FDAC, and we will proceed with 
above procedures. 


Sharon Serzan 
23 Oct 89 






DEPARTMENT OF THE AIR FORCE 

USAF ENVIRONMENTAL TECHNICAL APPLICATIONS CENTER (MAC) 
SCOTT AIR FORCE BASE, ILLINOIS 62225-5458 


FROM: LDD (DSN 576-2625) 19 OCT 90 

SUBJECT: Change In document status 

TO: DTIC-FDAC 

1. The distribution statement for four technical reports written by Air 
Weather Service personnel was changed to "Approved for Public Release; 
Distribution Unlimited." Photocopies of the title page, Review and Approval 
Statement, and Report Documentation Page for each are attached (Atch 1-4). 

2. Please delete 7WW/TN-82/002, AD A 123352 from the DTIC database. It has 
been declared obsolete. The same photocopied pages are attached (Atch 5). 


SUSAN A. TARBELL 5 Atch 

Cataloging Librarian 1. AD B 099412L 

Air Weather Service Technical Library 2. AD B 130894 

3. AD B 066230L 

4. AD B 082564L 

5. AD A 123353 












Table of Fnntvalents i£HL Outcnlr.t.i r? Cubic Footage of Records 


stter size material and any material 
filed in letter size folders . . . . 


1 full file drawer = 1.5 cubic feet 

or 

1 linear foot = .75 cubic foot 


Legal size material and any material 
filed in legal size folders. . . . 


1 full file drawer = 2 cub.z feet 

or 

1 linear foot = 1 cubic foot 


Cards in vertical files: 

3 by 5 inch cards. 

4 by 6 inch & tabulating size curds. 


1 linear foot = 
1 linear foot = 


. 1 cubic foot 
.17 cubic foot 


5 by 8 inch cards.. . 1 linear foot = .25 cubic foot 

■Cards in visible index trays: 

3 by 5 inch cards.1,000 full pockets = ..C/7 cubic foot 


4 by 6 inch 4 tabulating size cards. . '•, CCO full pockets = .11 cubic foot 

5 by S inch cards.1,CG0 full pockets = .18 cubic foot 

S by 10 1/2 inch cards. 1,000 full pockets = .36 cubic foot 


EAM listings, 14 by 17 inches or larger. 

HIM listings, small sizes. . .. . 

Records Shipment and Storage Boxes 

(15 by 12 by 10 inches). 

Maps and Drawings (Mixed sizes). 

(Since a wide range of sizes are 
frequently interfiled, using these 
average conversion factors is suf¬ 
ficiently accurate. In the case 


1 linear foot =1.5 cubic feet 
I linear foot = 1 cubic foot 


i box = 1 cubic foot - ' 

a 1-inch"stuck = .5 cubic foot 
or 

1,000 documents = 3 cubic feet 


menta are mostly a single size, a 
more accurate conversion factor may 
be determined ar.u used. ) 


ADDS mavn- 


;upe 


10 reels - 1 cubic foot 


Film; 

16 millimeter, ICO 1 reel 
1o millimeter, 10C0 1 red 
35 millimeter 
35 millimeter, 1G00 1 reel 


.... ICO reels = 1 cubic foot 
..... 12 reels = 1 cubic foot 

ICO’ reel.50 reels = 1 cubic foot 

.• 6 reels = 1 cubic foot 


Other records 


Intimate dimensions and compute cubic 
measurement 






















DEPARTMENT OF THE ARMY 

THE JUDGE ADVOCATE GENERAL’S SCHOOL 
CHARLOTTESVILLE. VIRGINIA 22903-1781 



JAGS-DD (310) 


4 September 1990 


MEMORANDUM FOR Defense Logistics Agency, Defense Technical 

Information Center, ATTN: Mr. Norm Walton, 
Cameron Station, Building 5, Alexandria, VA 
22304-6145 

SUBJECT: Request for Inclusion in and Deletion From DTIC 

1. Please add the following items to DTIC: 

Legal Assistance Guide: Soldiers' and Sailors' Civil Relief 
Act, JA-260-90 

Model Tax Assistance Program, JA 275-90 

2. Please delete the following publications from DTIC: 

AD B142445, Legal Assistance Guide: Soldiers' and Sailors' 
Civil Relief Act, JA-260-90, dated January 1990 

AD B124120, Model Tax Assistance Program/JAGS-ADA-88-2 

3. Point of contact for this action is Ms. Eva Skinner, <V T 274- 
7110, extension 972-6396; commercial (804) 972-6396. 



Enclosures 


MATTHEW E. WINTER 
CPT, JA 

Publications Officer 







